Abstract It is well known that certain fractional diffusion equations can be solved by the densities of stable Lévy motions. In this paper we use the classical semigroup approach for Lévy processes to define semi-fractional derivatives, which allows us to generalize this statement to semistable Lévy processes. A Fourier series approach for the periodic part of the corresponding Lévy exponents enables us to represent semi-fractional derivatives by a Grünwald-Letnikov type formula. We use this formula to calculate semi-fractional derivatives and solutions to semi-fractional diffusion equations numerically. In particular, by means of the Grünwald-Letnikov type formula we provide a numerical algorithm to compute semistable densities.