fortran三次插值。reference:丁月蓉,天文数据处理(forinterpolationanddifferential);徐士良,FORTRAN常用算法程序集(forcalculus)