SUBROUTINE EIGEN ( A, R, N, MV )
Argument Definitions (+ indicates altered content)
DOUBLE PRECISION  +A(*),      +R(*)
INTEGER            MV,         N
Description
a  = original matrix (symmetric), destroyed in computation. Resultant 
     eigenvalues are developed in diagonal of matrix (a) in descendin 
     order. 
r  = resultant matrix of eigenvectors (stored columnwise, in same 
     sequence as eigenvalues). 
n  = order of matrices (a) and (r). 
mv = 0  compute eigenvalues and eigenvectors. 
   = 1  compute eigenvalues only ((r) need not be dimensioned but mus 
        still appear in calling sequence).
Source file:eigen.f
Intrinsic Functions Called
DOUBLE PRECISION   DABS,       DBLE,       DSQRT
Local Variables (+ indicates altered content)
DOUBLE PRECISION  +ANORM,     +ANRMX,     +COSX,      +COSX2
DOUBLE PRECISION  +SINCS,     +SINX,      +SINX2,     +THR
DOUBLE PRECISION  +X,         +Y
INTEGER           +I,         +IA,        +IJ,        +IL
INTEGER           +ILQ,       +ILR,       +IM,        +IMQ
INTEGER           +IMR,       +IND,       +IQ,        +J
INTEGER           +JQ,        +K,         +L,         +LL
INTEGER           +LM,        +LQ,        +M,         +MM
INTEGER           +MQ