!234567 implicit none external steffen integer i,nobs,nint parameter(nobs=10,nint=40) real xobs(nobs),yobs(nobs),xint(nint),yint(nint) do i=1,nobs xobs(i)=i-1 yobs(i)=xobs(i)**1.2 enddo do i=1,nint xint(i)=0.25*(i-1) enddo call steffen(nobs,xobs,yobs,nint,xint,yint) print*,' xint ',' analytical solution ', &' interpolated solution' do i=1,nint print*, xint(i),xint(i)**1.2,yint(i) enddo end