cos(x)+ i sin(x)=ρ_n

 投稿者:たろさ  投稿日:2017年 6月13日(火)11時07分37秒
 
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

 

戻る