“…Many applied problems with memory can be successfully modeled via a fractional system of differential and/or integral equations, for example, semi-conductor devices [24], population dynamics [25], and identification of memory kernels in heat conduction [26]. During last years, several numerical techniques have been applied for solving fractional systems of differential and integral equations, for example, fractional power Jacobi spectral method [27], Bernoulli wavelets method [28], Haar wavelets method [29], Müntz-Legendre wavelets method [30], Block pulse functions method [31], finite difference method [32], spline collocation method [33][34][35], and spectral method [36].…”