SUBROUTINE RKDUMB(VSTART,NVAR,X1,X2,NSTEP) PARAMETER (nmax=100000) real*8 vstart,v,dv,xx,y,x1,x2,x,h COMMON /PATH/ XX(500),Y(nmax,500) DIMENSION VSTART(NVAR),V(NMAX),DV(NMAX) DO 11 I=1,NVAR V(I)=VSTART(I) Y(I,1)=V(I) 11 CONTINUE XX(1)=X1 X=X1 H=(X2-X1)/NSTEP DO 13 K=1,NSTEP CALL DERIVS(nvar,X,V,DV) CALL RK4(V,DV,NVAR,X,H,V) IF(X+H.EQ.X) write(6,*) 'Stepsize not significant in RKDUMB.' X=X+H XX(K+1)=X DO 12 I=1,NVAR Y(I,K+1)=V(I) vstart(i) = v(i) 12 CONTINUE 13 CONTINUE RETURN END