作者FlandreUFF (フランちゃんうふふ)
看板Fortran
標題[問題] 泰勒級數求pi
時間Wed Sep 28 13:16:42 2011
自己寫了一個泰勒級數展開式
如果正確的話應該是逼近pi=3.14159.....
也就是pi=4*(1-(1/3)+(1/5)-(1/7)...)這樣下去
不過算出來反而是個某名其妙的負數
請問是哪裡出錯了呢？
程式碼如下：
implicit none
integer i
real a
real b
real sum
sum=0.
c=-1
b=float(i)
do i=1, 100
a=4*(((-1)**(b-1))/(2*b-1))
sum=sum+a
end do
write(*,*) sum
stop
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.36.238
1F：推 jubilee2:把 b=float(i) 放到迴圈裏面 09/28 13:56
2F：→ YCTzeng:應該不能編譯吧? implicit none了，可是c=-1？c沒定義? 09/29 18:03