We consider a two- or three-dimensional body of arbitrary shape containing foreign inclusions. These inclusions are in perfect electromagnetic contact with the matrix. A boundary condition of the third kind is imposed on the surface of the body. We propose a method of combined application of the fundamental solution of the nonstationary of heat-conduction equation, near-boundary and contact elements, and a step-by-step time scheme of common initial condition for the construction of integral transformations of the components of the vectors of electromagnetic-field intensity at an arbitrary point of space and time. In the case of perfect electromagnetic contact between the zones, the application of contact elements on the interface of the media (instead of boundary or near-boundary elements) enables one to automatically satisfy the first condition of contact (the equality for the components of the vector of electric-field intensity). As a result, the number of unknown fictitious current sources decreases as compared with the case of application of traditional indirect methods of boundary and near-boundary elements.