added (library) code for interpolation from http://orion.math.iastate.edu/burkardt/f_src/f_src.html