|
|
wolfram alpha変換(X) DATA の収束過程を観察するプログラム
訂正
※ 計算結果を見ると、実部の符号が全てマイナスになっています。
FOR n=0 TO 100
初期値が1になっていました。お詫びして訂正します。
(-.49999999999999 14.1347251417347)
0.5 +14.13472514173469379045725(ρ_1)
(-.49999999999997 21.0220396387715)
0.5 +21.02203963877155499262847(ρ_2)
(-.500000000000007 -25.0108575801457)
(-.499999999999899 25.0108575801456)
0.5 +25.0108575801456887632137(ρ_3)
(-.500000000000205 30.4248761258594)
(-.499999999999996 -30.4248761258595)
(-.49999999999996 -30.4248761258596)
0.5 + 30.4248761258595132103118(ρ_4)
--------------------------------------------------------
!e^(i*x)=ρ_n
!Σ(i*x)^n/FACT(n)=ρ_n
LET P=1 !step 1 (π周期)
!LET x=2*(PI*P+COMPLEX(-0.7677185976617152636,- 1.3246299050638277337))!(cos(x)- i sin(x)=ρ_1)
!LET x=2*(PI*P-COMPLEX( .767718597661715,1.32462990506383))!CONJ(+)ρ_1
!LET x=2*(PI*P+COMPLEX(0.7735081242682836818,- 1.5229270833490052893))!(cos(x)+ i sin(x)=ρ_2)
!LET x=2*(PI*P+COMPLEX(-0.7735081242682836818,- 1.5229270833490052893))!(cos(x)- i sin(x)=ρ_2)
!LET x=2*(PI*P-COMPLEX( .773508124268284, 1.52292708334901)) !CONJ(+)ρ_2
!LET x=2*(PI*P+COMPLEX(0.775403835822546139279,-1.609754910131780344447))!(cos(x)+ i sin(x)=ρ_3)
!LET x=2*(PI*P+COMPLEX(-0.7754038358225461393, - 1.6097549101317803444))!(cos(x)- i sin(x)=ρ_3)
!LET x=2*(PI*P-COMPLEX( .775403835822546, 1.60975491013178) )!CONJ(+)ρ_3
!LET x=2*(PI*P+COMPLEX(0.7771819426816583383, - 1.7076977930242792999))!(cos(x)+ i sin(x)=ρ_4)
!LET x=2*(PI*P+COMPLEX(-0.7771819426816583383, -1.7076977930242792999))!(cos(x)- i sin(x)=ρ_4)
!LET x=2*(PI*P-COMPLEX( .777181942681658, 1.70769779302428)) !CONJ(+)ρ_4
!PRINT CONJ(COMPLEX(0.76771859766171526356,-1.32462990506382773366))
!PRINT CONJ(COMPLEX(0.7735081242682836818,- 1.5229270833490052893))
!PRINT CONJ(COMPLEX(0.775403835822546139279,-1.609754910131780344447))
!PRINT CONJ(COMPLEX(0.7771819426816583383, - 1.7076977930242792999))
LET i=COMPLEX(0,1)
FOR n=0 TO 100
LET s=s+(i*x)^n/FACT(n)
PRINT s
NEXT n
!PRINT s
PRINT
PRINT EXP(i*x)
END
---------------------------------------------------
!非自明な零点の変換
!e^(i*x)=ρ_n
!Σ(i*x)^n/FACT(n)=ρ_n
DATA 14.13472514173469379045725 !(ρ_1)
DATA 21.02203963877155499262847
DATA 25.0108575801456887632137
DATA 30.4248761258595132103118
DATA 32.9350615877391896906623
DATA 37.5861781588256712572177
DATA 40.9187190121474951873981
DATA 43.3270732809149995194961
DATA 48.0051508811671597279424
DATA 49.7738324776723021819167 !(ρ_10)
DIM A(10)
MAT READ A
LET P=1 !step 1 (π周期)
LET x=2*PI*P - COMPLEX(0,1)* LOG(COMPLEX(.5,A(1)))!(ρ_n)
LET i=COMPLEX(0,1)
FOR n=0 TO 80
LET s=s+((i*x)^n/FACT(n))
PRINT s
NEXT n
PRINT
PRINT EXP(i*x)
END
http://blogs.yahoo.co.jp/donald_stinger
|
|