Formulation involving the electric field integral equation (EFIE) and the related numerical solution via Method of Moments (MoM) variants requires tedious calculation of various double surface integrals arising from the use of vector triangular basis functions. The approach to evaluate such an integral strongly depends on the distance between source and observation points, respectively. In the case of a rather close interaction between the elements, due to a strong singularity, a purely numerical approach could fail resulting in the incorrect, spurious results. Hence, a combined analytical-numerical approach to extract the mentioned singularities is necessary. The present paper deals with an implementation of such an approach to the calculation of certain integrals types arising in the EFIE formulation for electromagnetic scattering from dielectric objects.