Speaking of 2nd-order linear ODEs w/ const factor (a.k.a. mass-spring-damper systems), I wrote a blog post[0] deriving all possible general solutions in a concise matrix that makes it easily implementable in code. The following is the complete solution in Lua: function sprung_response(t,pos,vel,k,c,m) local decay = c/2/m local omega = math.sqrt(k/m) local resid = decay*decay-omega*omega local scale = math.sqrt(math.a…
My grad school was basically working through F=Ma-Cv-Kx in all it's expanding varied glory, working up to full blown finite element analysis. Ultimately it all comes down to choosing the most convenient basis functions for the questions you're answering.
If I'm understanding correctly, you want another = there for a standard damped spring: F = ma = -Cv - Kx.