牛頓法與擬牛頓法

牛頓法

牛頓法求算術平方根

float SqrtNewton(float x) {
    float x0 = x + 0.25, x1;
    for (;;)
    {
        x1 = (x0 * x0 + x) / (2 * x0);
        if (abs(x1 - x0) < FLT_EPSILON)
        {
            break;
        }
        x0 = x1;
    }
    return x1;
}

擬牛頓法