牛顿法.txt

数值计算方法所用的上机实验C++源程序,个人所编.

#include <stdio.h>
#include <math.h>
#define eps 0.00001  /* 容许误差 */
float f(float x)  /* 定义函数f(x) */
{  return((-3*x+4)*x-5)*x+6;  }
float f1(float x)  /* 定义函数f(x)的导数 */
{  return (-9*x+8)*x-5;  }
main()
{  float  x0,x1=1.0;
   do
   {  x0=x1;   /* 准备下一次迭代的初值 */
      x1=x0-f(x0)/f1(x0);  /* 牛顿迭代 */
   }while(fabs(x1-x0)>eps); /*满足精度,输出近似根*/
   printf("x=%f\n",x1);
}
程序运行结果：
x=1.265328