function x_dot = f_3springs(t,x)
%simple masses w/ spring dampers: 
global m b k l0 f_ext;