“…Various powerful methods for obtaining the exact solutions for the nonlinear partial differential equations have been proposed (see, for example ). Among of these methods are Bäcklund transformation method [12,[31][32][33], the inverse scattering transform method , extended tanh method [9][10][11][12][13], Adomian Pade approximation [28], variational iteration method [5,6,9,17,22,30], homotopy perturbation method [18][19][20]23,40,41], various Lindstedt-Poincare methods [21,22], Adomian decomposition method [1][2][3][4]8,10,11,16,24,26,27,40], F-expansion method [34], exp-function method [35], the Jacobi elliptic function expansion method [37], the complex hyperbolic function method [39] and so on. The objective of this paper is to use the modified variational iteration method (MVIM) to solve the generalized nonlinear Hirota Satsuma coupled KdV equations [14,15,37,38] …”