function dot = simplependulum(t,input)

global m J lc g B

hdot = -B*sign(input(1))-m*g*lc*sin(input(2));
thetadot = input(1)/J;

dot = [hdot;thetadot];