7.3 求解非线性方程一个实根的埃特金迭代法 atkn.c

来自「许士良常用算法程序集C语言,包括c++一些常用算法代码」· C语言 代码 · 共 24 行

#include "math.h"
int atkn(x,eps,js,f)
int js;
double eps,*x,(*f)();
{ 
	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.0*u+x0);
    }
    *x=x0; l=js-l;
    return(l);
}