No abstract

Physical laws governing population dynamics are generally expressed as differential equations. Research in recent decades has incorporated fractional-order (non-integer) derivatives into differential models of natural phenomena, such as reaction-diffusion systems. In this paper, we develop a method to numerically solve a multi-component and multi-dimensional space-fractional system. For space discretization, we apply a Fourier spectral method that is suited for multidimensional partial differential equation systems. Efficient approximation of time-stepping is accomplished with a locally one dimensional exponential time differencing approach. We show the effect of different fractional parameters on growth models and consider the convergence, stability, and uniqueness of solutions, as well as the biological interpretation of parameters and boundary conditions.