PDF/spdf/aurpr.F

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