64 DOUBLE PRECISION FUNCTION dlapy2( X, Y )
79 parameter( zero = 0.0d0 )
81 parameter( one = 1.0d0 )
84 DOUBLE PRECISION w, xabs, yabs, z
87 INTRINSIC abs, max, min, sqrt
98 dlapy2 = w*sqrt( one+( z / w )**2 )
double precision function dlapy2(X, Y)
DLAPY2 returns sqrt(x2+y2).