FUNCTION EPS(T,ANE,ALAM,ION,N)
C     ==============================
C
C     NLTE PARAMETER EPSILON (COLLISIONAL/SPONTANEOUS DEEXCITATION)
C     AFTER  KASTNER, 1981, J.Q.S.R.T. 26, 377
C
      INCLUDE 'PARAMS.FOR'
      DATA CK0,CK1 /7.75E-8, 2.58E-8/
      X=1.438E8/ALAM/T
      XKT=12390./ALAM
      TT=0.75*X
      T1=TT+1.
      A=4.36E7*XKT*XKT/(1.-EXP(-X))
      IF(ION.EQ.1) GO TO 10
      B=1.1+LOG(T1/TT)-0.4/T1/T1
      C=X*B*SQRT(T)/XKT/XKT*ANE
      IF(N.EQ.0) C=CK0*C
      IF(N.NE.0) C=CK1*C
      GO TO 20
   10 C=2.16/T/SQRT(T)/X**1.68*ANE
   20 EPS=C/(C+A)
      RETURN
      END