非线性方程的数值计算方法实验(17)

#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[

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库非线性方程的数值计算方法实验(17)在线全文阅读。