10个重要的算法C语言实现源代码：拉格朗日,牛顿插值,高斯,龙贝格.docVIP

10个重要的算法C语言实现源代码：拉格朗日，牛顿插值，高斯，龙贝格~~ 关键字: 拉格朗日，牛顿插值，高斯，龙贝格 1.拉格朗日插值多项式 ，用于离散数据的拟合 C/C++ code #include stdio.h #include conio.h #include alloc.h float lagrange(float *x,float *y,float xx,int n)???? /*拉格朗日插值算法*/ { int i,j; ?? float *a,yy=0.0;??? /*a作为临时变量，记录拉格朗日插值多项式*/ ?? a=(float *)malloc(n*sizeof(float)); ?? 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]; ?? } free(a); return yy; } main() { int i,n; float x[20],y[20],xx,yy; printf(Input n:); scanf(%d,n); if(n=20) {printf(Error!The value of n must in (0,20).); getch();return 1;} if(n=0) {printf(Error! The value of n must in (0,20).); getch(); return 1;} for(i=0;i=n-1;i++) { printf(x[%d]:,i); ??? scanf(%f,x[i]); } printf(\n); for(i=0;i=n-1;i++) { printf(y[%d]:,i);scanf(%f,y[i]);} printf(\n); printf(Input xx:); scanf(%f,xx); yy=lagrange(x,y,xx,n); printf(x=%f,y=%f\n,xx,yy); getch(); } 2.牛顿插值多项式，用于离散数据的拟合 C/C++ code #include stdio.h #include conio.h #include alloc.h void difference(float *x,float *y,int n) { float *f; int k,i; f=(float *)malloc(n*sizeof(float)); for(k=1;k=n;k++) { f[0]=y[k]; ??? for(i=0;ik;i++) ????? f[i+1]=(f[i]-y[i])/(x[k]-x[i]); ??? y[k]=f[k]; } return; } main() { int i,n; float x[20],y[20],xx,yy; printf(Input n:); scanf(%d,n); if(n=20) {printf(Error! The value of n must in (0,20).); getch(); return 1;} if(n=0) {printf(Error! The value of n must in (0,20).);getch(); return 1;} for(i=0;i=n-1;i++) { printf(x[%d]:,i); ??? scanf(%f,x[i]); } ?? printf(\n); for(i=0;i=n-1;i++) { printf(y[%d]:,i);scanf(%f,y[i]);} printf(\n); difference(x,(float *)y,n); printf(Input xx:); scanf(%f,xx); yy=y[20]; for(i=n-1;i=0;i--) yy=yy*(xx-x[i])+y[i]; printf(NewtonInter(%f)=%f,xx,yy); getch(); } 3.高斯列主元消去法,求解其次线性方程组 C/C++ code #includestdio.h #include math.h #define N 20 int main() { int n,i,j,k; int mi,tmp,mx; float a[N][N],b[N],x[N]; printf(\nInput n:); scanf