作者latinboy (昵称)
看板Fortran
标题Re: [问题] 泰勒级数求pi
时间Wed Sep 28 21:36:10 2011
※ 引述《jubilee2 (3321)》之铭言：
: ※ 引述《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
program oo
implicit none
integer i
real*8 sum = 0
do i = 1, 10000000, 2
sum = sum + 1.d0/(4*i*i-1)
end do
write(*,*) sum * 8
end
将这个公式的数学形式对fortran最佳化
提出可化简的加减乘除步骤
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 222.92.110.76
1F：推 madokayui:请问为什麽我贴上得到Error executing df.exe. 09/28 23:25