Numerical solution of Fractional PDE using(OHAM Method)

i want expert mathmatica or matlab to make a code for 2 example in fractional pde using OHAM method. i want to get optimal values (c1,c2) and figuers the solutions using different values of (alpha)and exact solution in one figures and get the value of error and make comparison with exact solutions to prove the convergence of OHAM methodas the same in the attached PDF
Hint: (it is not my project it is for example what i want)