We present a new integral equation formulation of the polarizable continuum model (PCM) which allows one to treat in a single approach dielectrics of different nature: standard isotropic liquids, intrinsically anisotropic medialike liquid crystals and solid matrices, or ionic solutions. The present work shows that integral equation methods may be used with success also for the latter cases, which are usually studied with three-dimensional methods, by far less competitive in terms of computational effort. We present the theoretical bases which underlie the method and some numerical tests which show both a complete equivalence with standard PCM versions for isotropic solvents, and a good efficiency for calculations with anisotropic dielectrics. (C) 1997 American Institute of Physics.