103 SUBROUTINE slae2( A, B, C, RT1, RT2 )
111 REAL a,
b, c, rt1, rt2
118 parameter( one = 1.0e0 )
120 parameter( two = 2.0e0 )
122 parameter( zero = 0.0e0 )
124 parameter( half = 0.5e0 )
127 REAL ab, acmn, acmx, adf, df, rt, sm, tb
141 IF( abs( a ).GT.abs( c ) )
THEN
149 rt = adf*sqrt( one+( ab / adf )**2 )
150 ELSE IF( adf.LT.ab )
THEN
151 rt = ab*sqrt( one+( adf / ab )**2 )
158 IF( sm.LT.zero )
THEN
165 rt2 = ( acmx / rt1 )*acmn - (
b / rt1 )*
b
166 ELSE IF( sm.GT.zero )
THEN
173 rt2 = ( acmx / rt1 )*acmn - (
b / rt1 )*
b
subroutine slae2(A, B, C, RT1, RT2)
SLAE2 computes the eigenvalues of a 2-by-2 symmetric matrix.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i