Esope source files

   1 : C RIMT6     SOURCE    CHAT      05/01/13    03:00:43     5004           
   2 :       SUBROUTINE RIMT6(RE,E,RAI,X13,X23,Y13,HM,A)                               
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION RE(6,1),E(3,1),D(3,6)                                           
   6 :       AD=0.25*RAI*HM/A                                                          
   7 :       DO1I=1,5,2                                                                
   8 :       D(2,I)=0.                                                                 
   9 :     1 D(1,I+1)=0.                                                               
  10 :       D(1,1)=Y13                                                                
  11 :       D(1,3)=-D(1,1)                                                            
  12 :       D(1,5)=0.                                                                 
  13 :       D(2,2)=-X23                                                               
  14 :       D(2,4)=X13                                                                
  15 :       D(2,6)=-X13+X23                                                           
  16 :       D(3,1)=D(2,2)                                                             
  17 :       D(3,2)=D(1,1)                                                             
  18 :       D(3,3)=D(2,4)                                                             
  19 :       D(3,4)=D(1,3)                                                             
  20 :       D(3,5)=D(2,6)                                                             
  21 :       D(3,6)=0.                                                                 
  21 :       CALL PROSY(6,3,RE,E,D)                                                    
  22 :       DO2I=1,6                                                                  
  23 :       DO2J=I,6                                                                  
  24 :     2 RE(I,J)=RE(I,J)*AD                                                        
  25 :       RETURN                                                                    
  26 :       END                                                                       
  27 :