function dx=equation(t,x)
m=10; h=1; k=294*10^3; g=9.81;
om=sqrt(k/m);
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=g-om^2*x(1);