“The method of lines coupled with BVMs for the numerical approximation of elliptic PDEs”