牛顿迭代法求根

问题描述:
使用牛顿迭代法求f(x)=x6-x-1=0在x0=1.5附近的根,误差范围为10-8
问题分析:
牛顿迭代法是指由x0处的切线方程改写成xk+1=xk - f(xk)/f '(xk)且x0=初值带入迭代,k=0,1,2,…

由题目f(x)求导可得f ’(x)=6x5-1;带入迭代公式可得x - (x6-x-1)/(6x5-1) 即可求得结果

牛顿迭代法
#include<cstdio>
#include<cmath>
double fact(double x) {
 return x-(pow(x,6)-x-1)/(6*pow(x,5)-1); //迭代公式
}
int main() {
 double x0=1.5;//初始值
 double x1=fact(x0);
 while(fabs(x1-x0)>=10e-9) { //迭代开始,误差范围10-9
  printf("%.8lf\n",x1);
  x0=x1;
  x1=fact(x0);
 }
}

运行结果:
牛顿迭代法求根

上一篇:一本通例题埃及分数—题解&&深搜的剪枝技巧总结


下一篇:【记录】写题一时爽,一直写题一直爽