A new modification of homotopy analysis method (HAM) is proposed in this paper. The auxiliary differential operator is specifically chosen so that more than one secular term must be eliminated. The proposed method can capture asymmetric and period-2 solutions with satisfactory accuracy and hence can be used to predict symmetry-breaking and period-doubling bifurcation points. The variation of accuracy is investigated when different number of frequencies are considered.