MINICERN/mathlib/gen/g/disvav.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 DISVAV(X,I)
  11 C
  12 C     DISVAV(X,I) COMPUTES FOR I=0 THE VALUE OF THE DENSITY, FOR I=1
  13 C     THE VALUE OF THE CUMULATIVE DISTRIBUTION FUNCTION OF THE
  14 C     VAVILOV DISTRIBUTION AT THE POINT X
  15 C
  16       COMMON /VAVILI/ T0,T1,T,OMEGA
  17       COMMON /VAVILO/ A(155),B(155),N
  18 C
  19       IF(X .LT. T0) GO TO 3
  20       IF(X .GT. T1) GO TO 4
  21 C
  22       Y=X-T0
  23 C
  24       TRR=0.
  25       IF(I .EQ. 1) TRR=Y/T
  26 C
  27       Z=OMEGA*Y-3.1415926535898
  28       COF=2.0*COS(Z)
  29       A1=0.
  30       A0=A(1)
  31       NN=N-1
  32       DO 1 K = 2,N
  33       A2=A1
  34       A1=A0
  35     1 A0=A(K)+COF*A1-A2
  36 C
  37       TRR=TRR+0.5*(A0-A2)
  38       A1=0.
  39       A0=B(1)
  40       DO 2 K = 2,NN
  41       A2=A1
  42       A1=A0
  43     2 A0=B(K)+COF*A1-A2
  44 C
  45       DISVAV=TRR+A0*SIN(Z)
  46       RETURN
  47 C
  48     3 DISVAV=0.0
  49       RETURN
  50 C
  51     4 DISVAV=I
  52       RETURN
  53 C
  54       END