《C实例编程(用牛顿插值第一步为求差商).docx》由会员分享,可在线阅读,更多相关《C实例编程(用牛顿插值第一步为求差商).docx(3页珍藏版)》请在金锄头文库上搜索。
1、 C实例编程(用牛顿插值第一步为求差商)0 阶差商 n 个,一阶差商 n-1个 . 直到 n 阶差商 1个。求差商代码为:bool CShuzhijisuanDlg:Diff()int cur;cur=m_length*(m_length+1)/2+1;m_diff=new doublecur;for(int i=0;im_length;i+)/先存放0阶差商m_diffi=parmyi;m_diffi;/存放二阶以上差商cur=m_length;int temp;for(int internal=1;internalm_length;internal+)for(i=internal;im_l
2、ength;i+)temp=cur-m_length+internal-1;m_diffcur+=(m_difftemp-m_difftemp+1)/(parmxi-internal-parmxi);m_isdiff=true;return true;差商求完我们开头求牛顿插值了:double CShuzhijisuanDlg:Newton(double xx)if(m_isdiff=false)Diff();short int cur=m_length;double y,z=parmy0;for(int i=1;im_length;i+)y=1;for(int j=0;ji;j+)y=y*(xx-parmxj);z+=m_diffcur*y;cur=cur+m_length-i;return z;