PDF/spdf/aurpi3.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:03  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16 C
  17       SUBROUTINE AURPI3(I,NDRV,X,S,ANS)
  18 #if defined(CERNLIB_DOUBLE)
  19       DOUBLE PRECISION
  20 #endif
  21 #if defined(CERNLIB_SINGLE)
  22       REAL
  23 #endif
  24      +       CALCPIO(8,20,25)
  25       COMMON/W5051I6/CALCPIO
  26       REAL   F1(25),F2(25)
  27       DATA DELTA/.10/
  28       ANS=0.
  29       IF(X.GT.0.9985) RETURN
  30       IF(I.EQ.3.AND.X.GT.0.95) RETURN
  31       IF(I.EQ.8.AND.X.GT.0.95) RETURN
  32       IS=S/DELTA+1
  33       IS1=IS+1
  34       DO 1 L=1,25
  35       KL=L+NDRV*25
  36       F1(L)=CALCPIO(I,IS,KL)
  37       F2(L)=CALCPIO(I,IS1,KL)
  38     1 CONTINUE
  39       A1=AUGETFV(X,F1)
  40       A2=AUGETFV(X,F2)
  41       S1=(IS-1)*DELTA
  42       S2=S1+DELTA
  43       ANS=A1*(S-S2)/(S1-S2)+A2*(S-S1)/(S2-S1)
  44       RETURN
  45       END