[T Y] = ode45('Osc1',[0 20.*pi],[0; 1]);
plot(T,Y);