[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / rkstpcod.inc

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1996/04/01 15:02:16  mclareni
* Mathlib gen
*
*
*
* rkstpcod.inc
*
      DIMENSION Y(*),W(N,*)

      PARAMETER (Z1 = 1, HF = Z1/2, HS = Z1/6)

      IF(N .LE. 0) RETURN
      H2=HF*H
      H6=HS*H
      XH=X+H
      XH2=X+H2
      CALL SUB(X,Y,W(1,1))
      DO 1 J=1,N
    1 W(J,2)=Y(J)+H2*W(J,1)
      CALL SUB(XH2,W(1,2),W(1,3))
      DO 2 J=1,N
      W(J,1)=W(J,1)+2*W(J,3)
    2 W(J,2)=Y(J)+H2*W(J,3)
      CALL SUB(XH2,W(1,2),W(1,3))
      DO 3 J=1,N
      W(J,1)=W(J,1)+2*W(J,3)
    3 W(J,2)=Y(J)+H*W(J,3)
      CALL SUB(XH,W(1,2),W(1,3))
      DO 4 J=1,N
    4 Y(J)=Y(J)+H6*(W(J,1)+W(J,3))
      X=XH
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1996/04/01 15:02:16 mclareni
	6	* Mathlib gen
	7	*
	8	*
	9	*
	10	* rkstpcod.inc
	11	*
	12	DIMENSION Y(),W(N,)
	13
	14	PARAMETER (Z1 = 1, HF = Z1/2, HS = Z1/6)
	15
	16	IF(N .LE. 0) RETURN
	17	H2=HF*H
	18	H6=HS*H
	19	XH=X+H
	20	XH2=X+H2
	21	CALL SUB(X,Y,W(1,1))
	22	DO 1 J=1,N
	23	1 W(J,2)=Y(J)+H2*W(J,1)
	24	CALL SUB(XH2,W(1,2),W(1,3))
	25	DO 2 J=1,N
	26	W(J,1)=W(J,1)+2*W(J,3)
	27	2 W(J,2)=Y(J)+H2*W(J,3)
	28	CALL SUB(XH2,W(1,2),W(1,3))
	29	DO 3 J=1,N
	30	W(J,1)=W(J,1)+2*W(J,3)
	31	3 W(J,2)=Y(J)+H*W(J,3)
	32	CALL SUB(XH,W(1,2),W(1,3))
	33	DO 4 J=1,N
	34	4 Y(J)=Y(J)+H6*(W(J,1)+W(J,3))
	35	X=XH
	36	RETURN
	37	END