x=c_1 e^{-\left(\lambda-\sqrt{\lambda^2-\omega_0^2}\right)t}
+ c_2 e^{-\left(\lambda+\sqrt{\lambda^2-\omega_0^2}\right)t}
