MATLAB numerical or analytical code for fractional ADE

I have the following ADE with caputo Time fractional derivative. I would like to obtain the concentration C, but i dont have any knowlege about the implementation of fractional derivative. But i know how i can solve and draw the curves without fractional derivative. If somebody can help me with a MATLAB code of analytical or numerical solution of the problem. Thé problem is explain in the attached pdf