function [t,y]=RK4_1Step(t,y,f,h,params)

K1 = f(t,y,params);
K2 = f(t+0.5*h,y+0.5*h*K1,params);
K3 = f(t+0.5*h,y+0.5*h*K2,params);
K4 = f(t+h,y+h*K3,params);

y = y + (h/6)*(K1+2*K2+2*K3+K4);

t = t + h;

end