埃特金迭代法求解方程的根.docVIP

埃特金迭代法求解方程的根 1．原理： 埃特金迭代法是利用两次迭代结果通过计算得到下一个x的值，这样加快了迭代式的收敛速度。 x’n+1=ψ(xn); x’’n+1=ψ(x’n+1); x n+1=x’’n+1-(x’’n+1-x’n+1)2/(x’’n+1-2x’n+1+xn); 利用这种方法求解方程的根会更快地得到结果。 2．源程序： #includemath.h #includeiostream.h int js,n; double eps,*x; double f(double x); int atkn(double* x,double eps,int js) { int flag,l; double u,v,x0; l=0; x0=*x; flag=0; while((flag==0)(l!=js)) { l=l+1; u=f(x0); v=f(u); if(fabs(u-v)eps) { x0=v; flag=1; } else x0=v-(v-u)*(v-u)/(v-2*u+x0); cout迭代第l次\txn=x0endl; } *x=x0; l=js-l; return(l); } double f(double x) { double s; s=6-x*x; return(s); }