A sequential design for approximating the Pareto front using the expected Pareto improvement function