计算机上机报告

《计算机上机报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、计算方法上机实验报告上课时间：2014-2015学年秋学期，6～14周一．拉格朗日插值------------------------------------------------------1二．牛顿插值------------------------------------------------------------3三．改进欧拉法---------------------------------------------------------5四．四阶龙格-库塔-----------------------------------------------------7五．牛

2、顿迭代------------------------------------------------------------9六．复化Simpson公式------------------------------------------------11七．Romberg算法------------------------------------------------------14八.Seidel迭代法------------------------------------------------------17九.Gauss列主元消去法----------------------

3、------------------------20一.拉格朗日插值1.程序代码#includevoidLagrange(){inti=0;doublea[10],b[10],L,L1,L2,L3,L4,x;cout<<"x=";for(i=0;i<4;i++){cin>>a[i];}cout<<"y=";for(i=0;i<4;i++){cin>>b[i];}cout<<"x=";cin>>x;L1=(x-a[1])*(x-a[2])*(x-a[3])*b[0]/(a[0]-a[1])/(a[0]-a[2])/(a[0]-a[3]);L2=(x-a[0])*(

4、x-a[2])*(x-a[3])*b[1]/(a[1]-a[0])/(a[1]-a[2])/(a[1]-a[3]);L3=(x-a[0])*(x-a[1])*(x-a[3])*b[2]/(a[2]-a[0])/(a[2]-a[1])/(a[2]-a[3]);L4=(x-a[0])*(x-a[1])*(x-a[2])*b[3]/(a[3]-a[0])/(a[3]-a[1])/(a[3]-a[2]);L=L1+L2+L3+L4;cout<<"L="<

5、stream.h>#includevoidmain(){intn,i,j;doubleA[50][50],*x,*y;cout<<"请输入插值节点数：";cin>>n;x=newdouble[n];y=newdouble[n];cout<<"请输入这"<>x[i]>>y[i];doubleK=1,xx,N=0,P;for(i=0;i<=n-1;i++){A[i][0]=x[i];A[i][1]=y[i];}for(j=2;j<=n;j++)for(i=1;i<=n-1;i

6、++){A[i][j]=(A[i][j-1]-A[i-1][j-1])/(A[i][0]-A[i-j+1][0]);}for(i=0;i<=n-1;i++)cout<<"输出第"<>xx;for(i=0;i#inclu

7、dedoublefun(doublex,doubley){return(-0.9*y/(1+2*x));}voidmain(){doublea,b,*y,h,*x,yp,yc;intn,k;cout<<"常微分方程为y'=-0.9*y/(1+2*x)"<>a>>b;co