global b; 
b=1;
v0=1;
tspan = [0:0.01:5];
[t,x] = my_ode('f1',tspan,v0);