“…The above expression can, in turn, be written as Eu n =Āu n +Bc n + F (u n ) (4) where F : R n → R n represents the nonlinear term. The stiffness matrix Ā ∈ M n×n , the control matrix B ∈ M n×m c , as well as the mass matrix E ∈ M n×n , are assembled as usual, by computing the element integrals with Gauss-Legendre formulas.…”