An alternative polynomial expansion for electromagnetic field estimation inside three-dimensional dielectric scatterers is presented in this article. In a continuation with the previous work of authors, the Tensor-Volume Integral Equation (TVIE) is solved by using the Galerkin-based moment method (MoM) consisting of a combination of entire-domain and sub-domain basis functions including three-dimensional polynomials. Instead of using trivial power polynomials, Legendre polynomials are adopted for electromagnetic fields expansion in this study. They have the advantage of being a set of orthogonal functions, which allows the use of high-order basis functions without introducing an ill-condition MoM matrix. The accuracy of such approach in MoM is verified by comparing its numerical results with that of exact analytical method such as Mie theory and conventional procedures in MoM. Besides, it is also confirmed that the condition number of the MoM matrix obtained with the proposed approach is lower than that of the previous approaches.