t0=[0 400];
for i=1:10
    x0=[7*rand 7*rand];
     %x0=[rand rand];
    [t,p]=ode45('gyak4b', t0,x0);
    plot(p(:,1), p(:,2));
    hold on
end
%x0=[0.5 6];
%[t,p]=ode45('gyak4b', t0,x0);
%plot(p(:,1), p(:,2));
xlim([0 7]);
ylim([0 7]);
