Esope source files

   1 : C LINEQ1    SOURCE    CHAT      05/01/13    01:18:44     5004           
   2 :       SUBROUTINE LINEQ1(B,D,X)
   3 : C===================================================
   4 : C                                                  =
   5 : C     RESOLUTION DES EQUATIONS LINEAIRES           =
   6 : C       B*D=X                                      =
   7 : C     OUTPUT: X                                    =
   8 : C                                                  =
   9 : C===================================================
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)
  12 :       DIMENSION B(2,2),D(2),X(2)
  13 :       DATA EPSI/1.D-9/
  14 :       DET=B(1,1)*B(2,2)-B(1,2)*B(2,1)
  15 :       X(1)= (D(1)*B(2,2)-D(2)*B(1,2))/DET
  16 :       X(2)=-(D(1)*B(2,1)-D(2)*B(1,1))/DET
  17 :       IF(ABS(DET).LT.EPSI) PRINT*,'LINEQ1: |DET|<1.D-9'
  18 :       END
  19 :  
  20 :