Microstructural evolution in three dimensions of nucleation and growth transformations is simulated by means of cellular automata (CA). In the simulation, nuclei are located in space according to a heterogeneous Poisson point processes. The simulation is compared with exact analytical solution recently obtained by Rios and Villa supposing that the intensity is a harmonic function of the spatial coordinate. The simulated data gives very good agreement with the analytical solution provided that the correct shape factor for the growing CA grains is used. This good agreement is auspicious because the analytical expressions were derived and thus are exact only if the shape of the growing regions is spherical.