SUBROUTINE RUKUTT(X,Y,H) C FUEHRT EINEN RUNGE-KUTTA-SCHRITT AUS DIMENSION Z(2),Y(2),ST(5),R(2),E(2) COMMON Z E(1)=0 E(2)=0 R(1)=Y(1) R(2)=Y(2) ST(1)=H/2. ST(2)=H/2. ST(5)=H/2. ST(3)=H ST(4)=H D=X DO 2 I=1,4 CALL FKT(D,R) D=X+ST(I) DO 2 K=1,2 R(K)=Y(K)+ST(I)*Z(K) 2 E(K)=E(K)+ST(I+1)*Z(K)/3. X=D Y(1)=Y(1)+E(1) Y(2)=Y(2)+E(2) RETURN END