We study a numerical method of calculating the near field of ensembles of arbitrary spheres by extending Mie theory. A recursive method based on the orders of scattering is presented. This method represents a concise way to calculate the near field of aggregates of any number of arbitrary spheres, Numerical examples are given to show its validity.