function xp = bksys(t,x,k,b);
xp = -(k/b)*x;