EP matrices are a wide class of matrices which, among other things, can be characterized through factorizations. In this paper we are using two factorization algorithms in order to compute and factorize the Moore-Penrose inverse of a singular EP matrix. For the implementation of the algorithms we make use of a Computer Algebra System such as Maple. The results given by the proposed algorithms are faster and more accurate than the built-in Maple function in both symbolic and numerical tests, using matrices of different dimensions.