Esope source files

   1 : C PROSY     SOURCE    CHAT      05/01/13    02:35:04     5004           
   2 :       SUBROUTINE PROSY(I,J,A,B,C)                                               
   3 : C  S.P.EFFECTUANT A=T(C)*B*C ,AVEC A(I,J),B(I,J),C(J,I)---(1/2 A SUP)--         
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   6 :       DIMENSION A(I,1),B(J,1),C(J,1)                                            
   7 :       DO 1 L=1,I                                                                
   8 :       DO 2 K=1,L                                                                
   9 :       A(K,L)=0.D0                                                               
  10 :    2  CONTINUE                                                                  
  11 :       DO 1 M=1,J                                                                
  12 :       S=0.D0                                                                    
  13 :       DO 3 N=1,J                                                                
  14 :       S=S+B(M,N)*C(N,L)                                                         
  15 :    3  CONTINUE                                                                  
  16 :       DO 1 K=1,L                                                                
  17 :       A(K,L)=A(K,L)+C(M,K)*S                                                    
  18 :    1  CONTINUE                                                                  
  19 :       RETURN                                                                    
  20 :       END                                                                       
  21 :