Try this: dsl = DSolveValue[{0 == -Fmax (1 - Cos[t]) + k x[t] + d Derivative[1][x][t] + m Derivative[2][x][t], x[0] == xStart, x'[0] == vStart}, x[t], t]; dsl1 = dsl /. d -> 2*Sqrt[m*k] + \[Epsilon] // Simplify; lim=Limit[dsl1, \[Epsilon] -> 0] yielding the following: (* (1/(k m (k + m)^2))E^(-((k t)/Sqrt[ k m])) (E^((k t)/Sqrt[k m]) Fmax k^2 m - 3 Fmax k m^2 + 2 E^((k t)/Sqrt[k m]) Fmax k m^2 - Fmax m^3 + E^((k t)/Sqrt[k m]) Fmax m^3 - Fmax m^2 Sqrt[k m] t - Fmax (k m)^(3/2) t + k^3 m t vStart + 2 k^2 m^2 t vStart + k m^3 t vStart + k^3 m xStart + 2 k^2 m^2 xStart + k m^3 xStart + k^3 Sqrt[k m] t xStart + 2 k (k m)^(3/2) t xStart + m (k m)^(3/2) t xStart - E^((k t)/Sqrt[k m]) Fmax k (k - m) m Cos[t] - 2 E^((k t)/Sqrt[k m]) Fmax (k m)^(3/2) Sin[t]) *) looking like this Plot[lim /. {Fmax -> 1, m -> 1, k -> 1, xStart -> 0, vStart -> 1}, {t, 0, 5}] [![enter image description here][1]][1] Is it, what you are looking for? [1]: https://i.sstatic.net/DL488.jpg