Boundary-matrix construction of an ordinary differential Green function ($G=G_0-h^T{M}^{-1}b$)

Let $h_1,\dots ,h_n$ be a basis of homogeneous solutions of a linear ordinary differential operator $L$ of order $n$. Let $B_1,\dots ,B_n$ be independent homogeneous linear boundary conditions, which may be coupled. Suppose a fundamental solution $G_0(x,\xi )$ satisfies $L_x{G}_0=\delta (x-\xi )$. Define the boundary matrix $M_{ij}=B_i{h}_j$ and $b_i(\xi )=B_i{G}_0(\cdot ,\xi )$. If $M$ is invertible, thenThe added terms are homogeneous solutions, so the Dirac delta source is preserved. Applying the boundary operators gives $b-M{M}^{-1}b=0$. Any two such Green functions differ, for fixed source position, by a homogeneous solution with coefficient vector in the null space of $M$; invertibility therefore proves uniqueness. A singular boundary matrix signals a nontrivial homogeneous solution satisfying the endpoint conditions, where a compatibility condition and a nonunique inverse may arise. For $L=D^2$, one convenient fundamental solution is $(x-\xi )H(x-\xi )$.