The theory of differential inclusions has played a central role in many areas as biological systems, physical problems and population dynamics. The principle aim of our work is to compute explicitly the discrete approximate solution of a differential inclusion including a maximal monotone operator. Also we present a numerical application of our results for showing how to compute the discrete approximate solution of its corresponding differential inclusion.

In this paper, we present a method for finding the zero of the sum of finite family of maximal monotone operators on real Hilbert spaces. In the case where the number of maximal monotone operators is three, we define a function such that its fixed points are solutions of our problem. Some illustrative examples are given at the end of this paper.

Our aim is to calculate the discrete approximate solution of di⁄erential inclusion with normal cone and prox-regular set, the question is how to calculate this solution? We use the discrete approximation property of a new variant of nonconvex sweeping processes involving normal cone and a nite element method. Knowing that The majority of mathematicians have proved only the existence and uniqueness of the solution for this type of inclusions, like: Mordukhovich, Thibault, Aubin, Messaoud,
...etc.