Тексты программ:
program lab5
real x,An,eps,sum,xc,f,T
integer n,Nmax/100/
open(1, file='An.txt')
open(2, file='sum.txt')
open(3, file='f.txt')
open(4, file='error.txt')
write(*,*) 'input x, eps and press Enter:'
read(*,*) x,eps
f=asin(1.)-atan(x)
xc=-x**2.
n=0
An=-x
sum=asin(1.)-x
do while(abs(An)>eps)
T=((2.*n+1)/(2.*n+3))*xc
An=An*T
sum=sum+An
n=n+1
if(n>Nmax) then
write(4,*) 'Ошибка при x=',x,' n=',n,' sum=',sum,' f=',f
endif
write(1,*) n,An
write(2,*) n,sum
write(3,*) n,f
enddo
end
program lab5
real x,Xn,Xk,step,An,eps,sum,xc,f,T
integer n,Nmax/100/
open(1, file='results.txt')
write(*,*) 'input Xn, Xk, step, eps and press Enter:'
read(*,*) Xn,Xk,step,eps
write(1,*) 'Xn=',Xn,' Xk=',Xk,' step=',step,' eps=',eps
write(1,10)
do x=Xn,Xk,step
f=asin(1.)-atan(x)
xc=-x**2.
n=0
An=-x
sum=asin(1.)-x
do while(abs(An)>eps)
T=((2.*n+1)/(2.*n+3))*xc
An=An*T
sum=sum+An
n=n+1
if(n>Nmax) then
write(1,*) 'Ошибка при x=',x,' n=',n,' sum=',sum,' f=',f
endif
enddo
write(1,12) x,f,sum,n
enddo
write(1,11)
10 Format(&
' ___________________________________________ '/&
' | x | станд |Сумма ряда |Членов ряда|'/&
' |-------------------------------------------|')
11 Format(&
' |______|____________|___________|___________|')
12 Format(1x,'|',F5.2,2(' |',F11.7),'|',I6,5x,'|')
end
Графики:
Зависимость величины суммы ряда Тейлора от номера члена ряда - sum(n)
Значение функции f(x)
Зависимость значения члена ряда от его номера – An(n)
Вывод рекуррентной формулы
Санкт-Петербургский Государственный Университет
Информационных технологий,
Механики и оптики
Отчет
К лабораторной работе по информатике
Вариант № 8
Студент:
Преподаватель:
2005 г.
Результаты:
Xn= -9.000000E-01 Xk= 9.000000E-01 step= 1.000000E-01 eps= 1.000000E-07
___________________________________________
| x | станд |Сумма ряда |Членов ряда|
|-------------------------------------------|
| -.90 | 2.3036120 | 2.3036110| 54 |
| -.80 | 2.2455370 | 2.2455380| 27 |
| -.70 | 2.1815220 | 2.1815220| 18 |
| -.60 | 2.1112160 | 2.1112160| 13 |
| -.50 | 2.0344440 | 2.0344440| 10 |
| -.40 | 1.9513030 | 1.9513030| 7 |
| -.30 | 1.8622530 | 1.8622530| 6 |
| -.20 | 1.7681920 | 1.7681920| 4 |
| -.10 | 1.6704650 | 1.6704650| 3 |
| .00 | 1.5707960 | 1.5707960| 0 |
| .10 | 1.4711280 | 1.4711280| 3 |
| .20 | 1.3734010 | 1.3734010| 4 |
| .30 | 1.2793390 | 1.2793400| 6 |
| .40 | 1.1902900 | 1.1902900| 7 |
| .50 | 1.1071490 | 1.1071490| 10 |
| .60 | 1.0303770 | 1.0303770| 13 |
| .70 | .9600703 | .9600703| 18 |
| .80 | .8960553 | .8960553| 27 |
|______|____________|___________|___________|