function x_dot = f3(t,x)
global m  k;
x1=x(1);
x2=x(2);
x_dot = [x2; -k*x1/m];