Based on the homogeneous balance method, the Jacobi elliptic expansion method and the auxiliary equation method, the first elliptic function equation is used to get a new kind of solutions of nonlinear evolution equations. New exact solutions to the Jacobi elliptic function of MKdV equations and Benjamin-Bona-Mahoney (BBM) equations are obtained with the aid of computer algebraic system Maple. The method is also valid for other (1+1)-dimensional and higher dimensional systems.Nonlinear partial differential equations (NPDEs) are widely used to describe phenomena in physics, biology, chemistry, and so on. Some of the most interesting features of physical systems are hidden in their nonlinear behavior, and can only be studied with appropriate methods. Over the past decades, many authors have mainly paid attention to construct various methods such as Bäcklund transformation (see [1,5,10,11,36]), Darboux transformation [36] , Inverse scattering method [10] , Hirota's bilinear method [11] , the tanh-function method [24] , the sine-cosine method [38], the homogeneous balance method [37] , symmetry reduction method (see [2, 19, 21-23, 26, 27, 29]), the geometrical method [6] , the formal variable separation approach [4] , the multilinear variable separation approach (MLVSA) [17,20] , the functional variable separation approach (FVSA) (see [30,31,39]) and the derivative-dependent functional variable separation approach (DDFVSA) [40] . Recently, an extended tanh-function method [7] and symbolic computation are used in [8] for solving the new coupled modified KdV equations to obtain four kinds of soliton solutions.The aim of this paper is to use a direct method to solve two different types of equations such as the MKdV and BBM equations.