*-- Author : S.Burke / J.V. Morris
      SUBROUTINE FKSUB(C1,C2,C3)
**********************************************************************                                        
*                                                                    *                                        
* Subtracts two 5*5 symmetric matrices                               *                                        
*                                                                    *                                        
* Output is C3 = C1 - C2 (lower half only)                           *                                        
*                                                                    *                                        
**********************************************************************                                        
                                                                        
      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)                                       
      C3(2,1) = C1(2,1) - C2(2,1)                                       
      C3(3,1) = C1(3,1) - C2(3,1)                                       
      C3(4,1) = C1(4,1) - C2(4,1)                                       
      C3(5,1) = C1(5,1) - C2(5,1)                                       
                                                                        
      C3(2,2) = C1(2,2) - C2(2,2)                                       
      C3(3,2) = C1(3,2) - C2(3,2)                                       
      C3(4,2) = C1(4,2) - C2(4,2)                                       
      C3(5,2) = C1(5,2) - C2(5,2)                                       
                                                                        
      C3(3,3) = C1(3,3) - C2(3,3)                                       
      C3(4,3) = C1(4,3) - C2(4,3)                                       
      C3(5,3) = C1(5,3) - C2(5,3)                                       
                                                                        
      C3(4,4) = C1(4,4) - C2(4,4)                                       
      C3(5,4) = C1(5,4) - C2(5,4)                                       
                                                                        
      C3(5,5) = C1(5,5) - C2(5,5)                                       
                                                                        
      RETURN                                                            
      END                                                               
*