C JACOD3    SOURCE    CHAT      05/01/13    00:48:32     5004           
      subroutine jacod3(a,idim,d)
C======================================================================         
C     OBJET                                                                     
C     -----                                                                     
C     DIAGONALISATION D UNE MATRICE 3*3 SYMETRIQUE                              
C     sans sortir les vecteurs propres
C                                                                               
C     ENTREES                                                                   
C     -------                                                                   
C     A(3,3) = MATRICE SYMETRIQUE                                               
C     IDIM   = 2 OU 3  SI 2 ON NE S OCCUPE QUE DE A(2,2)                        
C                      SI 3                    DE A(3,3)                        
C     SORTIES                                                                   
C     -------                                                                   
C     D(3)   = VALEURS PROPRES ORDONNEES D(1)>D(2)>D(3)
C
C===============================================================                
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      dimension a(3,3),d(3),x(3,3)
      if (idim.ne.3) then
        call jacob2(a,d,x)
        return
      endif
      c3=-1
      c2=a(1,1)+a(2,2)+a(3,3)
      c1=-(a(1,1)*a(2,2)+a(2,2)*a(3,3)+a(3,3)*a(1,1))+
     >  a(1,3)**2+a(1,2)**2+a(2,3)**2
      c0=2*a(1,2)*a(1,3)*a(2,3)-a(1,1)*a(2,3)**2-
     >  a(2,2)*a(1,3)**2-a(3,3)*a(1,2)**2+
     >  a(1,1)*a(2,2)*a(3,3)
      c0=c0/c3
      c1=c1/c3
      c2=c2/c3
      call degre3(c0,c1,c2,d1,XI1,d2,XI2,d3,XI3)
      d(1)=max(d1,d2,d3)
      d(3)=min(d1,d2,d3)
      d(2)=d1+d2+d3-d(1)-d(3)
      end
*