求平方根或立方根

1.求平方根：

命f(x) = x^2 - a，求解f(x) = x^2 - a = 0。

利用泰勒公式展开，即f(x)在xo处的函数值为：f(x) = f(xo) +f'(xo)(x-xo) = xo^2-a+2xo(x-x0) = 0，

解之得：x = (x+a/xo) / 2。

#include <math.h>
double sqr(double a)
{
    double x = 1.0;
    while(fabs(x * x - a) > 1e-9)
	{
        x = (x + a / x) / 2;
    }
    return x;
}

2.求立方根：

命f(x) = x^3 - a，求解f(x) = x^3 - a = 0。

利用泰勒公式展开，即f(x)在xo处的函数值为：f(x) = f(xo) +f'(xo)(x-xo) = xo^3-a+3xo^2(x-x0) = 0，

解之得：x = xo - (xo^3 - a) / (3xo^2)。

double fun(double a) {
    double x = 1.0;
    while(fabs(x*x*x - a) > 1e-9)
        x = x - ((x*x*x - a) / (3*x*x));
    return x;
}