Re: [問題] 泰勒級數求pi

[latinboy (暱稱)]

※ 引述《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
推 madokayui:請問為什麼我貼上得到Error executing df.exe.