  The electric field used in most electrorheological (ER) experiments is
usually quite high, and nonlinear ER effects have been theoretically predicted
and experimentally measured recently. A direct method of measuring the
nonlinear ER effects is to examine the frequency dependence of the same
effects. For a sinusoidal applied field, we calculate the ac response which
generally includes higher harmonics. In is work, we develop a multiple image
formula, and calculate the total dipole moments of a pair of dielectric
spheres, embedded in a nonlinear host. The higher harmonics due to the
nonlinearity are calculated systematically.
