Abstract This paper is devoted to the numerical analysis of non-smooth ensemble optimal control problems governed by Liouville (continuity) equation that have been originally proposed R.W. Brockett with purpose determining an efficient and robust strategy for dynamical systems. A methodology solving these presented based on a Lagrange optimization framework where controls are characterized as s...