added efficient solve_ode working with harmonic oscillator for variable parameters