Dear Spencer,

How can I implement or modify the EulerCFE code (only 1D is needed) to obtain a variant with mass and stiffness matrices with variable coefficients (i.e. M_A = \int_{-1}^{1} L_q L_p A(\xi,t) d\xi), as well as setting up similar matrices with variable coefficients for setting up source terms (e.g. M_s S_p).

Best regards,

Helmut