]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |