Numerical simulation of the modified Korteweg-de Vries equation