1 SUBROUTINE dlae2( A, B, C, RT1, RT2 )
9 DOUBLE PRECISION A, B, C, RT1, RT2
57 parameter( one = 1.0d0 )
59 parameter( two = 2.0d0 )
61 parameter( zero = 0.0d0 )
63 parameter( half = 0.5d0 )
66 DOUBLE PRECISION AB, ACMN, ACMX, ADF, DF, RT, SM, TB
80 IF( abs( a ).GT.abs( c ) )
THEN
88 rt = adf*sqrt( one+( ab / adf )**2 )
89 ELSE IF( adf.LT.ab )
THEN
90 rt = ab*sqrt( one+( adf / ab )**2 )
104 rt2 = ( acmx / rt1 )*acmn - ( b / rt1 )*b
105 ELSE IF( sm.GT.zero )
THEN
112 rt2 = ( acmx / rt1 )*acmn - ( b / rt1 )*b
subroutine dlae2(A, B, C, RT1, RT2)