Apply variation of parameters to the underdamped oscillator \ddot x + 2\dot x + 5x = e^(-t)sin t