数值计算C语言常用小程序.doc

《数值计算C语言常用小程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1、秦九韶算法2、二分法3、拉格朗日插值4、埃特金算法5、复化梯形法6、复化辛甫生算法7、二阶龙格库塔方法8、四阶龙格库塔方法　9、改进的欧拉方法　10、迭代法　11、埃特金加速方法：12、牛顿迭代法　13、追赶法　14、雅克比迭代　15、蛋白质设计：17高斯消去法: 1、秦九韶算法P11.3利用秦九韶算法求多项式，在x=3时的值。程序：#include#includevoidmain(){floata[100],v,x;intn,i,k;scanf("%d%f",&n,&x);for(i=0;i<=n;i++)scanf("%f",&a[

2、i]);v=a[n];k=1;do{v=x*v+a[n-k];k=k+1;}while(k<=n);printf("v=%f",v);} 运行结果： 2、二分法P11.1用二分法求方程法x*x*x-x-1=0在[1,2]内的近似根，要求误差不超过#include#includefloatfun(float);voidmain(){floata,b,c,x,y,y1;scanf("%f%f%f",&a,&b,&c);y1=fun(a);do{x=(a+b)/2;y=fun(x);{if(y*y1>0)a=x;elseb=x;}}while((

3、b-a)>=c);printf("%f,%f",x,y);}floatfun(floatm){floatn;n=m*m*m-m-1;return(n);} 运行结果：  3、拉格朗日插值程序：#includemain(){floata,b,t,x[100],y[100];intn,i,j,k;scanf("%f%d",&a,&n);for(i=0;i<=n;i++)scanf("%f%f",&x[i],&y[i]);k=0;b=0;for(k=0;k<=n;k++){t=1;for(j=0;j<=n;j++){if(j!=k)t=t*(a-x[j])/

4、(x[k]-x[j]);}b=b+t*y[k];}printf("%f",b);}  4、埃特金算法程序：#include#includemain(){floata,b,c,x[100],y[100];inti,j,n,k;scanf("%d%f",&n,&a);for(i=0;i<=n;i++)scanf("%f%f",&x[i],&y[i]);for(k=1;k<=n;k++){for(i=k;i<=n;i++)y[i]=y[k-1]+(y[i]-y[k-1])*(a-x[k-1])/(x[i]-x[k-1]);}printf("%

5、f",y[n]);}5、复化梯形法P95.9设,用复化梯形法求积分的近似值程序：#include#includedoublefun(double);voidmain(){doublea,b,h,s,x,y;intn,k;scanf("%lf%lf%d",&a,&b,&n);h=(b-a)/n;s=0;x=a;y=fun(x);for(k=1;k<=n;k++){s=s+fun(x);x=x+h;s=s+fun(x);}s=(h/2)*s;printf("s=%lf",s);}doublefun(doublem){doublen;n=

6、exp(-m)*sin(4*m)+1;return(n);}运行结果：  6、复化辛甫生算法P95.9设,用复化辛甫生法求积分的近似值程序：#include#includedoublefun(double);voidmain(){doublea,b,h,s,x,y;intn,k;scanf("%lf%lf%d",&a,&b,&n);h=(b-a)/n;s=0;x=a;y=fun(x);for(k=1;k<=n;k++){s=s+fun(x);x=x+h/2;s=s+4*fun(x);x=x+h/2;s=s+fun(x);}s=(h/6)*s

7、;printf("s=%lf",s);}doublefun(doublem){doublen;n=exp(-m)*sin(4*m)+1;return(n);}运行结果：   7、二阶龙格库塔方法求解初值问题： 取h=0.2程序：#include#includefloatfun(float,float);voidmain(){floath,x0,y0,x1,y1,k1,k2;intn,N;scanf("%f%f%f%d",&x0,&y0,&h,&N);n=1;for(n=1;n<