In this paper, the computation of multiple (including two dimensional and three dimensional) Cauchy principal integral with generalized composite rectangle rule was discussed, with the density function approximated by the middle rectangle rule, while the singular kernel was analytically calculated. Based on the expansion of the density function, the asymptotic expansion formulae of error functional are obtained. A series is constructed to approach the singular point, then the extrapolation algorithm is presented, and the convergence rate is proved. At last, some numerical examples are presented to validate the theoretical analysis.