% initial values
x0 = 1
y0 = 0
% solver
w0 = [x0, y0]
[t, w] = ode45('sys',[0,6],w0);
a= w(:,1);
b= w(:,2);
plot(t,a)
hold on
plot(t,b)
% phase plot
% plot(a,b)