PDF/spdf/aurpr.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:27:37  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:04  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE AURPR(I,NDRV,X,S,ANS)
  14 #include "pdf/expdp.inc"
  15      +       CALCO(8,20,25)
  16       COMMON/W5051I3/CALCO
  17       REAL   F1(25),F2(25)
  18       DATA DELTA/.10/
  19       ANS=0.
  20       IF(X.GT.0.9985) RETURN
  21       IF(I.EQ.3.AND.X.GT.0.95) RETURN
  22       IF(I.EQ.8.AND.X.GT.0.95) RETURN
  23       IS=S/DELTA+1
  24       IS1=IS+1
  25       DO 1 L=1,25
  26       KL=L+NDRV*25
  27       F1(L)=CALCO(I,IS,KL)
  28       F2(L)=CALCO(I,IS1,KL)
  29     1 CONTINUE
  30       A1=AUGETFV(X,F1)
  31       A2=AUGETFV(X,F2)
  32       S1=(IS-1)*DELTA
  33       S2=S1+DELTA
  34       ANS=A1*(S-S2)/(S1-S2)+A2*(S-S1)/(S2-S1)
  35 C      ANS=EXP(ANS)
  36       RETURN
  37       END