数值分析实验答案

《数值分析实验答案》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验0截断误差与舍入误差#include"stdio.h"#include"math.h"constdoubleln2=0.693147190546;constdoublee1=5e-6;voidmain(){intsign;doubles;longi;s=0.0;sign=1;i=1;while(fabs(ln2-s)>=e1){s+=(1.0/i)*sign;sign=-sign;i++;}printf("n=%ld",i-1);getch();}实验1拉格朗日插值法编写拉格朗日插值法通用子程序，并用以下函数表来上机求,。x0.00.10.1950.3

2、0.4010.5f(x)0.398940.396950.391420.381380.368120.35206#includemain(){staticfloatLx[10],Ly[10];intn,i,j;floatx,y,p;printf("entern=");scanf("%d",&n);printf("enterxi");for(i=0;i

3、scanf("%f",&x);/*n=6;Lx[0]=0;Lx[1]=0.1;Lx[2]=0.195;Lx[3]=0.3;Lx[4]=0.401;Lx[5]=0.5;Ly[0]=0.39894;Ly[1]=0.39695;Ly[2]=0.39142;Ly[3]=0.38138;Ly[4]=0.36812;Ly[5]=0.35206;x=0.15;*/for(i=0;i

4、tch();}实验2最小二乘法测得铜导线在温度(℃)时的电阻如下表，求电阻R与温度T的近似函数关系。i0123456(℃)19.125.030.136.040.045.150.076.3077.8079.2580.8082.3583.9085.10#include#includefloatgs(floata[20][20],floatb[20],intn){inti,j,k,l;floats;k=1;while(k!=n+1){if(a[k][k]!=0){for(i=k+1;i<=n+1;i++){a[i][k]=a[i][

5、k]/a[k][k];b[i]=b[i]-a[i][k]*b[k];for(j=k+1;j<=n+1;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j];}}k=k+1;}for(k=n+1;k>=1;k--){s=0;for(l=k+1;l<=n+1;l++)s=s+a[k][l]*b[l];b[k]=(b[k]-s)/a[k][k];}return0;}main(){staticfloatb[20],Lx[20],Ly[20],c[20][20],ct[20][20],a[20][20];intm,n,i,j,k=0,l;floats,r

6、tn;printf("enterm=");scanf("%d",&m);printf("entern=");scanf("%d",&n);printf("enterxi");for(i=1;i<=m;i++)scanf("%f",&Lx[i]);printf("enteryi");for(i=1;i<=m;i++)scanf("%f",&Ly[i]);for(i=1;i<=m;i++){c[i][1]=1;for(j=2;j<=n+1;j++)c[i][j]=Lx[i]*c[i][j-1];}for(i=1;i<=m;i++)for(j=1;j<=n+1

7、;j++)ct[j][i]=c[i][j];for(i=1;i<=n+1;i++)for(j=1;j<=n+1;j++)a[i][j]=0;for(i=1;i<=n+1;i++)b[i]=0;for(i=1;i<=n+1;i++)for(k=1;k<=n+1;k++)for(j=1;j<=m;j++)a[i][k]=a[i][k]+ct[i][j]*c[j][k];for(i=1;i<=n+1;i++)for(j=1;j<=m;j++)b[i]+=ct[i][j]*Ly[j];gs(a,b,n);printf("Theresultis:");for(i=1;

8、i<=n+1;i++){