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