Esope source files

   1 : C RESYST    SOURCE    CHAT      05/01/13    02:55:04     5004           
   2 :       SUBROUTINE RESYST(A,B,M,X,LU)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 :       DIMENSION A(M,M),B(M),X(M)
   7 : C
   8 : C      PRINT*,A(LU,LU)
   9 :       X(LU)=B(LU)/A(LU,LU)
  10 : C
  11 :       DO 40 I=LU-1,1,-1
  12 : C
  13 :       TEMP=0.D0
  14 : C
  15 :       DO 50 J=I+1,LU
  16 :       TEMP=TEMP+A(I,J)*X(J)
  17 :    50 CONTINUE
  18 : C
  19 :       X(I)=(B(I)-TEMP)/A(I,I)
  20 : C   
  21 :    40 CONTINUE
  22 : C
  23 :       RETURN
  24 :       END
  25 : 
  26 :  
  27 :