[問題] 泰勒級數求pi

[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

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.115.36.238
推 jubilee2:把 b=float(i) 放到迴圈裏面
→ YCTzeng:應該不能編譯吧? implicit none了，可是c=-1？c沒定義?