I am working on numerical solutions of some nonlinear reaction diffusion Riesz space fractional equations. I wondering if there is a better way to code the coefficient matrix obtained from the discretization of the diffusion term (finite difference). I am aware of the m.file provided by Podlubny et at but there is no much to modify when using that.

Any suggestions or link to sample code would be highly appreciated.