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

#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，便于在下面输出时有

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