!###########################################################
! page.9
! リチャードソン加速のプログラム
!###########################################################
PROGRAM prog02
IMPLICIT REAL(8)(a-h,o-z)
IMPLICIT INTEGER(i-n)

INTEGER,PARAMETER :: nn=10,mm=nn
REAL(8),DIMENSION(nn,0:mm) :: S

S=0.d0

n=1;m=0

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

DO m=1,mm
WRITE(*,*) 
DO n=1,nn-m
S(n,m)=(dble((n+1)**m)*S(n+1,m-1)-dble(n**m)*S(n,m-1))/dble((n+1)**m-n**m)
WRITE(*,'("S(",i3,",",i3,")=",F12.7)') n,m,S(n,m) 
END DO
END DO

STOP
END PROGRAM prog02
!###########################################################