These codes require an ASCII input file interp.dat of the following form:



N: Number of Polynomial Interpolation Points (Small)

First Sample (x1,y1)



Second Sample (x2,y2)



...



Nth Sample (xN,yN)



N1: Number of Error Evaluation Points (Large)



First Sample (x1,y1)



Second Sample (x2,y2)



...



N1th Sample (xN1,yN1) -These codes will produce two ASCII output files:

lagrange.out or newton.out. The first section contains the coefficients of the interpolation polynomial. The second section contains the interpolation polynomial evaluated at the interpolation points, including the error at each interpolation point. The third section contains the overall error of the interpolation polynomial as computed from the error evaluation points.

lagrange.plt or newton.plt. The interpolation polynomial evaluated at the error evaluation points.