]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
7ef50f50 | 5 | * Revision 1.4 2000/09/18 10:02:35 hristov |
6 | * Makefile added to PDF8 | |
7 | * | |
21886bb6 | 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) | |
7ef50f50 | 18 | #if defined(CERNLIB_DOUBLE) |
19 | DOUBLE PRECISION | |
20 | #endif | |
21 | #if defined(CERNLIB_SINGLE) | |
22 | REAL | |
23 | #endif | |
21886bb6 | 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 |