function x_dot = f_onespring(t,x)
%simple mass-spring damper: 
global m b k l0 f_ext;