A Weierstrass elliptic function expansion method and its algorithm are developed in this paper. The method changes the problem solving nonlinear evolution equations into another one solving the correspondingsystem of nonlinear algebraic equations. With the aid of symbolic computation (e.g. Maple), the method is applied to the combined KdV-mKdV equation and (2+1)-dimensional coupled Davey-Stewartson equation. As a consequence, many new types of doubly periodic solutions are obtained in terms of the Weierstrass elliptic function. Jacobi elliptic function solutions and solitary wave solutions are also given as simple limits of doubly periodic solutions.