Using 4th order Runge-Kutta to approximate a damped harmonic oscillator. Not sure of which function to use. Mathematica

Question

The module I have below has worked for multi order differentials up to this point I believe. However, whenever I try to use some variant of an equation used to approximate t