% òúðñèì ðåøåíèå ñ óñëîâèÿ à = õ(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)




