Pseudo-spectral matrices as a numerical tool for dealing BVPs, based on Legendre polynomials’ derivatives