private static double Pythag( double a, double b )
{
double at = System.Math.Abs( a ), bt = System.Math.Abs( b ), ct, result;
if ( at > bt )
{
ct = bt / at;
result = at * System.Math.Sqrt( 1.0 + ct * ct );
}
else if ( bt > 0.0 )
{
ct = at / bt;
result = bt * System.Math.Sqrt( 1.0 + ct * ct );
}
else
{
result = 0.0;
}
return result;
}
}