The Skorokhod embedding problem aims to represent a given probability measure on the real line as the distribution of Brownian motion stopped at a chosen stopping time. In this paper, we consider an extension of the optimal Skorokhod embedding problem in Beiglböck, Cox & Huesmann [1] to the case of finitely-many marginal constraints1. Using the classical convex duality approach together with th...