
% Dušan Kokosík, 317 484, 2. ročník Bc štúdia, 3ZEH20
function semestralna_praca_3_dusan_kokosik_317484(f,x,y,a,b,h)

% !!!!! AK SA ZADÁVA SÚSTAVA DIFERENCIÁLNYCH ROVNÍC !!!!!
% na vstup používať stĺpcové vektory pre y,f typu: y = [0;1;...]
% f = @(x,y) [y(1)+y(2);-y(1)+y(2)]

interval = a:h:b;
for i = 1:length(interval)-1
    k1 = f(x,y)
    k2 = f(x+0.5*h, y+0.5*h*k1)
    k3 = f(x+0.5*h, y+0.5*h*k2)
    k4 = f(x+h, y+h*k3)

    x = h*i
    y = y + h/6*(k1 + 2*k2 + 2*k3 + k4)
    %disp(x);
    %disp(y);
end

disp(x)
disp(y);