  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.
