Esope source files

   1 : C RESOLU    SOURCE    CHAT      05/01/13    02:53:07     5004           
   2 :       SUBROUTINE RESOLU (AAA,FF,XX,KERRE)
   3 : c  Resolution d'un systeme 2*2 
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 : C
   7 : C=================================================================
   8 : C
   9 : C ENTREES:
  10 : C --------
  11 : C
  12 : C  AAA    : Matrice 
  13 : C  FF     : 2nd membre 
  14 : C 
  15 : C SORTIES:
  16 : C --------
  17 : C
  18 : C XX      : Solution 
  19 : C
  20 : C================================================================
  21 :       DIMENSION AAA(2,2),XX(2),FF(2) 
  22 :       PARAMETER (EPSILO=1.D-10,UN=1.D0)
  23 :       KERRE = 0
  24 : C
  25 : C Determinant
  26 : C
  27 :       XDET = AAA(1,1)*AAA(2,2) - AAA(1,2)*AAA(2,1)
  28 : C
  29 :       IF (ABS(XDET).LT.EPSILO) THEN 
  30 :         KERRE = 2
  31 :         RETURN
  32 :       ELSE
  33 :         XDET = UN/XDET
  34 :       ENDIF
  35 : C
  36 :       XX(1) = XDET*(AAA(2,2)*FF(1)-AAA(2,1)*FF(2))
  37 :       XX(2) = XDET*(-AAA(1,2)*FF(1)+AAA(1,1)*FF(2))
  38 : C===========================================================
  39 : C                     Fin
  40 : C===========================================================
  41 :       RETURN
  42 :       END
  43 :  
  44 :