**Abstract** : A multi-objective differentiable optimization algorithm had been proposed to solve problems presenting a hierarchy in the cost functions, {fj(x)} (j = 1,...,M ≥ 2; x ∈ Ωa ⊆ Rn). The first cost functions for which j ∈ {1, . . . , m} (1 ≤ m < M ) are considered to be of preponderant importance; they are referred to as the “primary cost functions” and are subject to a “prioritized” treatment, in contrast with the tail ones, for which j ∈ {m+1, . . . , M }, referred to as the “secondary cost functions”. The problem is subject to the nonlinear constraints, ck(x) = 0 (k = 1, . . . , K). The cost functions {fj(x)} and the constraint functions {ck(x)} are all smooth, say C2 (Ωa). The algorithm was first introduced in the case of two disciplines (m = 1, M = 2), and successfully applied to optimum shape design optimization in compressible aerodynamics concurrently with a secondary discipline [5] [9]. More recently, the theory has been enhanced in both framework and established results, and the new developments will be presented elsewhere in details. In short, an initial admissible point x⋆A that is Pareto-optimal with respect to the sole primary cost functions (subject to the constraints) is assumed to be known. Subsequently, a small parameter ε ∈ [0,1] is introduced, and it is established that a continuum of Nash equilibria {x ̄ε} exists for all small enough ε. The continuum originates from x⋆A (x ̄0 = x⋆A). Along the continuum: (i) the Pareto- stationarity condition exactly satisfied by the primary cost functions at x⋆A is degraded by a term O(ε2) only, whereas (ii) the secondary cost functions initially decrease, at least linearly with ε with a negative derivative provided by the theory. Thus, the secondary cost functions are reduced while the primary cost functions are maintained to quasi Pareto-optimality. In this report, we firstly recall the definition of the different steps in the computational Nash-game algorithm assuming the functions all have known first and second derivatives (here without proofs). Then we show how, in the absence of explicitly known derivatives, the continuum of Nash equilibria can be calculated approximately via the construction of quadratic surrogate functions. Numerical examples are provided and commented