作者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