Collocation method using quintic B-splines finite element have been developed for solving numerically the Hirota-Satsuma coupled MKdV equation. Accuracy of the proposed method is shown numerically by calculating conservation laws, L_2 and L_∞ norms on studying of a soliton solution. It is shown that the collocation scheme for solutions of the MKdV equation gives rise to smaller errors and is quite easy to implement. Numerical experiments support these theoretical results.