推 FlandreUFF:我算出來了，不過要到10^9個點數字才比較好看 09/28 14:51
※ 引述《FlandreUFF (フランちゃんうふふ)》之銘言：
: 自己寫了一個泰勒級數展開式
: 如果正確的話應該是逼近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
program Console5
implicit none
integer i
real*8 a,b,c,sum
sum=0.d0
c=-1.0d0
do i=1, 100
c=c*-1
b=float(2*i-1)
sum=sum+4.0*(c/b)
end do
write(*,*) sum
end
我改成這樣子
不過算出來是3.13159
那個你原本的話最好都加上小數點
應該就會好一點
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.77.208