《数值分析不动点》由会员分享,可在线阅读,更多相关《数值分析不动点(3页珍藏版)》请在金锄头文库上搜索。
1、不动点#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);e
2、lseprintf(%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);bre
3、ak;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);