Evaluates an interpolated Runge-Kutta result 'rk_res' at 't'.
For details and examples, c.f. 'odeRK()'.