Uncertainty quantification of coupled odes with stochastic galerkin methods using adaptive higher order Runge–Kutta methods