For the time discretization we use the generalized mid-point rule. This is an implicit method. The time marching scheme can be written in concise form, where the matrices are formed by the component matrices of the previous system of equations.

---