  We provide a new approach to obtain solutions of certain evolution equations
set in a Banach space and equipped with nonlocal boundary conditions. From this
approach we derive a family of numerical schemes for the approximation of the
solutions. We show by numerical tests that these schemes are numerically robust
and computationally efficient.
