Am trying to compute the analytic solution of a boundary value problem for the time fractional inverse advection equation. Some where along the way when inserting the boundary conditions, I encounter an inversion of a Mittag-Leffler function. The function I have is E_{\alpha,1}(-\pho^2t^\alpha). At the same time I also have to compute the product between two Mittag-functions.

Please help me am stack