You can reduce any (linear) n-th order ordinary differential equation into n coupled first order ODEs simply enough:
Your third order ODE is:
a.x''' + b.x'' + c.x' + dx + e = 0
Introduce variables q, r:
y = x'
z = x'' (= y')
Then your equation becomes:
a.z' + b.z + c.y + d.x + e = 0
y - x' = 0
z - y' = 0
... three coupled first-order ODEs

You can solve these problems easy enough.
Incidentally, this procedure is the basis of solving partial differential equations using (traditional) numerical methods. You break down the PDE into a massive system of linearlised ODEs, and solve using Runge-Kutta or similar.