A New Collocation Algorithm for Solving Even-Order Boundary Value Problems via a Novel Matrix Method

This paper is dedicated to presenting and analyzing a numerical algorithm for the solution of even-order boundary value problems. The proposed solutions are spectral and they depend on introducing a new matrix of derivatives of certain shifted Legendre polynomial basis, along with the application of the collocation method. The nonzero elements of the introduced matrix are expressed in terms of the well-known harmonic numbers. Numerical examples provide favorable comparisons with other existing methods and ascertain the efficiency and applicability of the proposed algorithm.