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

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:27:35  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:02  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
      SUBROUTINE AURGAM(I,NDRV,X,S,ANS)
      DIMENSION F1(32),F2(32),F3(32)
      DIMENSION AF(3),AS(3)
      DIMENSION CALCO(8,20,32)
      COMMON/W5051I7/CALCO
      DATA DELTA/0.8000E-01/
      ANS=0.
      IF(X.GT.0.9985) RETURN
      N=3
      IS=S/DELTA+1
      IF(IS.GE.17) IS=17
      IS1=IS+1
      IS2=IS1+1
      DO 1 L=1,32
      KL=L+32*NDRV
      F1(L)=CALCO(I,IS,KL)
      F2(L)=CALCO(I,IS1,KL)
      F3(L)=CALCO(I,IS2,KL)
    1 CONTINUE
      AF(1)=AFGETFV(X,F1)
      AF(2)=AFGETFV(X,F2)
      AF(3)=AFGETFV(X,F3)
      AS(1)=(IS-1)*DELTA
      AS(2)=AS(1)+DELTA
      AS(3)=AS(2)+DELTA
      CALL AFPOLIN(AS,AF,N,S,AANS,DY)
      ANS=AANS
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:27:35 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:02 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	SUBROUTINE AURGAM(I,NDRV,X,S,ANS)
	14	DIMENSION F1(32),F2(32),F3(32)
	15	DIMENSION AF(3),AS(3)
	16	DIMENSION CALCO(8,20,32)
	17	COMMON/W5051I7/CALCO
	18	DATA DELTA/0.8000E-01/
	19	ANS=0.
	20	IF(X.GT.0.9985) RETURN
	21	N=3
	22	IS=S/DELTA+1
	23	IF(IS.GE.17) IS=17
	24	IS1=IS+1
	25	IS2=IS1+1
	26	DO 1 L=1,32
	27	KL=L+32*NDRV
	28	F1(L)=CALCO(I,IS,KL)
	29	F2(L)=CALCO(I,IS1,KL)
	30	F3(L)=CALCO(I,IS2,KL)
	31	1 CONTINUE
	32	AF(1)=AFGETFV(X,F1)
	33	AF(2)=AFGETFV(X,F2)
	34	AF(3)=AFGETFV(X,F3)
	35	AS(1)=(IS-1)*DELTA
	36	AS(2)=AS(1)+DELTA
	37	AS(3)=AS(2)+DELTA
	38	CALL AFPOLIN(AS,AF,N,S,AANS,DY)
	39	ANS=AANS
	40	RETURN
	41	END