global m b k; 
m=1;
b=1;
k=10;
x0=[2;1];
tspan = [0:0.01:5];
[t,x] = my_ode('f2',tspan,x0);