The main purpose of this paper is to discuss the numerical procedures of optimization with respect to unknown distribution functions which are based on ideas of duality. Dual problem is formulated as minimax type problem without concavity of "inner" problem of maximization. Numerical procedure allowed to avoid difficulties concerning "inner" problem solution is proposed.