  We use a numerical-analytic technique to construct a sequence of successive
approximations to the solution of a system of fractional differential
equations, subject to Dirichlet boundary conditions. We prove the uniform
convergence of the sequence of approximations to a limit function, which is the
unique solution to the boundary value problem under consideration, and give
necessary and sufficient conditions for the existence of solutions. The
obtained theoretical results are confirmed by a model example.
