resolu
C RESOLU SOURCE CHAT 05/01/13 02:53:07 5004 SUBROUTINE RESOLU (AAA,FF,XX,KERRE) c Resolution d'un systeme 2*2 IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) C C================================================================= C C ENTREES: C -------- C C AAA : Matrice C FF : 2nd membre C C SORTIES: C -------- C C XX : Solution C C================================================================ DIMENSION AAA(2,2),XX(2),FF(2) PARAMETER (EPSILO=1.D-10,UN=1.D0) KERRE = 0 C C Determinant C XDET = AAA(1,1)*AAA(2,2) - AAA(1,2)*AAA(2,1) C IF (ABS(XDET).LT.EPSILO) THEN KERRE = 2 RETURN ELSE XDET = UN/XDET ENDIF C XX(1) = XDET*(AAA(2,2)*FF(1)-AAA(2,1)*FF(2)) XX(2) = XDET*(-AAA(1,2)*FF(1)+AAA(1,1)*FF(2)) C=========================================================== C Fin C=========================================================== RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales