牛顿下山法

牛顿下山法计算高次方程真是利器啊,举个栗子

#include <limits>

#include <iostream>

#include <math.h>

double f(double x)

{

return pow(x, 3) + pow(x, 2) + 2;

}

double f1(double x)

{

return 3 * pow(x , 2) + 2 * x;

}

int main()

{

double a0 = 5.0;

std::cin >> a0;

double da = 0.0;

do

{

double a = a0 - f(a0) / f1(a0);

da = a - a0;

a0 = a;

} while (abs(da) > 0.001);

std::cout << a0 + da << std::endl;

system("pause");

return -1;

}