Thanks, Enrique. Yes, I know the work by Seybold and Hilfer. However, their calculation of the inverse Mittag-Leffler function involves the numerical solution of a non-linear equation (let's say by bisection, Newton, or any similar methods). I need to embed the inverse Mittag-Leffler calculations within an optimisation algorithm and process a massive amount of data, for which the mentioned approach would be terribly inefficient.
Pablo, thanks for the reference too. However, their series only converges when 0
Thanks a lot, Enrique. Actually I'm currently using Roberto's code (available in Matlab FileExchange), particularised and vectorised for calculations on the real axis.