In this paper we use the theory of Faber polynomials for solving N-dimensional linear initial value
problems. In particular, we use Faber polynomials to approximate the evolution operator creating the
so-called exponential integrators. We also provide a consistence and convergence analysis. Some tests
where we compare our methods with some Krylov exponential integrators are finally shown.
