#include<float.h>//为了调用FLT_EPSILON，防止出现分母为零的情况。//
#define N 1000//保证有足够多的运算次数能达到比较精确的结果，且便于修改。//
#define pre 0.0000000001//使循环能够在达到某一误差精度后停下来。//
double main() {
double f(double t);//声明f函数。// double r(double t);//声明r函数。//
double t[N],err[N],relerr[N],dif[N];//使结果能够以数组形式输出，便于分析。//
char ch1[]="t",ch2[]="dif",ch3[]="err",ch4[]="relerr";//对输出进行标注。//
int k,i;//定义i，以便在达到所需精度之后，对所有结果以及对最终结果进行输出。//
t[0]=5.0,t[1]=15.0;//给出任意合法的初始值。// err[0]=fabs(t[1]-t[0]),relerr[0]=err[0]/fabs(t[1]+FLT_EPSILON),dif[0]=t[1]-t[0];//对各个结果的初始值进行运算。//
for(k=1;k<N;k++) {
t[k+1]=t[k]-(f(t[k])*(t[k]-t[k-1]))/(f(t[k])-f(t[