把f(x)在x0附近展开成泰勒级数f(x) = f(x0)+f'(x0)(x-x0)+f''(x0)/2!*(x-x0)^2+...然后取其线性部分,作为非线性方程f(x) = 0的
近似方程,即泰勒展开的前两项,则有
f(x) = f'(x0)x - x0*f'(x0) + f(x0) = 0
f'(x0)x = x0*f'(x0) - f(x0)
x = x0 - f(x0)/f'(x0)
得到牛顿的一个迭代序列:
->x(n+1) = x(n)-f(x(n))/f'(x(n))
近似方程,即泰勒展开的前两项,则有
f(x) = f'(x0)x - x0*f'(x0) + f(x0) = 0
f'(x0)x = x0*f'(x0) - f(x0)
x = x0 - f(x0)/f'(x0)
得到牛顿的一个迭代序列:
->x(n+1) = x(n)-f(x(n))/f'(x(n))
例:求方程f(x) = 2*x^3-4*x^2+3*x-6在x = 1.5附近的根,求出总的迭代次数和每次的近似根。
#include <iostream> #include <cmath> using namespace std; double fun1(double x) { return 2*x*x*