function dzdt = ode3(t,z);
% ODE3: computes dzdt. 
w=1;
dzdt = [z(2); -w^2*sin(z(1))];