function metgadjanja

    [x,y1]=ode45(@nehomogen, [0 0.5], [2 0]);

    [x,y2]=ode45(@homogen, [0 0.5], [1 2]);

    c=-y1(end,1)/y2(end,1)

end



function yp=nehomogen(x,y)
    yp=[y(2); -y(1)+x];
end

function yp=homogen(x,y)
    yp=[y(2); -y(1)];
end 