! ! $Header$ ! SUBROUTINE eigen( e,d) IMPLICIT NONE INCLUDE "dimensions.h" real :: e( iim,iim ), d( iim ) real :: asm( iim ) integer :: im,i,j im=iim c DO i = 1,im asm( i ) = d( im-i+1 ) END DO DO i = 1,iim d( i ) = asm( i ) END DO c c PRINT 70,d 70 FORMAT(5x,'Valeurs propres',/,8(1x,8f10.4,/),/) print * c DO i = 1,im DO j = 1,im asm( j ) = e( i , im-j+1 ) END DO DO j = 1,im e( i,j ) = asm( j ) END DO END DO RETURN END