Computation of Green's functions through algebraic decomposition of operators

In this article we use linear algebra to improve the computational time for the obtaining of Green’s functions of linear differential equations with reflection (DER). This is achieved by decomposing both the `reduced’ equation (the ODE associated to a given DER) and the corresponding two-point boundary conditions.