% 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).