A numerical solution to electromagnetic scattering by three-dimensional nonlinear objects in free space using a statistical cooling procedure