In this paper, we propose a numerical method to solve the classic $L^2$-optimal transport problem. Our algorithm is based on use of multiple shooting, in combination with continuation procedure, boundary value problem associated Based viewpoint Wasserstein Hamiltonian flow initial and target densities, our reflects structure underlying exploits it discretization. Several examples are presented ...