rtnewt.c

来自「一些常用数值算法的 C程序源代码 对初级学数值计算很有帮助」· C语言 代码 · 共 24 行

#include <math.h>
#define JMAX 20

float rtnewt(void (*funcd)(float, float *, float *), float x1, float x2,
	float xacc)
{
	void jswlerror(char error_text[]);
	int j;
	float df,dx,f,rtn;

	rtn=0.5*(x1+x2);
	for (j=1;j<=JMAX;j++) {
		(*funcd)(rtn,&f,&df);
		dx=f/df;
		rtn -= dx;
		if ((x1-rtn)*(rtn-x2) < 0.0)
			jswlerror("Jumped out of brackets in rtnewt");
		if (fabs(dx) < xacc) return rtn;
	}
	jswlerror("Maximum number of iterations exceeded in rtnewt");
	return 0.0;
}
#undef JMAX