[u/mrichter/AliRoot.git] / PDF / spdf / aurpi3.F

*
* $Id$
*
* $Log$
* 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:03  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
C
      SUBROUTINE AURPI3(I,NDRV,X,S,ANS)
#include "pdf/expdp.inc"
     +       CALCPIO(8,20,25)
      COMMON/W5051I6/CALCPIO
      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)=CALCPIO(I,IS,KL)
      F2(L)=CALCPIO(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)
      RETURN
      END
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:03 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	C
	14	SUBROUTINE AURPI3(I,NDRV,X,S,ANS)
	15	#include "pdf/expdp.inc"
	16	+ CALCPIO(8,20,25)
	17	COMMON/W5051I6/CALCPIO
	18	REAL F1(25),F2(25)
	19	DATA DELTA/.10/
	20	ANS=0.
	21	IF(X.GT.0.9985) RETURN
	22	IF(I.EQ.3.AND.X.GT.0.95) RETURN
	23	IF(I.EQ.8.AND.X.GT.0.95) RETURN
	24	IS=S/DELTA+1
	25	IS1=IS+1
	26	DO 1 L=1,25
	27	KL=L+NDRV*25
	28	F1(L)=CALCPIO(I,IS,KL)
	29	F2(L)=CALCPIO(I,IS1,KL)
	30	1 CONTINUE
	31	A1=AUGETFV(X,F1)
	32	A2=AUGETFV(X,F2)
	33	S1=(IS-1)*DELTA
	34	S2=S1+DELTA
	35	ANS=A1(S-S2)/(S1-S2)+A2(S-S1)/(S2-S1)
	36	RETURN
	37	END