We consider a nonlinear model that describes the dynamics of a phytoplankton population with aggregation and competition between individual cells. A finite difference method is developed for approximating the solution of this partial differential equation. The convergence of this approximation to a unique bounded variation solution of the model is proved. Numerical results showing the accuracy of this scheme are presented.