Start with the equations of motion
$$ m \ddot{x} = -k x - d \dot{x} \tag{1}$$
where $m$ is mass, $k$ is stiffness (force/distance) and $d$ is damping (force/speed).
Parametrize the problem with $k = m\, \omega_n^2$ and $d =2 \zeta m\, \omega_n$ where instead of specifying $m$, $k$ and $d$, you specify a system with $\omega_n$ undamped natural frequency and $\zeta$ damping ratio.
The equations of motion are
$$ \ddot{x} = -\omega_n^2\, x - 2 \zeta \omega_n\, \dot{x} \tag{2}$$
with the well known solution of the type
$$ x= X \exp\left(-\lambda t\right) \sin \left( \omega t \right)\tag{3}$$
with $\lambda = \zeta \omega_n$ and $\omega = \omega_n \sqrt{1-\zeta^2}$.
Expanded this is
$$ x= X \exp\left(-\tfrac{d}{2 m} t\right) \sin \left( t\sqrt{\tfrac{k}{m}-\tfrac{d^2}{4 m^2}} \right)\tag{4}$$
provided that $ d \leq 2 \sqrt{k m}$.
The solution for higher damping is
$$ x= X \exp\left(-\tfrac{d}{2 m} t\right) \sinh \left( t\sqrt{\tfrac{d^2}{4 m^2}-\tfrac{k}{m}} \right)\tag{6}$$