function [DF,F]=rk_func(sw,t_anf,t_end,model,plott)
t=t_anf;
h=0.0001;
x=sw;
A=eye(length(sw));
t1=t_end-h/2;
if plott==1
    lo1(1)=x(1);
    lo2(1)=x(2);
    i=1;
    T(1)=t_anf;
end
while t<t1
    if model==1
        k1=lin_rw(t,x);
        k2=lin_rw(t+h/2,x+(h/2)*k1);
        k3=lin_rw(t+h/2,x+(h/2)*k2);
        k4=lin_rw(t+h,x+h*k3);
        K1=lin_rw_mat(t,A,x);
        K2=lin_rw_mat(t+h/2,A+(h/2)*K1,x);
        K3=lin_rw_mat(t+h/2,A+(h/2)*K2,x);
        K4=lin_rw_mat(t+h,A+h*K3,x);
    end
    if model==2
        k1=harm_osc(t,x);
        k2=harm_osc(t+h/2,x+(h/2)*k1);
        k3=harm_osc(t+h/2,x+(h/2)*k2);
        k4=harm_osc(t+h,x+h*k3);
        K1=harm_osc_mat(t,A,x);
        K2=harm_osc_mat(t+h/2,A+(h/2)*K1,x);
        K3=harm_osc_mat(t+h/2,A+(h/2)*K2,x);
        K4=harm_osc_mat(t+h,A+h*K3,x);
    end
    if model==3
        k1=nl(t,x);
        k2=nl(t+h/2,x+(h/2)*k1);
        k3=nl(t+h/2,x+(h/2)*k2);
        k4=nl(t+h,x+h*k3);
        K1=nl_mat(t,A,x);
        K2=nl_mat(t+h/2,A+(h/2)*K1,x);
        K3=nl_mat(t+h/2,A+(h/2)*K2,x);
        K4=nl_mat(t+h,A+h*K3,x);
    end
    if model==4
        k1=nl2(t,x);
        k2=nl2(t+h/2,x+(h/2)*k1);
        k3=nl2(t+h/2,x+(h/2)*k2);
        k4=nl2(t+h,x+h*k3);
        K1=nl2_mat(t,A,x);
        K2=nl2_mat(t+h/2,A+(h/2)*K1,x);
        K3=nl2_mat(t+h/2,A+(h/2)*K2,x);
        K4=nl2_mat(t+h,A+h*K3,x);
    end
    x=x+(h/6)*(k1+2*k2+2*k3+k4);
    A=A+(h/6)*(K1+2*K2+2*K3+K4);
    t=t+h;
    if plott==1
        i=i+1;
        lo1(i)=x(1);
        lo2(i)=x(2);
        T(i)=t;
    end
end
if plott==1
    plot(T,lo1);
    figure
    plot(lo1,lo2)
end
F=x;
DF=A;