The aim of this paper is to establish convergence, properties and error bounds for the fully discrete solutions of a class of nonlinear systems of reaction-diffusion nonlocal type with moving boundaries, using the finite element method with polynomial approximations of any degree. A coordinate transformation which fixes the boundaries is used. Some numerical tests to compare our Matlab code with a moving finite element method are investigated.