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