function dy=hw4_3_fun(t,Y)
dy=zeros(2,1); %makes output a 2d column vector
dy(1)=Y(1).^2-Y(1)-Y(2);
dy(2)=-Y(2)+Y(1);