function [t,y] = setequationsimnode
tspan=[0,5];
for k=-10:10;
  for l=-10:20:10; 
    init=[k;l];
    [t,y]=ode45(@f,tspan,init);
    plot(y(:,1),y(:,2));
    hold on;
  end
end

for k=-10:10;
  for l=-10:20:10;
    init=[l;k];
    [t,y]=ode45(@f,tspan,init);
    plot(y(:,1),y(:,2));
    hold on;
  end
end

%Plot annotation 
xlabel('x1')
ylabel('x2')
title('AN IMPROPER NODE');
grid on;

function yprime=f(t,y)
yprime=zeros(2,1);
yprime(1)=-y(1)+y(2);
yprime(2)=-y(2);
end
clc
hold off
end