Calculate collocation coefficients using long double to avoid cancellation errors, #948
IRKIntegrator: added construction of discrete forward dynamics, #948