10个重要及算法C语言实现源代码.docVIP

窗体顶端 10个重要的算法C语言实现源代码：拉格朗日，牛顿插值，高斯，龙贝格，牛顿迭代，牛顿-科特斯，雅克比，秦九昭，幂法，高斯塞德尔 1.拉格朗日插值多项式?，用于离散数据的拟合  ?1#include?stdio.h?2?#include?conio.h?3?#include?alloc.h?4?float?lagrange(float?*x,float?*y,float?xx,int?n)?????/**//*拉格朗日插值算法*/?5?{?int?i,j;?6???float?*a,yy=0.0;????/**//*a作为临时变量，记录拉格朗日插值多项式*/?7???a=(float?*)malloc(n*sizeof(float));?8???for(i=0;i=n-1;i++)?9???{?a[i]=y[i];10?????for(j=0;j=n-1;j++)11?????if(j!=i)?a[i]*=(xx-x[j])/(x[i]-x[j]);12?????yy+=a[i];13???}14?free(a);15?return?yy;16}17main()18{?int?i,n;19?float?x[20],y[20],xx,yy;20?printf(Input?n:);21?scanf(%d,n);22?if(n=20)?{printf(Error!The?value?of?n?must?in?(0,20).);?getch();return?1;}23?if(n=0)?{printf(Error!?The?value?of?n?must?in?(0,20).);?getch();?return?1;}24?for(i=0;i=n-1;i++)25?{?printf(x[%d]:,i);26????scanf(%f,x[i]);27?}28?printf(\n);29?for(i=0;i=n-1;i++)30?{?printf(y[%d]:,i);scanf(%f,y[i]);}31?printf(\n);32?printf(Input?xx:);33?scanf(%f,xx);34?yy=lagrange(x,y,xx,n);35?printf(x=%f,y=%f\n,xx,yy);36?getch();37}38 2.牛顿插值多项式，用于离散数据的拟合 ?1#include?stdio.h?2#include?conio.h?3#include?alloc.h?4void?difference(float?*x,float?*y,int?n)?5{?float?*f;?6?int?k,i;?7?f=(float?*)malloc(n*sizeof(float));?8?for(k=1;k=n;k++)?9?{?f[0]=y[k];10????for(i=0;ik;i++)11??????f[i+1]=(f[i]-y[i])/(x[k]-x[i]);12????y[k]=f[k];13?}14?return;15}16main()17{?int?i,n;18?float?x[20],y[20],xx,yy;19?printf(Input?n:);20?scanf(%d,n);21?if(n=20)?{printf(Error!?The?value?of?n?must?in?(0,20).);?getch();?return?1;}22?if(n=0)?{printf(Error!?The?value?of?n?must?in?(0,20).);getch();?return?1;}23?for(i=0;i=n-1;i++)24?{?printf(x[%d]:,i);25????scanf(%f,x[i]);26?}27???printf(\n);28?for(i=0;i=n-1;i++)29?{?printf(y[%d]:,i);scanf(%f,y[i]);}30?printf(\n);31?difference(x,(float?*)y,n);32?printf(Input?xx:);33?scanf(%f,xx);34?yy=y[20];35?for(i=n-1;i=0;i--)?yy=yy*(xx-x[i])+y[i];36?printf(NewtonInter(%f)=%f,xx,yy);37?getch();38}39 3.高斯列主元消去法,求解其次线性方程组  ??1#includestdio.h??2#include?math.h?