function y1 = rk2step(fun, h, t, y)
% RK2STEP one integration step of RK2
	k1 = h*fun(t, y);
	k2 = h*fun(t + h, y + k1);
    y1 = y + (k1 + k2)/2;
end
