The Coulomb interaction for a sphericalâ€”deformed interacting pair is derived assuming realistic nuclear charge distributions. The effect of a finite diffuseness parameter is described either by the folding product of spherical or deformed sharp-surface distribution and a spherical short-range function or by using a Fermi two-parameter distribution function. The approximate solutions obtained using these categories of charge distributions are then compared to the numerical solution obtained within the framework of the double-folding model. We found that the finite surface diffuseness parameter affects slightly the inner region of the total Coulomb potential, while it produces large errors in calculating the Coulomb form factors used frequently in nuclear reactions and fusion numerical codes.