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

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:27:33  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:00  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
      SUBROUTINE AFPOLIN(XA,YA,N,X,Y,DY)
      PARAMETER (NMAX=10)
      DIMENSION XA(NMAX),YA(NMAX),C(NMAX),D(NMAX)
      Y=0.
      IF(N.GT.NMAX) RETURN
      NS=1
      DIF=ABS(X-XA(1))
      DO 11 I=1,N
        DIFT=ABS(X-XA(I))
        IF (DIFT.LT.DIF) THEN
          NS=I
          DIF=DIFT
        ENDIF
        C(I)=YA(I)
        D(I)=YA(I)
11    CONTINUE
      Y=YA(NS)
      NS=NS-1
      DO 13 M=1,N-1
        DO 12 I=1,N-M
          HO=XA(I)-X
          HP=XA(I+M)-X
          W=C(I+1)-D(I)
          DEN=HO-HP
C         IF(DEN.EQ.0.)PAUSE
          DEN=W/DEN
          D(I)=HP*DEN
          C(I)=HO*DEN
12      CONTINUE
        IF (2*NS.LT.N-M)THEN
          DY=C(NS+1)
        ELSE
          DY=D(NS)
          NS=NS-1
        ENDIF
        Y=Y+DY
13    CONTINUE
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:27:33 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:00 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	SUBROUTINE AFPOLIN(XA,YA,N,X,Y,DY)
	14	PARAMETER (NMAX=10)
	15	DIMENSION XA(NMAX),YA(NMAX),C(NMAX),D(NMAX)
	16	Y=0.
	17	IF(N.GT.NMAX) RETURN
	18	NS=1
	19	DIF=ABS(X-XA(1))
	20	DO 11 I=1,N
	21	DIFT=ABS(X-XA(I))
	22	IF (DIFT.LT.DIF) THEN
	23	NS=I
	24	DIF=DIFT
	25	ENDIF
	26	C(I)=YA(I)
	27	D(I)=YA(I)
	28	11 CONTINUE
	29	Y=YA(NS)
	30	NS=NS-1
	31	DO 13 M=1,N-1
	32	DO 12 I=1,N-M
	33	HO=XA(I)-X
	34	HP=XA(I+M)-X
	35	W=C(I+1)-D(I)
	36	DEN=HO-HP
	37	C IF(DEN.EQ.0.)PAUSE
	38	DEN=W/DEN
	39	D(I)=HP*DEN
	40	C(I)=HO*DEN
	41	12 CONTINUE
	42	IF (2*NS.LT.N-M)THEN
	43	DY=C(NS+1)
	44	ELSE
	45	DY=D(NS)
	46	NS=NS-1
	47	ENDIF
	48	Y=Y+DY
	49	13 CONTINUE
	50	RETURN
	51	END