数值分析代码

《数值分析代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、拉格朗口插值#include#includeusingnamespacestd;intlagrange(double*px,double*py,double*pcoeff,intiNum){double*pk,*ptemp_coeff;inti,j,k,kk;pk二newdouble[iNum];ptemp_coeff二newdouble[iNum];for(i=0;i

2、um;j++){if(i!=i){pk[i]*=px[i]-px[j];}}}for(i=0;i0;kk—){ptemp_coeff[kk]-=ptemp_coeff[kk-l]*px[j];}}}for(j=0;j

3、elagcaculate(double*pcoeff,doubledxx,intiNum){doubledsum;inti;dsum=0.0;for(i=0;i>iNum;cout«n是否进行自动测试(测试函数为正弦函数)，是输入1,否输入0：“；cin>>iisrandom;px

4、二newdouble[iNum];py=newdouble[iNum];pcoeff二newdouble[iNum];if(iisrandom二二0){cout«n请输入iNum个已知点的横坐标:"«endl;for(i=0;i>px[i];}cout«"请输入iNum个己知点的纵坐标:"«endl;for(i=0;j>py[i];}else{for(i=0;i

5、"输入测试值数量”；cin»testnum;double*pxvalue;pxvalue二newdouble[testnum];double*pyvalue;pyvalue二newdouble[testnum];cout«"输入"«testnum«"个需要测试的值u«endl;for(i=0;i>pxvalue[i];}for(i=0;i

6、'“VV"标准函数值为“<#includevoidinain(){intn,ij；doubleA[50][50];doublex[50J,y[50J;doubleK=1玉=aN=0,P;coutvv"请输入所求均差阶数:”；cin»n;fbr(i=O;i<=n;i-H-){coutvv“潇俞入x”vviv

7、v“=”；cin»x[i];coutvv“请输入y”vvivv”=”；cin»y[i];A[i][0]=x[i];A[i][l]=y[i];}for(j=2;j<=n4-lj++){fbr(i=l;i<=n;i-H-){})ibr(i=0;iv二n;i卄){coutvv“输出第Wvivv邛介均差为H«A[i][i+l]«cndl;}coutvv“请所要代入计算的x的值：X』；cin»X;for(i=0;i