function f=dy4(t,x)
f=zeros(2,1);
f(1)=x(1)-x(1)*x(2);
f(2)=-x(2)+x(1)*x(2);