function y=rungekutta4step(t,x,h)
% one step of RK4 method
k1 = ydot(t, x);
k2 = ydot(t + h/2, x + h/2 * k1);
k3 = ydot(t + h/2, x + h/2 * k2);
k4 = ydot(t + h, x + h * k3);
y = x + h/6 * (k1 + 2*k2 + 2*k3 + k4);