Computation of Green’s functions for boundary value problems with Mathematica

This paper is devoted to construct an algorithm that allows us to calculate the explicit expression of the Green’s function related to a nnth-order linear ordinary differential equation, with constant coefficients, coupled with two-point linear boundary conditions. We develop this algorithm by making a Mathematica package.