This paper deals with the numerical simulation of the nucleation and crystal growth processes involved in the batch precipitation of photographic emulsions by the simultaneous addition of solutions of silver nitrate and sodium bromide to an aqueous solution of sodium bromide and gelatin. These processes have been successfully modelled in a previous paper (Muhr et at, 1995, Chem. Engng Sci. 50, 345-355). Two methods used to solve the population balance equation in order to follow the crystal size distribution are described. The partial differential equation which represents the dynamics of the particle size density function psi is transformed into a series of ordinary differential equations by means of the method of classes or by the direct discretization of psi. After a comparison of both methods, the second option was selected for the calculations. To this series of equations ale added balances for the total volume, moles of silver and moles of bromide. The aim of the present paper is to point out several details of the simulation. Results from example simulations are reported; the total number of crystals, the average crystal size, the standard deviation of the crystal size, and the crystal size density function are calculated. Special emphasis is set on the problems of grid spacing and representation of the primary nucleation process.