function linearODE %linearODE.m allows for the visualization of solutions to xdot=A*x %Specify A and the starting position x0 A=[1.24 -1.45 3.80 -1.22]; x0=[1;0]; L=5; %Trajectory length parameter n=500; %finesse of the animation [t,x]=ode45(@der,linspace(0,L,n),x0,[],A); fig=figure; set(fig,'color',[1 1 1]) b=line('xdata',[],'ydata',[],'color',[0 0 1]); a1=line('xdata',[min(x(:,1)) max(x(:,1))],'ydata',[0 0],'color',[0 0 0]) a2=line('xdata',[0 0],'ydata',[min(x(:,2)) max(x(:,2))],'color',[0 0 0]) xlim([min(x(:,1)) max(x(:,1))]) ylim([min(x(:,2)) max(x(:,2))]) grid on for k=1:n set(b,'xdata',x(1:k,1),'ydata',x(1:k,2)); pause(0.01) end function xp=der(t,x,A) xp=A*x;