Finite volume simulation of a semi-linear Neumann problem (Keller-Segel model) on rectangular domains

In this study, the finite volume method is implemented for solving the problem of the semilinear equation: $-d \delta u+ u=u^q (d, q>0)$with a homogeneous Neumann boundary condition. This problem is equivalent to the known stationary Keller-Segel model, which arises in chemotaxis.After discretization, a nonlinear algebraic system is obtained and solved on the platform Matlab. As a result, many single peaked and multi-peaked shapes in 3D and contour plots can be drawn depending on the parameters d and q.