!###########################################################
! page.8
! エイトケン加速のプログラム
!###########################################################
PROGRAM prog01
IMPLICIT REAL(8)(a-h,o-z)
IMPLICIT INTEGER(i-n)

INTEGER,PARAMETER :: nn=9,jj=nn/2-1+mod(nn,2)
REAL(8),DIMENSION(nn,0:jj) :: S

S=0.d0

n=1;j=0

S(1,0)=1.d0-1.d0/3.d0
WRITE(*,'("S(",i3,",",i3,")=",F10.7)') n,j,S(n,j) 
DO n=2,nn
S(n,0)=S(n-1,0)+(-1.d0)**n/(2.d0*dble(n)+1.d0)
WRITE(*,'("S(",i3,",",i3,")=",F10.7)') n,j,S(n,j) 
END DO

DO j=1,jj
WRITE(*,*) 
DO n=1+j,nn-j
S(n,j)=(S(n+1,j-1)*S(n-1,j-1)-S(n,j-1)**2)/(S(n+1,j-1)+S(n-1,j-1)-2.d0*S(n,j-1))
WRITE(*,'("S(",i3,",",i3,")=",F10.7)') n,j,S(n,j) 
END DO
END DO

STOP
END PROGRAM prog01
!###########################################################