计算方法及程序实现

《计算方法及程序实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、一、对分法1、#include"math.h"main(){floata[3]={-3.0,1.0,0.0},b[3]={-2,0,l}，c[3I;floatf(floatx);inti;for(i=0;i<3;i++){do{c[i]=(a[i]+b[i])/2.0;if(f(c[i])==O)exit(O);elseif(f(c[i])*f(a[j])<0)b[i]=c[i];elsea[i]=c[";}while((b[i]-a[i])>le-5);}c[i]=a[i];printf("therootsare:"

2、)；for(i=0;i<3;i++)printf(n%f"，c[i]);printf("2.000000n);}floatf(floatx){floaty;y=x*x*x-2*x*x-4*x-7;return(y);}3、对分部分函数调用(题目要求如2)//include"math.h"floatfffloatx){floaty;y=x*x*x-2*x*x-4*x-7;return(y);}floatfl(floata'floatb){floatc;do{c=(a+b)/2;if(f(c)==O)exit(O);els

3、eif(f(a)*f(c)<0)b=c;elsea=c;}while((b-a)>le-5);return(a);}main(){floata=3.0,b=4.0,s;s=fl(a,b);printf("therootis%f",s);}2、用对分法求岀方程x3-2x2-4x-7=0在区间【3,4】内的根，精度要求为105。#include"math.h"main(){floata=3.0,b=4.0,c;floatf(floatx);do{c=(a+b)/2.0;if(f(c)==O){printf("theroot

4、is%f",c);exit(O);/*找到方程的根*/}elseif(f(c)*f(a)<0)b=c;elsea=c;}while((b-a)>le-5);printfC'therootis%f",a>;}floatf(floatx){floaty;y=x*x*x-2*x*x-4*x-7;return(y);}对分法的算法：扫描法的算法：4、对分法和扫描结合求方程x4-5x2+x+2=0的实根的上、下界，实现根的隔离，并用对分法求出所有的实根，精度要求为105。此方程的4个实根分别为：root=2.246979root

5、=0.554956root=0.801941root=2.000000include"math.h"floatf(floatx){floaty;y=x*x*x*x-5*x*x+x+2;returny;}main(){floata,b,c;inti=0;floatx,h=0.001,p[5】,q[5】,n;scanf("%f%f"/&a,&b)x=a;while(x

6、+h;n=i;}for(i=0;ile-5){c=(a+b)/2;if(f(c)==O){printf("%f",c);exit(0);}elseif(f(a)*f(c)<0)b=c;elsea=c;}printf("%f",c)}}二、秦九韶算法使用秦九韶算法计算多项式的值aOxn+alxln+alnx+an例如计算3x2+2x+l,当x=-l时值为2。main(){floataa[20],y,x;inti,n;printf("inputduox

7、iangshicishun:")scanf("%d",&n);printf("inputduoxiangshidexishu");for(i=0;i<=n;i++)scanf(H%f",&aa[i]);printf(ninputx");scanf("%fn,&xhy=aa[0];for(i=l;i<=n;i++)y=y*x+aa[i];printf(,,pn=°/of,,zy);}秦九韶算法三、牛顿法1、编写牛顿法求方程根的通用程序，精度要求为105。(x3-x2-2x-3=0初值x0=2)#inclu

8、de"math.h"main(){floatf(floatx);floatfd(floatx);floatx,x0,xl;xl=2.0;do{x0=xl;xl=x0-f(x0)/fd(x0);}while(fabs(xl-x0)>le-5);printf("therootis%f",xl);}floatf(floatx){re