实验4 牛顿迭代法.doc

实验4 牛顿迭代法求方程的一个实根 一、算法描述 参数说明： x 双精度实型变量指针。在该指针指向的变量中存放迭代初值；返回时在该指针指向的变量中存放终值。 eps 双精度实型变量。控制精度要求。 js 整型变量。最大迭代次数 f 无类型函数指针变量。指向计算方程左端函数f（x）值极其导数值f’（x）的函数。该函数由用户自编，其形式为 void f(x, y) double x, y[2]; { y[0]=f(x)的表达示； y[1]=f’(x)的表达示； return; } 本函数返回一个整型标志值。若返回的标志值小于0，则表示在出现了f’(x)=0的情况，函数中将输出信息“err”；若返回的标志值等于最大迭代次数，则表示还未满足精度要求，此时返回的实根终值只作为参考；若返回的标志值大于或等于0且小于最大迭代次数，则表示正常返回。 #include “stdio.h” #include “math.h” int newt(x, eps, js, f) int js; double *x, eps; void (*f)(); { int k;l; double y[2], d, p, x0, x1; l=js; k=1; x0=*x; (*f)(x0, y); /*计算f(x0)与f’(x0)*/ d=eps+1.0; while ((d=eps)(l!=0)) { if (fabs(y[1])+1.0==1.0) { printf(“err\n”); return(-1); } x1=x0-y[0]/y[1]; /*迭代计算x1=x0-f(x0)/f’(x0)*/ (*f)(x1, y); /*计算f(x1)与f’(x1)*/ d=fabs(x1-x0); p=fabs(y[0]); if (pd) d=p; x0=x1;l=l-1; } *x=x1; k=js-l; return(k); } 二、用牛顿迭代法求方程 f(x)=x3-x2-1=0 在x0=1.5附近的一个实根。取ε=0.000001，最大迭代次数为60。 其中f’(x)=3x2-2x。 主函数程序以及计算f(x)、f’(x)的函数程序如下： #include “math.h” #include “stdio.h” #include “newt.c” main() { int js, k; double x, eps; void newtf(); eps=0.000001; js=60; x=1.5; k=newt(x, eps, js, newtf); if (k=0) printf(“k=%d x=%13.7e\n”, k, x); printf(“\n”); } void newtf(x, y) double x, y[2]; { y[0]=x*x*(x-1.0)-1.0; y[1]=3.0*x*x-2.0*x; return; } #include “stdio.h” #include “math.h” int newt(double *x, double eps, int js,void (*f)()) { int k;l; double y[2], d, p, x0, x1; l=js; k=1; x0=*x; (*f)(x0, y); /*计算f(x0)与f’(x0)*/ d=eps+1.0; while ((d=eps)(l!=0)) { if (fabs(y[1])+1.0==1.0) { printf(“err\n”); return(-1); } x1=x0-y[0]/y[1]; /*迭代计算x1=x0-f(x0)/f’(x0)*/ (*f)(x1, y); /*计算f(x1)与f’(x1)*/ d=fabs(x1-x0); p=fabs(y[0]); if (pd) d=p; x0=x1;l=l-1; } *x=x1; k=js-l; return(k); } main() { int js, k; double x, eps; void newtf(); eps=0.000001; js=60; x=1.5; k=newt(x, eps, js, newtf); if (k=0) printf(“k=%d x=%13.7e\n”, k, x); printf(“\n”); } void newtf(x, y)