21、多项式插值，最小二乘法和牛顿迭代法

1、多项式插值
利用函数f (x)在某区间中已知的若干点的函数值，作出适当的特定函数，在区间的其他点上用这特定函数的值作为函数f (x)的近似值，这种方法称为插值法。如果这特定函数是多项式，就称它为插值多项式。

//*x、*fx分别为已知点数组（x0，f（x0））..（xn，f（xn））；n为已知点个数
// *z、*fz分别为待求点数组（z0，f（z0））..（zn，f（zn））；m为待求点个数
int interpol(const double *x, const double *fx, int n, double *z, double *pz, int m)
{
    double term, *table, *coeff;
    int i, j ,k;

    //table存临时差商，每组从一阶到n阶更新，对应上图的每行更新 ，
    //如 [x0 x1]，[x1x2],...，[xn-1 xn]；更新到[x0 x1 x2],...,[xn-2, xn-1, xn] 
    //coeff存最终差商 ，机每次取table更新的第一项，如[x0 x1]到[x0 x1...xn] 
    if((table = (double *)malloc(sizeof(double) * n)) ==