function y = rk4step(y, t, h, f)
%RK4STEP one integration step of RK4
	k1 = h*f(t, y);
	k2 = h*f(t + h/2, y + k1/2);
	k3 = h*f(t + h/2, y + k2/2);
    k4 = h*f(t + h, y + k3);
    y = y + (k1 + 2*k2 + 2*k3 + k4)/6;
end
