···65656666 /* Parts of the quadratic formula solution */
6767 const double b = u + 1.0;
6868- const double root = sqrt(b + 4 * u * v * v);
6868+ const double root = sqrt(b * b + 4 * u * v * v);
6969 const double two_a = 2 * v;
70707171 /* The two possible solutions. */