function dy = Osc1(t,y);
w = 4.*pi.^2;
dy = [y(2); -w.*y(1)];