subroutine timeloop2(y, n, maxsteps, step, time, nsteps)

integer n, step, nsteps, maxsteps
real*8 time, y(n,0:maxsteps-1)

y  : solution (all time steps), input and output
n  : no of solution components (2 in our example), input
maxsteps : max no of time steps, input
step   : no of current time step, input and output
time   : current value of time, input and output
nsteps : no of time steps to advance the solution