We develop computational methods for solving the martingale optimal transport (MOT) problem -a version of the classical optimal transport with an additional martingale constraint on the transport's dynamics. We prove that a general, multi-step multi-dimensional, MOT problem can be approximated through a sequence of linear programming (LP) problems which result from a discretization of the marginal distributions combined with an appropriate relaxation of the martingale condition. We further furnish two generic approaches for discretizing probability distributions, suitable respectively for the cases when we can compute integrals against these distributions or when we can sample from them. These render our main result applicable and lead to an implementable numerical scheme for solving MOT problems. Finally, specializing to the one-step model on real line, we provide an estimate of the convergence rate which, to the best of our knowledge, is the first of its kind in the literature. for insightful discussions and comments. † The second author also gratefully acknowledges support from St John's College, Oxford. MSC 2010 subject classifications: Primary 49M25, 60H99; secondary 90C08. 1 imsart-aap ver.