作者jubilee2 (3321)
看板Fortran
标题Re: [问题] 泰勒级数求pi
时间Wed Sep 28 13:48:17 2011
※ 引述《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
1F：推 FlandreUFF:我算出来了，不过要到10^9个点数字才比较好看 09/28 14:51