求解a的N次方根

牛顿迭代法

$f(x)=f(x_0)+f^{{}^{\prime }}(x_0)\ast (x-x_0)$

假设$f(x)=x^N-a=0$，取N=2，得到
$f(x)=x^2-a{f}^{{}^{\prime }}(x)=2x$
$f(x)=x_0^2-a+2x_0(x-x_0)=0$
$x=\frac{1}{2}(x_0+\frac{a}{x_0})$
通过迭代可以逼近算出真实值。

对应任意的N，计算公式为
$x_{n+1}=\frac{1}{N}((N-1)x_n+\frac{a}{x_n^{N-1}})$

分析：
Newton迭代法求单根时，收敛速度很快(平方收敛)。
但如果方程根是重根，则收敛速度较慢，且重数越高速度越慢

代码：

#include <cmath>
#include <iostream>
#include <iomanip>
using namespace std;

double Newton(double x, int N, double delta) {
    if (x <= 0 || N <= 0) {
        return 0;
    }
    double x0 = 1;
    double x1 = ((N - 1) + x) / N;;

    while (abs(x1 - x0) > delta) {
        x0 = x1;
        x1 = ((N - 1)*x0 + x / pow(x0, N - 1)) / N;
    }
    return x1;
}

const double DELTA = 0.0000001;
int main() {

    cout << setprecision(20) << Newton(2, 2, DELTA) << endl;
    cout << setprecision(20) << Newton(8, 1, DELTA) << endl;
    cout << setprecision(20) << Newton(8, 2, DELTA) << endl;
    cout << setprecision(20) << Newton(8, 3, DELTA) << endl;
    cout << setprecision(20) << Newton(8, 4, DELTA) << endl;
    return 0;
}