y^0 &=& (x_1^0, ..., \dot{x_1}^0, ...)^T \\ y^{k+1} &=& y^k + \Delta t (A y^k + b)