!STH NBU

!FOLLOWING PROGRAMME IS IN FORTRAN 77 (I DID SOME MISTAKE WHICH I CAN'T SLOVE KINDLY ANYONE SOLVE  PLZ... INBOX HERE 

!'***NEWTON'S FORWARD interpolation****'

       REAL f(0:5),Y(0:5)

                    Print*,'***NEWTON'S FORWARD interpolation****'

       x0=1

       X=1.5

          h=0.2

                n=5

          data (f(i),i=0,5)/2.317,2.425,2.522,2.609,2.689,2.762/

       p=(x-x0)/h

​

        do i=0,(n-1)

          diff(i)=(f(i+1)-f(i))

       enddo

       coeff=p

          yx=f(0)+coeff*diff(0)  

          do i=2,n

       coeff(i)=coeff(i)*(p-i+1)

       ENDDO

       do i=0,n-1

          diff(i)=diff(i+1)-diff(i)

          yx(i)=yx(i)+coeff(i)*diff(0)

       enddo 

​

       Print*,yx

       End