MINICERN/mathlib/gen/c/qfunct.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:00  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE QFUNCT(XK,XKPRIM,Q)
  11       XKPRIM=SQRT(1.-XK**2)
  12       A=MAX(ABS(XK),XKPRIM)
  13       EPSI=(1.-A)/(2.*A+4.*SQRT(A)+2.)
  14       E=EPSI**4
  15       Q=EPSI+EPSI*E*(2.+E*(15.+E*150.))
  16       IF(XKPRIM.LT.A) Q=EXP(9.869604401089359/LOG(Q))
  17       RETURN
  18       END