% Example: Node-Sink

clear all
tfinal=100;
[t xa]=ode23(@nodesink,[0 tfinal],[-3 4]);
[t xb]=ode23(@nodesink,[0 tfinal],[-2 4]);
[t xc]=ode23(@nodesink,[0 tfinal],[-1 4]);
[t xd]=ode23(@nodesink,[0 tfinal],[1 4]);
[t xe]=ode23(@nodesink,[0 tfinal],[2 4]);
[t xf]=ode23(@nodesink,[0 tfinal],[3 4]);

[t xg]=ode23(@nodesink,[0 tfinal],[-3 -4]);
[t xh]=ode23(@nodesink,[0 tfinal],[-2 -4]);
[t xi]=ode23(@nodesink,[0 tfinal],[-1 -4]);
[t xj]=ode23(@nodesink,[0 tfinal],[1 -4]);
[t xk]=ode23(@nodesink,[0 tfinal],[2 -4]);
[t xl]=ode23(@nodesink,[0 tfinal],[3 -4]);


plot(xa(:,1),xa(:,2),'b');
hold on
plot(xb(:,1),xb(:,2),'b');
hold on
plot(xc(:,1),xc(:,2),'b');
hold on
plot(xd(:,1),xd(:,2),'b');
hold on
plot(xe(:,1),xe(:,2),'b');
hold on
plot(xf(:,1),xf(:,2),'b');
hold on
plot(xg(:,1),xg(:,2),'b');
hold on
plot(xh(:,1),xh(:,2),'b');
hold on
plot(xi(:,1),xi(:,2),'b');
hold on
plot(xj(:,1),xj(:,2),'b');
hold on
plot(xk(:,1),xk(:,2),'b');
hold on
plot(xl(:,1),xl(:,2),'b');

xlabel('X1');
ylabel('X2');
title('Node Sink')
grid


% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009
% Copyright(c).