Underdamped solution to a damped driven oscillator

A damped, driven oscillator is described by the equation,

$$m\frac{d^2x}{dt^2}+b\frac{dx}{dt} + kx = F_0\cos(\omega t).$$

If $b^2-4km < 0$, the system is underdamped and the solution has the form,

$$x(t) = C \exp \left(-t/\tau\right)\cos (\omega_0 t+\delta)+\frac{F_0}{\rho}\cos (\omega t -\theta),$$

where

$$\rho = \sqrt{(k-m\omega^2)^2+\omega^2b^2},\qquad\tan\theta = \frac{\omega b}{k-m\omega^2},\qquad \omega_0 = \sqrt{\frac{k}{m}-\frac{b^2}{4m^2}},\qquad \tau=\frac{2m}{b},$$ $$\tan\delta = -\frac{1}{\omega_0}\frac{v_{x0} +\frac{\omega F_0}{\rho}\sin (-\theta)}{x_0-\frac{F_0}{\rho}\cos (-\theta)} -\frac{b}{2m\omega_0},\qquad C=\frac{x_0-\frac{F_0}{\rho}\cos (-\theta)}{\cos\delta}.$$

Here $x_0$ is the position at $t=0$ and $v_{x0}$ is the velocity at $t=0$.

$x$

$t$

$v_x$

$t$

$m=$ 1 [kg]

$F_0=$ 1 [N]

$b/\sqrt{4km}=$ 0.1 

$\omega=$ 1 [rad/s]

$k=$ 1 [N/m]

$x_0=$ 1 [m]

$v_{x0}=$ 1 [m/s]