This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / d / rkstpcod.inc
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