新しく発言する EXIT インデックスへ

　　εアルゴリズムによる無限級数の和の近似値 山中和義 2008/02/16 10:55:05 
　　└つづき 山中和義 2008/02/16 10:58:22 

Re: εアルゴリズムによる無限級数の和の近似値	返事を書く ノートメニュー
山中和義 <drdlxujciw> 2008/02/16 10:58:22
つづき !●収束が遅い場合 !log(x+1)=x-x^2/2+x^3/3-x^4/4+ …　※１ !!!LET n=10000 !第n項　※単純な加算は精度をあげる LET n=15 !第n項 MAT a=ZER FOR i=0 TO n !1,-1/2,1/3,-1/4,… LET a(i)=(-1)^i/(i+1) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT LOG(2) PRINT !※１より、1/2*log(2)=1/2*(1-1/2+1/3-1/4+ …)=1/2-1/4+1/6-1/8+ … LET n=15 !第n項 MAT a=ZER FOR i=1 TO n+1 !1/2,-1/4,1/6,-1/8,… LET a(i-1)=(-1)^(i+1)/(2*i) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT LOG(2)/2 PRINT !tan-1(x)=x-x^3/3+x^5/5-x^7/7+ … LET n=15 !第n項 MAT a=ZER FOR i=0 TO n !グレゴリ・ライプニッツ級数　1,-1/3,1/5,-1/7,… LET a(i)=(-1)^i/(2*i+1) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT PI/4 PRINT LET n=50 !第n項 MAT a=ZER FOR i=1 TO n+1 !1/(1*3),1/(5*7),1/(9*11),… LET a(i-1)=1/((4*i-3)*(4*i-1)) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT PI/8 PRINT LET n=200 !第n項 MAT a=ZER FOR i=1 TO n+1 !オイラー級数　1,1/4,1/9,1/16,1/25,…　自然数の2乗の逆数の和 LET a(i-1)=1/i^2 NEXT i PRINT EpsilonAlgorithm(a,n) PRINT PI^2/6 PRINT !●収束が速い場合 !e^x=1/0!+x/1!+x^2/2!+x^3/3!+ … LET n=15 !第n項 MAT a=ZER FOR i=0 TO n !自然対数の底　1/0!,1/1!,1/2!,1/3!,… LET a(i)=1/fact(i) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT EXP(1) PRINT !●精度が期待できない場合 LET n=20 !第n項 MAT a=ZER FOR i=0 TO n !(1/1-1/2)+(1/2-1/3)+(1/3-1/4)+… LET a(i)=1/(i+1)-1/(i+2) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT 1 PRINT LET n=200 !第n項 MAT a=ZER FOR i=1 TO n+1 !1/(1*2)+1/(2*3)+1/(3*4)+… LET a(i-1)=1/(i*(i+1)) NEXT i PRINT EpsilonAlgorithm(a,n) PRINT 1 PRINT END

インデックスへ

EXIT