  In this article, we deal with the efficient computation of the Wright
function in the cases of interest for the expression of solutions of some
fractional differential equations. The proposed algorithm is based on the
inversion of the Laplace transform of a particular expression of the Wright
function for which we discuss in detail the error analysis. We also present a
code package that implements the algorithm proposed here in different
programming languages. The analysis and implementation are accompanied by an
extensive set of numerical experiments that validate both the theoretical
estimates of the error and the applicability of the proposed method for
representing the solutions of fractional differential equations.
