A numerical solution of the modified Korteweg-de Vries (MKdV) equation is presented by using a nonstandard finite difference (NSFD) scheme with theta method which includes the implicit Euler and a Crank-Nicolson type discretization. Local truncation error of the NSFD scheme and linear stability analysis are discussed. To test the accuracy and efficiency of the method, some numerical examples are given. The numerical results of NSFD scheme are compared with the exact solution and a standard finite difference scheme. The numerical results illustrate that the NSFD scheme is a robust numerical tool for the numerical integration of the MKdV equation.