A hybrid numerical method is proposed involving a successful combination of the equivalent electrodes and boundary element method. Electric field and electric potential distribution are determined by using both equivalent electrodes and equivalent charges. Green's function for potential of the charge per unit length, placed in free space, on the boundary of two dielectrics, is used. Magnetic-field distribution is determined by applying equivalent currents through magnetic vector potential and substitution of Ampere's microscopic currents with total currents. The accuracy of the developed method is highlighted with some simple numerical examples.