Abstract: A variational iteration method (VIM) for a class of fractional predator-prey model is studied. A pair equivalent coupled integro-differential equation to the model is obtained by means of integral transformation. According to variational theory, the Lagrange multiplier is calculated and the VIM scheme is constructed. Finally, a convergence analysis of the scheme is given and a numerical simulation is carried out.The results verify the feasibility and effectiveness of the method.