収束加速法：Aitken 法

$n$ 個の数列に対する Aitken 加速のサブルーチンを使って、次の ような数列を加速してみる。

c     MAIN ROUTINE:EXAMPLE
c
      implicit real*8 (a-h,o-z)
      dimension a(20),b(20),c(20)
      integer i,n
c
      n=10
      do i=1,n
         b(i)=0.0d0
         c(i)=0.0d0
      end do
      do i=1,n
         a(i)= 0.7 + 0.5*(0.9)**i + 0.2*(0.38)**i
      end do
      call Aitken(a,b,n)
      call Aitken(b,c,n-2)
      do i=1,n
         write(*,1000) i,a(i),b(i),c(i)
      end do
 1000 format(i4,3e18.10)
      end

これは、2段の Aitken を加速した例である(入出力の format 文の書式も調べるこ と)。