* * $Id$ * * $Log$ * Revision 1.4 2000/09/18 10:02:35 hristov * Makefile added to PDF8 * * Revision 1.1.1.2 1996/10/30 08:27:37 cernlib * Version 7.04 * * Revision 1.1.1.1 1996/04/12 15:29:04 plothow * Version 7.01 * * #include "pdf/pilot.h" SUBROUTINE AURPR(I,NDRV,X,S,ANS) #if defined(CERNLIB_DOUBLE) DOUBLE PRECISION #endif #if defined(CERNLIB_SINGLE) REAL #endif + CALCO(8,20,25) COMMON/W5051I3/CALCO REAL F1(25),F2(25) DATA DELTA/.10/ ANS=0. IF(X.GT.0.9985) RETURN IF(I.EQ.3.AND.X.GT.0.95) RETURN IF(I.EQ.8.AND.X.GT.0.95) RETURN IS=S/DELTA+1 IS1=IS+1 DO 1 L=1,25 KL=L+NDRV*25 F1(L)=CALCO(I,IS,KL) F2(L)=CALCO(I,IS1,KL) 1 CONTINUE A1=AUGETFV(X,F1) A2=AUGETFV(X,F2) S1=(IS-1)*DELTA S2=S1+DELTA ANS=A1*(S-S2)/(S1-S2)+A2*(S-S1)/(S2-S1) C ANS=EXP(ANS) RETURN END