% Чертаем координатните оси
plot([-0.5,12],[0,0],'k');
hold on
plot([0,0],[-5,20],'k');

% Чертаем графиката на решението
[x1, y1] = ode45('pr1zad1func', -0.5, 12, [4; -4; -1]);
plot(x1(:, 1), y1(:, 1), 'g')
hold on
grid on

% Чертаем графиката на правата y = 2*x + 4
xxx = -0.5:0.1:12;
yyy = 2*xxx-4;
plot(xxx,yyy, 'r');
axis([-0.5,12,-5,20]);
