x0=linspace(0,1/4,20);
y0=-4*x0+1;

t=0:0.1:5;
fix=zeros(length(t),20);
fiy=zeros(length(t),20);

for i=1:20
    [t fi]=ode45('gyak4b',t,[x0(i),y0(i)]);
    fix(:,i)=fi(:,1);
    fiy(:,i)=fi(:,2);
    hold on;
    plot(fi(:,1),fi(:,2),'b');
end

%plot(fix,fiy);
xlim([-0.5 3]);
ylim([-0.5 3]);

x1=ones(1,20)*2.8;
y1=linspace(0.005,0.007,20);

t=0:-0.1:-5;
fix=zeros(length(t),20);
fiy=zeros(length(t),20);

for i=1:20
    [t fi]=ode45('gyak4b',t,[x1(i),y1(i)]);
    %fix(:,i)=fi(:,1);
    %fiy(:,i)=fi(:,2);
    hold on;
    plot(fi(:,1),fi(:,2),'r');
end



%{
function dydt=gyak4(t,y)
dydt=zeros(2,1);
dydt(1)=y(1)*(3-2*y(2)-y(1));
dydt(2)=y(2)*(1-2*y(1)-y(2)/2);
end
%}