function rown()

opcje=odeset('OutputFcn',@odephas2);
[t,y]=ode23(@dp,[0 12],[2 4],opcje);
hold on;
[t,y]=ode23(@dp,[0 12],[-4 0],opcje);
hold on;
 	[t,y]=ode23(@dp,[0 12],[-1 -4],opcje);
hold on;
 	[t,y]=ode23(@dp,[0 12],[4 -1],opcje);
set(gca,'XLim',[-5,5],'YLim',[-5,5],'XGrid','on','YGrid','on')
hold on;

function dydt = dp(t,y)
dydt = zeros(2,1);
dydt(1) = [(-1.5*y(1))+y(2)];
dydt(2) = [(0.25*y(1))-(1.5*y(2))];