%hazis
function dxydt=gyak4b(t, p)
dxydt=zeros(2, 1);
dxydt(1)=3*p(1)-2*p(1)*p(2)-p(1)^2;
dxydt(2)=p(2)-2*p(1)*p(2)-0.5*p(2)^2;
