There is a Generalized Variation of Parameters.
But again it involves finding all 3 independent solution to,
z''+z=0 (after reduction of order).
Are three linearly independent solutions.
Now! You can use the generalized variation of parameters techinique for linear differencial equations of order 3.