The Double Folding (DF) model calculation of the internuclear potential in heavy-ion interactions when the participant nuclei are deformed in their ground states involves a six-dimensional integral. Using the multipole expansion in these calculations, the DF six-dimensional integral reduce to the sum of the products of three single-dimensional integrals. In this paper we have presented a procedure for the calculation of the radius dependent functions in the multipole expansion of the nuclear density and their Fourier transforms. We have also reduced the DF model integrals to the sum of the single dimensional integrals using the obtained relations for the radius dependent functions in the multipole expansion and their Fourier transforms.
