SUBROUTINE FKMLT
*-- Author : S.Burke / J.V. Morris
      SUBROUTINE FKMLT(C1,C2,C3)
**********************************************************************                                        
*                                                                    *                                        
* Multiplies two 5*5 symmetric matrices                              *                                        
*                                                                    *                                        
* Output is C3 = C1.C2 (note that C3 is NOT symmetric!!!)            *                                        
*                                                                    *                                        
**********************************************************************                                        
                                                                        
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)                               
      DIMENSION C1(5,5),C2(5,5),C3(5,5)                                 
                                                                        
**********************************************************************                                        
                                                                        
      C3(1,1) = C1(1,1)*C2(1,1) + C1(2,1)*C2(2,1) + C1(3,1)*C2(3,1)     
     &        + C1(4,1)*C2(4,1) + C1(5,1)*C2(5,1)                       
      C3(2,1) = C1(2,1)*C2(1,1) + C1(2,2)*C2(2,1) + C1(3,2)*C2(3,1)     
     &        + C1(4,2)*C2(4,1) + C1(5,2)*C2(5,1)                       
      C3(3,1) = C1(3,1)*C2(1,1) + C1(3,2)*C2(2,1) + C1(3,3)*C2(3,1)     
     &        + C1(4,3)*C2(4,1) + C1(5,3)*C2(5,1)                       
      C3(4,1) = C1(4,1)*C2(1,1) + C1(4,2)*C2(2,1) + C1(4,3)*C2(3,1)     
     &        + C1(4,4)*C2(4,1) + C1(5,4)*C2(5,1)                       
      C3(5,1) = C1(5,1)*C2(1,1) + C1(5,2)*C2(2,1) + C1(5,3)*C2(3,1)     
     &        + C1(5,4)*C2(4,1) + C1(5,5)*C2(5,1)                       
                                                                        
      C3(1,2) = C1(1,1)*C2(2,1) + C1(2,1)*C2(2,2) + C1(3,1)*C2(3,2)     
     &        + C1(4,1)*C2(4,2) + C1(5,1)*C2(5,2)                       
      C3(2,2) = C1(2,1)*C2(2,1) + C1(2,2)*C2(2,2) + C1(3,2)*C2(3,2)     
     &        + C1(4,2)*C2(4,2) + C1(5,2)*C2(5,2)                       
      C3(3,2) = C1(3,1)*C2(2,1) + C1(3,2)*C2(2,2) + C1(3,3)*C2(3,2)     
     &        + C1(4,3)*C2(4,2) + C1(5,3)*C2(5,2)                       
      C3(4,2) = C1(4,1)*C2(2,1) + C1(4,2)*C2(2,2) + C1(4,3)*C2(3,2)     
     &        + C1(4,4)*C2(4,2) + C1(5,4)*C2(5,2)                       
      C3(5,2) = C1(5,1)*C2(2,1) + C1(5,2)*C2(2,2) + C1(5,3)*C2(3,2)     
     &        + C1(5,4)*C2(4,2) + C1(5,5)*C2(5,2)                       
                                                                        
      C3(1,3) = C1(1,1)*C2(3,1) + C1(2,1)*C2(3,2) + C1(3,1)*C2(3,3)     
     &        + C1(4,1)*C2(4,3) + C1(5,1)*C2(5,3)                       
      C3(2,3) = C1(2,1)*C2(3,1) + C1(2,2)*C2(3,2) + C1(3,2)*C2(3,3)     
     &        + C1(4,2)*C2(4,3) + C1(5,2)*C2(5,3)                       
      C3(3,3) = C1(3,1)*C2(3,1) + C1(3,2)*C2(3,2) + C1(3,3)*C2(3,3)     
     &        + C1(4,3)*C2(4,3) + C1(5,3)*C2(5,3)                       
      C3(4,3) = C1(4,1)*C2(3,1) + C1(4,2)*C2(3,2) + C1(4,3)*C2(3,3)     
     &        + C1(4,4)*C2(4,3) + C1(5,4)*C2(5,3)                       
      C3(5,3) = C1(5,1)*C2(3,1) + C1(5,2)*C2(3,2) + C1(5,3)*C2(3,3)     
     &        + C1(5,4)*C2(4,3) + C1(5,5)*C2(5,3)                       
                                                                        
      C3(1,4) = C1(1,1)*C2(4,1) + C1(2,1)*C2(4,2) + C1(3,1)*C2(4,3)     
     &        + C1(4,1)*C2(4,4) + C1(5,1)*C2(5,4)                       
      C3(2,4) = C1(2,1)*C2(4,1) + C1(2,2)*C2(4,2) + C1(3,2)*C2(4,3)     
     &        + C1(4,2)*C2(4,4) + C1(5,2)*C2(5,4)                       
      C3(3,4) = C1(3,1)*C2(4,1) + C1(3,2)*C2(4,2) + C1(3,3)*C2(4,3)     
     &        + C1(4,3)*C2(4,4) + C1(5,3)*C2(5,4)                       
      C3(4,4) = C1(4,1)*C2(4,1) + C1(4,2)*C2(4,2) + C1(4,3)*C2(4,3)     
     &        + C1(4,4)*C2(4,4) + C1(5,4)*C2(5,4)                       
      C3(5,4) = C1(5,1)*C2(4,1) + C1(5,2)*C2(4,2) + C1(5,3)*C2(4,3)     
     &        + C1(5,4)*C2(4,4) + C1(5,5)*C2(5,4)                       
                                                                        
      C3(1,5) = C1(1,1)*C2(5,1) + C1(2,1)*C2(5,2) + C1(3,1)*C2(5,3)     
     &        + C1(4,1)*C2(5,4) + C1(5,1)*C2(5,5)                       
      C3(2,5) = C1(2,1)*C2(5,1) + C1(2,2)*C2(5,2) + C1(3,2)*C2(5,3)     
     &        + C1(4,2)*C2(5,4) + C1(5,2)*C2(5,5)                       
      C3(3,5) = C1(3,1)*C2(5,1) + C1(3,2)*C2(5,2) + C1(3,3)*C2(5,3)     
     &        + C1(4,3)*C2(5,4) + C1(5,3)*C2(5,5)                       
      C3(4,5) = C1(4,1)*C2(5,1) + C1(4,2)*C2(5,2) + C1(4,3)*C2(5,3)     
     &        + C1(4,4)*C2(5,4) + C1(5,4)*C2(5,5)                       
      C3(5,5) = C1(5,1)*C2(5,1) + C1(5,2)*C2(5,2) + C1(5,3)*C2(5,3)     
     &        + C1(5,4)*C2(5,4) + C1(5,5)*C2(5,5)                       
                                                                        
      RETURN                                                            
      END                                                               
*