problem comes from "Numerical solution of fractional integro-differential , equations by collocation method , E.A. Rawashdeh, Department of Mathematics, Yarmouk University, Irbid 21110, Jordan"
$D^qy(t)=p(t)y(t)+f(t)+\int_{0}^{1}{K(t,s)y(s)\,ds} , t\in I=[0,1]$
I want to create a maple code to check if the results in given article is valid or not but I do not have any idea about collocation method!
any reference to collocation method solution are welcome!
See first how the collocation method is applied to simpler problems and then advance with your problem. See for example here where the method is applied to find the solution of a an ODE with boundary conditions. Here is paper on collocation method for solving integral equations.