C SIMP3     SOURCE    CHAT      05/01/13    03:18:45     5004           
      SUBROUTINE SIMP3(A,MP,NP,I1,K1,IP,KP)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC        
C                                                                     C
C     METHODE DU SIMPLEX : ESCLAVE DE SIMPLX                          C
C                                                                     C
C     P.PEGON 31/8/92                                                 C
C                                                                     C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION A(MP,NP)
C
      PIV=1.D0/A(IP+1,KP+1)
      IF(I1.GE.0.D0)THEN
        DO 12 II=1,I1+1
          IF(II-1.NE.IP)THEN
            A(II,KP+1)=A(II,KP+1)*PIV
            DO 11 KK=1,K1+1
              IF(KK-1.NE.KP)THEN
                A(II,KK)=A(II,KK)-A(IP+1,KK)*A(II,KP+1)
              ENDIF
 11         CONTINUE
          ENDIF
 12     CONTINUE
      ENDIF
C
      DO 13 KK=1,K1+1
        IF(KK-1.NE.KP)A(IP+1,KK)=-A(IP+1,KK)*PIV
 13   CONTINUE
      A(IP+1,KP+1)=PIV
C
      RETURN
      END