Информ. лабы комп / 20
.docЛабораторная работа №1 часть 1. Запись арифметических выражений на алгоритмическом языке. Вар20
-
(z**x)**(2*y)+(2*y)**(2*x)
-
3*cos(x/2)**2+sqrt(exp(-x)+2*z**2)
-
x**2/(2*3*4)+(x-4)**(1/7.)/(sin(x)**2+1)
-
x+(x**2+log(x)**2+0,3)/(x+x**2/(x+1))
часть 2. : Программирование
Программа на алгоритмическом языке:
print *,'vvedite x5,y'
read*,x5,y
u=3*x5-y
alpha=x5*atan(y)
Rez=((abs(u)+2*alpha)**(1/4.))/(cos(alpha)**4+3*exp(-alpha))
print*,'Rez=',rez
pause
end
Результаты счета:
vvedite x5,y
0 -16
Rez= 0.5000000
vvedite x5,y
-4.5 9.75
Rez= 7.9978682E-04
часть 3. : Вычисление корней квадратного уравнения
Program sq_ur
print*,'Vvedite znachenie parametra t'
read*,t
a = t-4
b = 4*t
c = t-5
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=4 yravnenie imeet odin koren x=6.2500000E-02
! npu t=0 yravnenie deistvitelnih kornei ne imeet d<0
! npu t=5 korni yravneniya: x1=0.0000000E+00, x2=-20.00000
Лабораторная работа №2. часть 1. Запись логических выражений на алгоритмическом языке
а) x.gt.-2.and.x.lt.5 b) (x.ge.-10.and.x.le.-3).or.( x.ge.9.and.x.lt.22)
часть2. Ветвящиеся алгоритмы.
Program labos_2_2
print*,'Vvedite x,y,z'
read*,x,y,z
if (abs(x).lt.3) then
print*,'znachenie x=',x**2
endif
if (abs(y).lt.3) then
print*,'znachenie y=',y**2
endif
if (abs(z).lt.3) then
print*,'znachenie z=',z**2 Вар20
endif
pause
enD
! Rezyltati scheta:
! Vvedite x,y,z
!2.5 4 0
! znachenie x= 6.250000
! znachenie z= 0.0000000E+00
Лабораторная работа 3. Определение наибольшего и наименьшего значения функции и построение графика.
Program labos_3_table
f(x)=2*x**4-39*x**3-0*x**2-70*x+4
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 = 6479.000 npu x= -5.0000000
naimenshee znachenie fynkcii = -3971.000 npu x= 5.000000
Program labos_4_summa
open(1,file='sum.dat')
print*,'vvedite n' Вар20
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
znak=1
S=0
do k=1,n
S=S+znak*u*k
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.0454|
--------------------
| 0.3000| 0.1134|
--------------------
| 0.4000| 0.1389|
--------------------
| 0.7000| 0.1920|
--------------------
| 1.0000| 0.2188|
--------------------
Program laba_5_dimensions
dimension P(16)
print*,'Vvedite m - kolichestvo elementov massiva P'
read*,m
print*,'Vvedite elementi massiva P'
read*,(p(i),i=1,m)
print 10,(p(i),i=1,m)
10 format (20x,'Massiv P:'/(10f8.2))
kp=0
km=0
kz=0
do 1 i=1,m
if(p(i).gt.0)then
kp=kp+1
end if
1 continue
do 2 i=1,m
if(p(i).lt.0)then
km=km+1
end if
2 continue
do 3 i=1,m
if(p(i).eq.0)then
kz=kz+1
end if
3 continue
print 20,'velichina p(i)>0:',kp
print 20,'velichina p(i)<0:',km
print 20,'velichina p(i)=0:',kz
pause
20 format(//2x,a,i2/2x,a,f10.3)
enD
!rezyltati shcheta
Vvedite m - kolichestvo elementov massiva P
5
Vvedite elementi massiva P
-2 3 7 0 1
Massiv P:
-2.00 3.00 7.00 0.00 1.00
velichina p(i)>0: 3
velichina p(i)<0: 1
velichina p(i)=0: 1
Лабораторная работа 6. Вычисление скалярного произведения. Алгоритм и ручной счет: Вар20
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 lin(1.,r,-1.,p,x,3)
call mult(b,x,y,3)
call mult(a,p,z,3)
print 6,scal(z,y,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= 9.050