资源预览内容
第1页 / 共3页
第2页 / 共3页
第3页 / 共3页
亲,该文档总共3页全部预览完了,如果喜欢就下载吧!
资源描述
不动点#include#include#define N 1000double fun1(double x)double y;y=(1/2)*(x+3/x);return y;double fun2(double x)double y;y=(1/2)*(1-3/(x*x);return y;main()double x0=2,e=0.001,x1,x,a,l;int k=0;l=fun2(x0);if(fabs(l)=1)printf(error);elsedok+;x1=fun1(x0);x=fabs(x1-x0);x0=x1;while(xe&kN)printf(too slowly);elseprintf(%f,%d”,x1,k);2#include#include#define C 1double fun1(double x)double y;y=pow(x,3)+2*pow(x,2)+10*x-20;return y;double fun2(double x)double y;y=3*pow(x,2)+4*x+10;return y;main()double x0=1.5,f0,f1,e=pow(10,-8),x1,a;int k,N=1000;for(k=1;k=N;k+)f0=fun1(x0);f1=fun2(x0);if(f1=0)printf(error);break;elsex1=x0-f0/f1;if(fabs(x1)C)a=fabs(x1-x0);elsea=fabs(x1-x0)/fabs(x1);if(aN)printf(too slowly);3#include#includedouble fun(double x)double y;y=pow(x,3)+2*pow(x,2)+10*x-20; return y;main()double x0,a=1.0,b=1.5,c;double e=pow(10,-8);int k=0;dok+;x0=(a+b)/2;c=fun(x0)*fun(a);if(ce);printf(%f,%dn”,x0,k);
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号