用牛顿插值第一步为求差商：
我们把 0 阶差商到 n 阶差商都放在一个一维数组里面就会发现其中的规律，也可以根据这些规律来指导我们编程。
0 阶差商 n 个，一阶差商 n-1个 .... 直到 n 阶差商 1个。
求差商代码为：
bool CShuzhijisuanDlg::Diff(){ int cur; cur=m_length*(m_length+1)/2+1; m_diff=new double[cur]; for(int i=0;i<m_length;i++){//先存放0阶差商 m_diff[i]=parmy[i]; m_diff[i]; } //存放二阶以上差商 cur=m_length; int temp; for(int internal=1;internal<m_length;internal++){ for(i=internal;i<m_length;i++){ temp=cur-m_length+internal-1; m_diff[cur++]=(m_diff[temp]-m_diff[temp+1])/(parmx[i-internal]-parmx[i]); } } m_isdiff=true; return true;}差商求完我们开始求牛顿插值了：
double CShuzhijisuanDlg::Newton(double xx){ if(m_isdiff==false) Diff();
short int cur=m_length; double y,z=parmy[0]; for(int i=1;i<m_length;i++){ y=1; for(int j=0;j<i;j++){ y=y*(xx-parmx[j]); } z+=m_diff[cur]*y; cur=cur+m_length-i; } return z;}
