テイラー展開のプログラムについて
テイラー展開のプログラムについて
テイラー展開のプログラムについて
プログラミングはまだ始めたばかりですが、宿題がいきなり出されました。
指数関数をテイラー展開する際に次のプログラムがありましたが、これを表している意味がよく分かりません。
(1)EPSとは??(2)Tnnは何?(3)format('EXP(', f10.5, ') ')の「10.5」は何か?
ここら辺が特によく分かりません。
また、同じようにsinxや対数関数でテーラ展開のプログラムのとき気をつけなくちゃいけない点なども教えていただければ嬉しいです。
subroutine expotn (x)
c
real Tn, ans, err
integer n
c
EPS = 1.0e-6
n = 0
Tn = 1.0
ans = Tn
write(6,10) x
10 format('EXP(', f10.5, ') ')
50 if( abs(Tn) .gt. EPS) then
n = n + 1
Tn = Tn * x / n
ans = ans + Tn
err = ans - exp(x)
write(6,11) n, ans, err
11 format(I3,' exp(X) = ',f12.7,' err = ',f15.10)
go to 50
end if
write(6,*) ' Exp( ',x,' )= ',ans
return
end
お礼
ご回答ありがとうございます。 URL、とても参考になります!頑張ってみます。