[t, y, s] = ode45('vdp1', [0 25], [3 3]');
s