Lie’s theory of symmetry groups plays an important role in analyzing and solving differential equations; for instance, by decreasing the order of equation. Moreover, there are some analytic methods to find the infinitesimal generators that span the Lie algebra of symmetries. In this paper, we first converted the problem of finding infinitesimal generators in to the problem of solving a system of polynomial equations in the context of computational algebraic geometry. Then, we used Gröbner basis a novel computational tool to solve this problem. As far as we know, when a differential equation contains some parameters, there is no linear algebraic algorithm up to our knowledge to deal with these parameters; so, we must apply the algorithms, which are based on Gröbner basis.