In this paper we present a systematic way of computing the polarization tensors,anisotropic as well as isotropic, based on the boundary integral method. We then use this method to compute the anisotropic polarization tensor for ellipses and ellipsoids. The computation reveals the pair of anisotropy and ellipses which produce the same polarization tensors.