]>
Commit | Line | Data |
---|---|---|
21886bb6 | 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 |