西电计算方法第三次作业.doc

文档介绍：班级:021113学号:02111283姓名::x01245y01646880先求差商表:kf()f(,)f()f()f()000167--32246213488-88450于是x(x-1)(x-2)(x-4):(x)=△------------+=,t=()=+=,t=()=:已知函数表:=:#include<>floatLagrange(floatx[],floaty[],floatxx,intn){ inti,j; float*a,yy=0; a=newfloat[n]; for(i=0;i<=n-1;i++) { a[i]=y[i]; for(j=0;j<=n-1;j++) if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } deletea; returnyy;}voidmain(){ floatx[4]={,,,}; floaty[4]={,,,}; floatxx=,yy; yy=Lagrange(x,y,xx,4); printf("x=%f,y=%f\n",xx,yy);}运行结果:x==:()和Nn()。