“…Here the power-like shape of the memory functions, sitting as kernels in the Caputo derivatives, turns out to be decisive to judge nonnegativity of obtained solutions and show the stable character of the subordination. Among another paths to generalize diffusion equations we mention giving up the constant character of the diffusion coefficient [67] or modeling diffusive processes using evolution equations with more than one fractional time derivative, including also the case of so-called distributed order fractional derivatives [20,21,37]. Models using the time dependent diffusion coefficient, inspired by experimental data, have been used to explain differences in behavior of MSDs for short and long times -relevant examples are provided by the scaled Brownian motion with diffusion coefficient D(t) ∝ t α−1 , where t ∈ R + and α > 0, for which the MSD grows linearly at short times, x 2 (t) ∝ t, whereas at long times behaves as x 2 (t) ∝ t α [17,42,68,69,81].…”