Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
41
Добавлен:
03.03.2015
Размер:
63.49 Кб
Скачать

Лабораторная работа №1 часть 1. Запись арифметических выражений на алгоритмическом языке. Вар31

  1. log(abs(x**2))-1/tan(x**3)-exp(x+z)

  2. (3*x*z**y-6*sin(x-2))/acos(2*x)

  3. ((x**4-1)**(1/5.)+sqrt(x+1))/(6*x*y*z)

  4. (1-x)/(x+3)-atan(abs(x))**2+6,9 часть 2. : Программирование

Программа на алгоритмическом языке:

print *,'vvedite s2,t'

read*,s2,t

x=s2-4*t

y=s2/t

Rez=(cos(x)**3+3*y)/(1+2*x+3*y)

print*,'Rez=',rez

pause

end

Результаты счета : режим отладки

vvedite s2,t

12 3

Rez= 1.000000

vvedite s2,t

-4.5 9.75

Rez= 7.8877676E-03

часть 3. : Вычисление корней квадратного уравнения

Program sq_ur

print*,'Vvedite znachenie parametra t'

read*,t

a = 2*t-3

b = -7*t

c = t-2

if(a.eq.0)then

x = -c/b

print*,'yravnenie imeet odin koren x=',x

pause

stop

endif

d = b**2 - 4*a*c

if(d.lt.0)then

print*,'Yravnenie deistvitelnih kornei ne imeet d<0'

else

x1 = (-b+sqrt(d))/(2*a)

x2 = (-b-sqrt(d))/(2*a)

print*,'korni yravneniya: x1=',x1,', x2=',x2

endif

pause

enD

! Rezyltati scheta:

! npu t=1.5 yravnenie imeet odin koren x=-4.7619049E-02

! npu t=-0.5 yravnenie deistvitelnih kornei ne imeet d<0

! npu t=2 korni yravneniya: x1=0.0000000E+00, x2=14.00000

Лабораторная работа №2. часть 1. Запись логических выражений на алгоритмическом языке

а) x.ge.-8.and.x.le.-2 b) (x.gt.-20.and.x.le.-1).or.( x.ge.2.and.x.lt.6)

часть2. Ветвящиеся алгоритмы.

Program labos_2_2

print*,'Vvedite x,y,z'

read*,x,y,z

if ((x.ge.1).and.(x.le.5)) then

print*,'znachenie x=',x

endif

if ((y.ge.1).and.(y.le.5)) then

print*,'znachenie y=',y

endif

if ((z.ge.1).and.(z.le.5)) then

print*,'znachenie z=',z

endif Вар31

pause

enD

! Rezyltati scheta:

Vvedite x,y,z

0 1 3

znachenie y= 1.000000

znachenie z= 3.000000

Лабораторная работа 3. Определение наибольшего и наименьшего значения функции и построение графика.

Program labos_3_table

f(x)=1*x**4+8*x**3-1*x**2+10*x+5

open(1, file='v1.dat')

print*,'Input a,b,n'

read*,a,b,n

xmin=a

xmax=a

fmin=f(a)

fmax=f(a)

do x=a,b,(b-a)/n

y=f(x)

if(y.gt.fmax)then

fmax=y

xmax=x

endif

if(y.lt.fmin)then

fmin=y

xmin=x

endif

write(1,*),x,y

enddo

print*,'naibolshee znachenie fynkcii =',fmax,' npu x=',xmax

print*,'naimenshee znachenie fynkcii =',fmin,' npu x=',xmin

pause

close(1)

enD

! Rezyltati scheta:

Input a,b,n

-5 5 20

naibolshee znachenie fynkcii = 1655.000 npu x= 5.000000

naimenshee znachenie fynkcii = -445.0000 npu x= -5.000000

Program labos_4_summa Вар31

open(1,file='sum.dat')

print*,'vvedite n'

read*,n

write(1,10)

10 format(5x,20('-')/5x,'|',4x,'x',4x,'|',4x,'S',4x,'|'/5x,20('-'))

