function [tout,yout] = RK4(F,tspan,y0,h);
    t0 = tspan(1);
    tfinal = tspan(2);
    t = t0;
    y = y0(:);
    tout = t;
    yout = y.';
    while t < tfinal
        s1 = F(t,y);
        s2 = F(t + 0.5*h, y + h*(0.5*s1));
        s3 = F(t + 0.5*h, y + h*(0.5*s2));
        s4 = F(t + h, y + h*s3);
        y = y + h * (s1/6 + s2/3 + s3/3 + s4/6);
        t = t + h;
        tout(end+1,1) = t;
        yout(end+1,:) = y.';
    end
end
% Simulaci jsem nalezl maximalni delku kroku
%pro zachovani staility hmax = 27.852935