The Coulomb potential between spherical-deformed and deformed-deformed nuclei has been calculated using the Monte Carlo simulation. The results obtained for the Coulomb potential in the 16O + 238U and 27Al+70Ge reactions are in good agreement with those obtained using the double-folding method. The simulation technique employed here has the ability of calculating the Coulomb potential taking into account the finite diffuseness parameter, all the possible deformation degrees of freedom, and different orientations of the symmetry axes of the target and the projectile nuclei with respect to each other. The accuracy of this simulation technique is high and the computer time taken to do these calculations is much less than those of the double-folding method.