# 4. A Numerical Approximation of the Sphere to Sphere Force Calculation

This approximation again uses a simple spreadsheet, basic trigonometry and basic algebra to determine the correct electrostatic force between two spheres using the corrected equation derived above for the point to sphere force calculations.

For the purposes of this approximation the symmetry of a circle is used to approximate a sphere.  Unit charge is also assumed for simplicity.

The illustration below demonstrates the approach of calculating the force between a sphere and a sphere by approximating a circle with many points, then simply summing the force between the point outside the circle and each of the many points on the circle using the corrected equation F_{e}=q_{1}q_{2}d/(d^2-r^2)^{3/2}) , not Coulomb’s Law, for the point to sphere force calculation.

Using a Cartesian plane to calculate the forces, Point 1, the centroid of Circle 1 lies outside the circle, the centroid of the circle is 0,0.  Point 2 moves around the circle allowing many force calculations between Point 1 and Point 2 as Point 2 revolves around the circle using the corrected equation for the point to point force calculation (Point 1 is a sphere).  Point 1 remains stationary.

Assuming 360 points on the circle, one for each degree, each point will have 1/360th of the charge of the entire circle.  Forces in the Y direction are ignored as the corresponding point on the opposite side of the circle (across the X axis) exerts an equal and cancelling force in the Y direction.

The summed individual forces of the numerical approximation of a sphere to sphere do not equal the force calculated by Coulomb’s Law where two spheres are approximated as a points. Figure 6 – Approach used for the numerical approximation of sphere to a sphere forces using the corrected equation.

Note r_{1} is used in the force equation and  used in the calculation of distance for the force.

Using the cosine law  c=(r_{2}^2 + d^2 -2r_{2}d \cdot cos(\alpha))^{0.5}

The force along c is F_{e}=q_{1}q_{2}d/((c^{2}-r_{1}^{2})^{(3/2)})

Using basic trigonometry determine the X axis component of the force for each point as follows;

\beta =\pi-\alpha h=r_{2}\sin(\beta) \theta=a\sin(h/c)

Resultant force along X axis = =F_{e} \cos(\theta))

Steps to reproduce point to sphere numerical force approximation using Coulomb’s Law in a spreadsheet;

For simplicity assume unit force on the point and the circle.

360 degrees are converted to rads purely to explain the concept in degrees.

Create a spreadsheet with nine columns;

Summing column nine gives the total force between the two spheres using the corrected equation in a numerical approximation.

A sample spreadsheet table is below.  The full table is in Appendix 3.