“…In our work, we will present a numerical method for calculating the approximate solution of a differential inclusion with normal cone for prox-regular sets, this problem has been studied extensively see, [18], [19], [17], for existence and uniqueness of the solutions see, [11], [12], [3], [9] and also gave some properties and theories that help us reach the proof of the proposed convergence on our side. We use also the same method (Discrete approximation solution), see, [14], [15], [7], [22], [20], [4], knowing that these methods were applied to the differential inclusion ẋ(t) ∈ F (x(t)) a.e t ∈ [0, T ] , x (0) = x 0 that was studied by Mordukhovich, see, [13] and also applied to the differential inclusion ẋ(t) ∈ −N C (x(t)) a.e t ∈ [0, T ] , x (0) = x 0 see, [16]. In this paper we will study the following problem.…”