function dotq = damped_vibration_global (t, q)
    global b;
    x = q(1);
    v = q(2);
    dotx = v;
    dotv = -b*v - 9*x;
    dotq = [dotx; dotv];
end