do i=1,5

print*,'vvedite ocherednoe znachenie x'

read*,x

u=x**2/2

znak=-1

S=0

do k=1,n

S=S+znak*u

znak=-znak

u=u*x/2

enddo

write(1,20)x,S

20 format(5x,2('|',f9.4),'|'/5x,20('-'))

enddo

close(1)

enD

! Rezyltati shceta:

--------------------

| x | S |

--------------------

| 0.1000| -0.0048|

--------------------

| 0.3000| -0.0391|

--------------------

| 0.4000| -0.0667|

--------------------

| 0.7000| -0.1815|

--------------------

| 1.0000| -0.3330|

--------------------

Program laba_5_dimensions

dimension A(10),B(10),T(10)

print*,'Vvedite kolichestvo elementov massivov'

read*,m

print*,'Vvedite elementi massiva  A'

read*,(A(i), i=1,m)

print*,'Vvedite elementi massiva  B'

read*,(B(i), i=1,m)

do 1 i=1,m

if(B(i).gt.0)then

T(i)=A(i)/B(i)

else

T(i)=A(i)*B(i)

end if

1 continue

print 2,(A(i),i=1,m)

2 format (10x,'massiv A',/(10f8.2))

print 3,(B(i),i=1,m)

3 format (10x,'massiv B',/(10f8.2))

print 4,(T(i),i=1,m)

4 format (10x,'massiv T',/(10f8.2))

pause

enD

!rezyltati shcheta

Vvedite kolichestvo elementov massivov

5

Vvedite elementi massiva? A

6 4 8 -3 6

Vvedite elementi massiva? B

-2 -8 6 7 -3

massiv A

6.00 4.00 8.00 -3.00 6.00

massiv B

-2.00 -8.00 6.00 7.00 -3.00

massiv T

-12.00 -32.00 1.33 -0.43 -18.00

Лабораторная работа 6. Вычисление скалярного произведения. Алгоритм и ручной счет: Вар31

Program laba_6_Matrixs

dimension a(3,3),b(3,3),p(3),q(3),r(3),x(3),y(3),z(3)

open(1,file='lab6.dat')

read(1,*)((a(i,j),j=1,3),i=1,3),((b(i,j),j=1,3),i=1,3)

print 3,((a(i,j),j=1,3),i=1,3)

3 format (4x,'matrix A'/(3f4.1))

print 4,((b(i,j),j=1,3),i=1,3)

4 format (4x,'matrix B'/(3f4.1))

read(1,*),p,q,r

print 5,p,q,r

5 format(5x,'vectori'/2x,'P:',3f6.1/2x,'Q:',3f6.1/2x,'R:',3f6.1)

call mult(a,p,x,3)

call lin(1.,r,1.,q,y,3)

call mult(b,y,z,3)

print 6,scal(x,z,3)

6 format(3x,'Rezyltat s=',f8.3)

close(1)

pause

enD

Subroutine mult(a,b,c,n)

dimension a(n,n),b(n),c(n)

do 1 i=1,n

c(i)=0

do 1 j=1,n

1 c(i)=c(i)+a(i,j)*b(j)

enD

Subroutine lin(al,a,bt,b,c,n)

dimension a(n),b(n),c(n)

do 1 i=1,n

1 c(i)=al*a(i)+bt*b(i)

enD

Function scal(a,b,n)

dimension a(n),b(n)

do 1 i=1,n

1 scal=scal+a(i)*b(i)

enD

!Rezyltati sheta:

matrix A

1.0 2.0 3.0

1.0 2.0 1.0

3.0 2.0 0.0

matrix B

4.0 1.0 2.0

0.0 4.0 3.0

1.0 1.0 1.0

vectori

P: 0.1 1.7 -1.5

Q: -1.6 0.8 1.1

R: -0.7 1.3 0.2

Rezyltat s= 33.170

Соседние файлы в папке Информ. лабы комп