We present a multilevel fast multipole solution method for volume integral equations dealing with scattering from arbitrarily shaped inhomogeneous dielectrics. The solution accuracy, convergence, computer time and memory savings of the method are de