v(t)={1 \over dt} \cdot [{- x_0 e^{-\beta t} \sin (\omega t+\alpha) \omega dt - \beta dt x_0 e^{-\beta t} \cos (\omega t+\alpha)]