欢迎来到天天文库
浏览记录
ID:8959574
大小:15.50 KB
页数:2页
时间:2018-04-13
《c语言算法锦集一元非线性方程求根》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、一元非线性方程求根:/*牛顿迭代法求函数跟*/#include#includeintmain(void){doublex,x1,eps=1e-6,f,f1;/*误差为eps*/x=1.0;/*x=1.0是初值*/do{x1=x;f=6-x1*(5-x1*(4-3*x1));/*f为f(x)函数*/f1=-5+x1*(8-9*x1);/*f1为f(x)的导函数*/x=x1-f/f1;f=6-x*(5-x*(4-3*x));}while(fabs(f)>=eps&&fabs(x-x1)>=eps);prin
2、tf("x=%f",x);}/*二分法求函数跟*/#include#includedoublef(doublex){return6-x*(5-x*(4-3*x));/*f(x)函数*/}intmain(void){doublea,b,c,x,eps=1e-6;do{scanf("%lf%lf",&a,&b);}while(f(a)*f(b)>0);if(fabs(f(a))<1e-6)x=a;elseif(fabs(f(b))<1e-6)x=b;else{c=(b+a)/2;while(fabs(f(c
3、))>eps&&fabs(b-a)>eps){if(f(a)*f(c)<0)b=c;elsea=c;c=(b+a)/2;}x=c;}printf("x=%f",x);}/*弦截法求函数跟*/c=(a*f(b)-b*f(a))/(f(b)-f(a));while(fabs(f(c))>eps){if(f(a)*f(c)<0)b=c;elsea=c;c=(a*f(b)-b*f(a))/(f(b)-f(a));}#includevoidf();intmain(void){intx,loop=0;do{for(x=1;x<5;
4、x++){intx=2;printf("%d",x);}printf("%d",x);f();loop++;}while(loop<1);getch();}voidf(){printf("%d",x++);}
此文档下载收益归作者所有