[t, y] = ode45('vdp1', [0 25], [3 3]', odeset('Stat', 'on'));