Re: [問題] 泰勒級數求pi
※ 引述《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
推
09/28 23:25, , 1F
09/28 23:25, 1F
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
