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

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