牛顿法求平方根

对于给定的输入N,求平方根的公式是 x(k+1) =   (xk + N/xk)/2;

需要设置一个 N的平方根的初始估计值。取N的值就可以。

#include <stdio.h>

double newton(double d)
{
        if( d < 0)
        {
                return -1;
        }

        double ret = d;

        while( ret * ret -d > 0.00001 || ret * ret - d < -0.00001)
        {
                ret = (ret + d/ret ) / 2;
        }

        return ret;
}

int main(void)
{
        double d;
        scanf("%lf", &d);
        double ret = newton(d);
        printf("%f\n", ret);
        return 0;
}