With the aid of symbolic computation system Maple and by using the extended hyperbolic function method, the exact solitary wave solutions of Zakharov equations are obtained which include bell-shaped soliton solutions, kink-shaped soliton solutions, envelop soliton solutions, singular soliton solutions and a new type of soliton solutions. The method can be also used to solve other nonlinear wave equations.