We consider the numerical integration of coupled self-adjoint non-autonomous partial differential systems. Under convergence conditions, the solution can be written as a series expansion where each of its terms correspond ...