% търсим решение с условия а = х(0) , b = y(0),
% като знаем че решенията зависят от t , за 0 < t < 4

[t,y] = ode45(@hlp12,[0,4],['a','b']);  
for k = 1:30
    plot (y(1:k,1),y(1:k,2));  % ????? digit ????
    hold on;
    plot (y(k,1),y(k,2),'+');  % отбелязваме местоположенията на точките
    axis ([0 500 0 500])
    MM(k)= getframe;
end
movie (MM)