|
|
> No.1933[元記事へ]
南澤さんへのお返事です。
> シャープの級数によって円周率を求めたいです。
> 例えば小数点以下200までとしたとき、
> テイラー級数をどの程度求めれば良いのでしょうか？
!シャープの級数 π=6*ArcTan(1/√3)
OPTION ARITHMETIC DECIMAL_HIGH !1000桁モード
LET X=SQR(3)/3 !1/√3
LET N=414 !第n項まで
LET S=ArcTan(X,N)
PRINT "ArcTan(X)≒";S
PRINT "その6倍";6*S
PRINT USING "πとの差=-%.###^^^^^^":6*S-PI
END
EXTERNAL FUNCTION ArcTan(X,N) !ArcTan(x)の級数を第n項まで計算する
OPTION ARITHMETIC DECIMAL_HIGH !1000桁モード
LET F=-X*X
LET T=X
LET S=0
FOR J=1 TO N
LET S=S+T/(2*J-1)
LET T=F*T
NEXT J
LET ArcTan=S
END FUNCTION
|
|