This paper studies the integration of inclusion of subdifferentials. Under various verifiable conditions, we obtain that if two proper lower semicontinuous functions $f$ and $g$ have the subdifferential of $f$ included in the $\gamma$-enlargement of the subdifferential of $g$, then the difference of those functions is $ \gamma$-Lipschitz over their effective domain.

Formulas for the Clarke subdifferential are always expressed in the form of inclusion. The equality form in these formulas generally requires the functions to be directionally regular. This paper studies the directional regularity of the general class of extended-real-valued functions that are directionally Lipschitzian. Connections with the concept of subdifferential regularity are also established.

We construct Lipschitz functions such that for all $s>0$ they are $s$-H\"older, and so proximally, subdifferentiable only on dyadic rationals and nowhere else. As applications we construct Lipschitz functions with prescribed H\"older and approximate subderivatives.

Given an integral functional defined on $L_p$, $1 \leq p <\infty$, under a growth condition we give an upper bound of the Clarke directional derivative and we obtain a nice inclusion between the Clarke subdifferential of the integral functional and the set of selections of the subdifferential of the integrand.

Recently, F.~H.~Clarke and Y.~Ledyaev established a multidirectional mean value theorem applicable to lower semi-continuous functions on Hilbert spaces, a result which turns out to be useful in many applications. We develop a variant of the result applicable to locally Lipschitz functions on certain Banach spaces, namely those that admit a ${\cal C}^1$-Lipschitz continuous bump function.