[問題] 泰勒級數求pi
自己寫了一個泰勒級數展開式
如果正確的話應該是逼近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
推
09/28 13:56, , 1F
09/28 13:56, 1F
→
09/29 18:03, , 2F
09/29 18:03, 2F
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
