New computational schemes, symbolic-numerical algorithms and programs implementing the high-accuracy finite element method (FEM) for the solution of quantum mechanical boundary-value problems (BVPs) are reviewed. The elliptic BVPs in 2D and 3D domains are solved using the multivariable FEM and Kantorovich method using parametric basis functions. We demonstrate and compare the efficiency of the proposed calculation schemes, algorithms, and software by solving the benchmark BVPs that describe the scattering on a barrier and a well, the bound states of a helium atom, and the quadrupole vibration in a collective nuclear model. © 2018, Allerton Press, Inc.