#include<math.h>
#include<float.h>//为了调用FLT_EPSILON，防止出现分母为零的情况。//
#define pre 0.00000000001//使循环能够在达到某一误差精度后停下来。//
#define N 1000//迭代次数，以保证能得到最精确的值。// double main() {
double fa(double x);//声明fa函数。// double fb(double x);//声明fb函数。// double fc(double x);//声明fc函数。//
double a[N],b[N],c[N],err[N];//定义结果数组，与误差数组。//
int k,i,j,m;
a[0]=1.0,b[0]=1.0,c[0]=2.0;//赋初始值，由书上给出。// for(k=0;k<N;k++) {
a[k+1]=a[k]-5.0*fa(a[k]);//加速牛顿-拉夫森公式，见书上P64。下同。//
err[k]=fabs(2.0-a[k+2]);//对绝对误差的运算。k加2，以保证下面输出时其有足够多的值输出。//
i=k;//将循环停止时的k值赋给i，便于在下面输出时有