It is well known that the discretization of fractional diffusion equations with fractional derivatives
, using the so‐called weighted and shifted Grünwald formula, leads to linear systems whose coefficient matrices show a Toeplitz‐like structure. More precisely, in the case of variable coefficients, the related matrix sequences belong to the so‐called generalized...