c语言数值方法上机.doc

第一次作业 拉格朗日差值计算 #includestdio.h #includemath.h float sort(int n,float x,float (*p)[2]); main() { int i,j,n; float x,y; float (*p)[2]; scanf (%f,x); printf(x=%f\n,x); scanf (%d,n); printf(n=%d\n,n); for(i=0;in;i++) for(j=0;j=1;j++) { scanf(%f,*(p+i)+j); printf(x[%d][%d]=%f\n,i,j,*(*(p+i)+j)); } y=sort( n,x,p); printf(the result is %f\n,y); } float sort(int n,float x,float (*p)[2]) { int i,j; float a,b,k,c,e,d,y,q; b=k=1; y=0; for(i=0;in;i++) { e=*(*(p+i)+1); b=1; k=1; for(j=0;jn;j++) { if(j!=i) { a=*(*(p+i))-*(*(p+j)); b=a*b; c=x-*(*(p+j)); k=k*c; } } q=k/b; d=q*e; y=y+d; } return(y); } 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 0 0.0953 0.1823 0.2624 0.3365 0.4055 0.47 0.5306 0.5878 0.6419 0.6431 计算结果 把X=1.54代入的结果 把X=1.98代入得结果 牛顿向前插值 #includestdio.h #includemath.h float sort(int n,int i,float (*p)[2]); main() { int i,j,n,b; float x,N,k,y,m,y0,t,u; float h; float (*p)[2]; float a[100]={0}; float *p1; scanf (%f,h); printf (h=%f\n,h); scanf (%f,x); printf(x=%f\n,x); scanf (%d,n); printf(n=%d\n,n); for(i=0;in;i++) for(j=0;j=1;j++) { scanf(%f,*(p+i)+j); printf(x[%d][%d]=%f\n,i,j,*(*(p+i)+j)); } t=x-*(*p); t=t/h; y0=*(*(p+0)+1); for(j=0;jn-1;j++) { b=n; b=b-1; printf(\n); for(i=0;ib;i++) {y=sort( n,i,p); printf( %f,y); if(i==0) a[j]=y; } } i=0; j=1; m=t; N=y0; p1=a[0]; while(in-1) { u=*p1; N=N+m*u; printf (\n???%f???\n,N); j=j*(i+2); m=m*(t-1)/j; i++; p1++; } printf(%12f\n,N); } float sort(int n,int i,float (*p)[2]) { int j; float a,y; y=0; j=i++; a=*(*(p+i)+1)-*(*(p+j)+1); *(*(p+j)+1)=a; y=a; return(y); } 计算结果 当X=1.54时，f（x）=lnx的值 当X=1.98时，f（x）=lnx的值 牛顿向后插值 #includestdio.h #includemath.h float sort(int n,int i,float (*p)[2]); main() { int i,j,n,b,kk; float x,N,k,y,m,yn,t,u; float h; float (*p)[2]; float a[100]={0}; float *p1; scanf