Sources Esope | Cast3M

xpuiss
C XPUISS    SOURCE    CHAT      05/01/13    04:14:52     5004
       FUNCTION XPUISS(X,A,AMAX)
C
      IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)
       IF (X.LE.0.0) THEN
       XPUISS = 0.0
       IF (A.LT.0.0) XPUISS = AMAX
       IF (A.EQ.0.0) XPUISS = 1.0
       ELSE
       XX = LOG(X)*A
       XM = LOG(AMAX)
       IF (XX.GT.XM) THEN
       XPUISS = AMAX
       ELSEIF (XX.LE.-XM) THEN
       XPUISS = 0.0
       ELSE
       XPUISS = EXP(XX)
       ENDIF
       ENDIF
       RETURN
       END