This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / g / vavfcn.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/04/01 15:02:47  mclareni
6 * Mathlib gen
7 *
8 *
9 #include "gen/pilot.h"
10       FUNCTION VAVFCN(X,RKA,BE2)
11 C
12       COMMON /VAVILI/ T0,T1,T,OMEGA
13       COMMON /FORFCN/ SS,LFCN
14       DATA E,PI,RG /5E-4, 3.1415926535898, 0.5772156649015/
15 C
16       GO TO (1,2,3), LFCN
17 C
18 C     FOR LFCN=1 FCN IS USED TO DETERMINE THE NUMBER N OF FOURIER
19 C     COEFFICIENTS IN DISVAV(X,I)
20 C
21     1 RN=5.
22       IF(RKA .GE. 0.07) RN=8.
23       VAVFCN=RKA*(2.+BE2*RG)-1.596312592
24      1   +BE2*RKA*LOG(OMEGA*X/RKA)-0.5*OMEGA*PI*X+2.302585093*RN
25       RETURN
26 C
27 C     FOR LFCN=2 FCN IS USED TO DETERMINE THE SS-QUANTILE OF THE
28 C     VAVILOV DISTRIBUTION
29 C
30     2 VAVFCN=DISVAV(X,1)-SS
31       RETURN
32 C
33 C     FOR LFCN=3 FCN IS USED TO DETERMINE THE RIGHT-HAND ENDPOINT T+
34 C     OF THE SUPPORT OF DISVAV(X,0)
35 C
36     3 VAVFCN=1.-BE2*(1.-RG)-LOG(E)/RKA-X+BE2*(LOG(ABS(X))+EXPINT(X))
37      1    -(1.-BE2)*EXP(-X)
38       RETURN
39 C
40       END