% function p2c_tst0 % needs to write in the workspace % Feb 2016, J. Gaspar A= [0 1; -1 -1]; % B= [0; 1]; % C= eye(2); % D= [0; 0]; th= (0:45:360)*pi/180; x0= [cos(th); sin(th)]; % collection of initial conditions x1= A*x0; % compute local evolution (direction) % plot the results quiver(x0(1,:), x0(2,:), x1(1,:), x1(2,:)) axis equal grid on