clear all;
close all;

x0 =[0 -1 1 0];
[t, x] = ode45('twoOsc',[0 100], x0);

figure;
subplot(2,1,1);
plot(t,x(:,1));
subplot(2,1,2);
plot(t,x(:,3));
figure;
plot(x(:,1),x(:,3))

