“…Moreover numerical solutions of conservation equations, with variable coefficients in the heterogeneous media, need to be matured. In order to core calculation in the nuclear reactors there is a new version of the neutron diffusion equation which is established on the fractional space derivatives (Maleki Moghaddam et al, 2014;Espinosa-Paredes et al, 2013).The neutron fractional diffusion equation, NFDE, uses the non-local gradient of Fick's law (J = ÀD a r a /) where the fractional derivative operator, r a , can be defined in the Riemann-Liouville, Caputo and Grünwald-Letnikov sense (Das, 2012;Oldham and Spanier, 1974;Hilfer, 2000;Miller and Ross, 1993) and a is the order of fractional derivatives. In the limit a ?…”