a=1;b=.1;[t,x]=ode45(@(t,x)[a b;-b a]*x,[0,10*pi],[1;1]);plot(x(:,1),x(:,2));xlabel('x1');ylabel('x2');
