**Applied Mathematics**

Vol.08 No.06(2017), Article ID:77059,15 pages

10.4236/am.2017.86064

The Gravitational Attraction between Hemispherical Masses

Jonathan D. Weiss^{ }

JSA Photonics LLC, 8019 Edith Blvd NE, Albuquerque, NM, USA

Copyright © 2017 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: May 11, 2017; Accepted: June 18, 2017; Published: June 21, 2017

ABSTRACT

This paper is a study of the gravitational attraction between two uniform hemispherical masses aligned such that the pair is cylindrically symmetric. Three variations are considered: flat side to flat side, curved side to curved side, and flat side to curved side. Expressions for the second and third variation are derived from the first, with the use of superposition and the well- known gravitational behavior of a spherical mass as equivalent to a point mass at its center. The study covers two masses of equal diameter and of different diameters, such that one is four times that of the other. Calculations are done for separations from zero to fifty times the radius of the larger of the two, which is effectively the asymptotic limit. It is demonstrated that at any separation, the force can be expressed as if the two hemispheres were point masses separated by a certain distance. Expressions for that distance and the location of the (fictitious) point masses within each hemisphere are presented. Unlike the case of two spherical masses, the location within their respective hemisphere is not necessarily the same for each point and both are dependent upon the separation between the two hemispheres. The calculation for the first variation is done in two ways. The first is a “brute force” multi-dimensional integral with the help of Wolfram Mathematica. The second is an axial expansion for the potential modified for off-axis locations by Legendre polynomials. With only a few terms in the expansion, the results of the second method are in extremely good agreement with those of the first. Finally, an interesting application to a split earth is presented.

**Keywords:**

Gravitational Attraction, Hemispherical Solids, Legendre Polynomials, Split Earth

1. Introduction

With at least one exception, a perusal of online publications of the (Newtonian!) gravitational attraction between solid masses yields discussions of a pedagogical nature concerning spherically symmetric solids and shells. In a few cases, the force between a point mass and a disk or a hemisphere is derived for the point mass placed along the center line of the non-spherical object. The exception is a calculation of the force between two hemispheres with their flat surfaces in contact, as in the case of the earth [1] . In addition to considering pairs of non- spherically symmetric masses, this paper generalizes the discussion in [1] to the force between two hemispheres of varying separation, such that the pair forms a cylindrically-symmetric configuration. Furthermore, two additional relative orientations are analyzed, along with the one involving both flat surfaces facing each other, as in [1] . These involve the curved surfaces facing each other, or a curved surface facing a flat surface. As a further generalization, two sets of sizes are also considered. The first set involves both hemispheres of the same diameter, while the second involves diameters that differ from each other by a factor of four. The situation is shown in Figure 1, where 2A and 2B merge when B = A. The symmetry dictates that the net force between the two hemispheres is in the vertical direction. Calculations of the force will be done for separations between the closest two points, S, such that 0 ≤ S ≤ 50 A, essentially the asymptotic limit.

Figure 2 illustrates how the force associated with the second two arrangements can be derived from the force associated with the flat-to-flat arrangement through linear superposition of forces. Thus only the first arrangement requires a nontrivial calculation. The arrangements involving a sphere and a hemisphere are equivalent to a point source outside the hemisphere and positioned along its center line. This calculation is comparatively easy to perform.

In addition to the force calculations, it will be shown that the force can be expressed in the form of Newton’s law of gravitation for point masses, GM_{A}M_{B}/ (R_{eff})^{2}, where R_{eff} is an effective separation between the two point masses, each having the same mass as the associated hemisphere. Only the magnitude of the force is considered here. This situation is shown in Figure 3 for the flat-to-flat configuration. The location of each of the two fictitious point masses is within its respective hemisphere at a position α_{A}A or α_{B}B from the flat surface. The two α’s

Figure 1. The four configurations analyzed in this paper, showing the axis of symmetry. The mass and radius of the lower (upper) hemisphere are M_{A} (M_{B}) and A (B), respectively.

Figure 2. Illustrations of how configurations 2A, 2B, and 3 can be determined from the first configuration, thereby simplifying the analysis.

Figure 3. Illustration of how the force can be expressed in terms of an effective separation between two point masses and how to determine the location of those point masses within the hemispheres.

refer to the fraction of the radius, A or B, where the point mass is located. If a curved surface faces down, the position will be measured from that curved surface. In all cases, the position will be calculated as a function of S and shown to approach the center of mass of the respective hemisphere at infinite separation. The separation of 50 A is fully adequate to demonstrate that point. This center-of-mass limit ((3/8) A (or B), measured from the flat surface, and (5/8) A (or B), measured from the curved surface) can be proved analytically and can be viewed as a test of the accuracy of the calculation. Their position at any separation will be calculated, as shown in the bottom of Figure 3, by replacing one of the hemispheres by its equivalent point mass at the unknown position α_{B}B and calculating the force between it and the other hemisphere. That force between the two will be equated to the force previously calculated between the two hemispheres, thus determining α_{B}B. Since
${\alpha}_{A}A+{\alpha}_{B}B+S={R}_{eff}$
, the other location is determined.

2. Calculations

2.1. Preliminary

First, let us be specific about the superposition of collinear forces, at least in one case. Others are similar. Assuming the form GM_{A}M_{B}/(R_{eff})^{2} and referring to Figure 2, panel 2A, we find that

$\frac{1}{{\left[{R}_{eff}\left(cf;{S}_{2}\right)\right]}^{2}}=\frac{2}{{\left[{R}_{ef}f\left(pf;{S}_{2}\right)\right]}^{2}}-\frac{1}{{\left[{R}_{eff}\left(ff;{S}_{2}+B\right)\right]}^{2}}$ (1)

where, in order of their appearance, “cf” refers to curved-flat; “pf” refers to point-flat (because the sphere acts as a point source), and “ff” refers to flat-flat, as also shown in the first panel. The separation for each R_{eff}, S_{2} and S_{2} + B, are also shown. The factor of 2 exists because the mass of the sphere is twice that of the hemisphere.

Next, we present the formulas for the magnitude of the force between a sphere and a hemisphere required in Figure 2(A), Figure 2(B), and 3 of Figure 2. In panels 2A and 2B, the flat surface of the hemisphere faces the sphere, but with A and B reversed. In panel 3, the curved side of the hemisphere faces the sphere. A simple integration yields the following for the force, F_{pf}, in Figure 2(A):

where we identify the multiplier of [GM_{A}(2M_{B})] with (1/R_{eff})^{2}, S = S_{2} + B, and M_{A} and M_{B} are the two hemispherical masses. In the case of panel 2B, A and B are reversed in the above formula. The equation for the force between the curved side facing the sphere and the sphere can be obtained by a similar integration, resulting in a similar formula. More simply, it can be obtained by subtracting from the force between two spheres the force between the sphere and the hemisphere (flat side facing the sphere), as shown in Figure 4. The force between the two spheres is simply proportional to
$4/{\left({S}_{3}+A+B\right)}^{2}$
.

2.2. Force between Two Hemispheres: Flat Surface Facing Flat Surface

As was stated earlier, this situation is the only one requiring a nontrivial calculation. It was performed in two ways. The first involves a direct integration of the force between all of the points in one hemisphere and all of the points in the

Figure 4. An illustration of the use of superposition to determine the force between a sphere and a hemisphere, curved side facing forward, needed to analyze configuration 3 in Figure 1.

other. Because integration over two volumes was required, it started out as a six-fold integration. However, azimuthal symmetry reduced it by one integral and another integral could be done analytically, thus reducing the problem to a still-unpleasant four-fold integration. The second way involved expansions of the potential created by one hemisphere along its axis for positions z ≤ A and z ≥ A. The off-axis values were then expressed in terms of the original expansion multiplied by Legendre polynomials, as in electrostatics.

3. Direct Integration Method

The direct integration method was performed with respect to Figure 5, where two hemispheres, along with the two sets of dimensions and the separation between them are shown.

The basic integral expression for the point-to-point (pp) force, F_{pp}, is:

${F}_{pp}=G{\rho}_{A}{\rho}_{B}{\displaystyle \iint \frac{\text{d}{V}_{A}\text{d}{V}_{B}\mathrm{cos}\phi}{{R}^{2}}}$ (3)

where V_{A} and V_{B} are the volume elements in their respective hemisphere; R is the separation between any two points in different hemispheres; φ is the angle between R and the vertical; and ρ_{A}_{(B)} is the mass density of hemisphere A(B), assumed to be different for the sake of generality. By symmetry, the force is in the vertical direction. This equation can be written in an explicit form:

Writing the equation in the form of the product of the two point masses divided by the square of an effective separation between them scaled by the radius A(1/R_{eff}_{,A})^{2}, we obtain the following complicated expression:

${\left(\frac{1}{{R}_{eff,A}}\right)}^{2}=\frac{9}{\text{\pi}}{\displaystyle {\int}_{0}^{1}\text{d}{z}_{b}}{\displaystyle {\int}_{0}^{{\left(1-{z}_{b}^{2}\right)}^{0.5}}{R}_{b}\text{d}{R}_{b}}{\displaystyle {\int}_{0}^{\text{\pi}}\text{d}{\phi}_{a}}{\displaystyle {\int}_{0}^{1}\text{d}{z}_{a}\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}\left[\right]$ (5)

Figure 5. The geometry needed for the direct integration method. The angles φ_{A} and φ_{B} are azimuthal, defining the direction between the horizontal vectors, r_{A} and r_{B}, and the dotted reference lines.

where

$\begin{array}{c}\left[\right]=\frac{{\left[{\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}^{2}+{\left(b{R}_{b}\right)}^{2}\right]}^{1/2}}{{\left(bR{r}_{b}\mathrm{sin}{\phi}_{a}\right)}^{2}+{\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}^{2}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{\left(b{R}_{b}\mathrm{cos}{\phi}_{a}\right){\left(1-{z}_{a}^{2}\right)}^{1/2}-\left({\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}^{2}+{\left(b{R}_{b}\right)}^{2}\right)}{\left({\left(b{R}_{b}\mathrm{sin}{\phi}_{a}\right)}^{2}+{\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}^{2}\right){X}^{\frac{1}{2}}}\end{array}$

And where

$X=\left(1-{z}^{2}\right)-\left(2b{R}_{b}\mathrm{cos}{\phi}_{a}\right){\left(1-{z}_{a}^{2}\right)}^{1/2}+\left({\left({z}_{a}+b{z}_{b}+{s}_{a}\right)}^{2}+{\left(b{R}_{b}\right)}^{2}\right)$

R_{eff}_{,A} = R_{eff}/A. All the other linear parameters in the expression are similarly scaled and use similar notation (e.g., s_{a} ≡ S/A; b = B/A; z_{a} = Z_{A}/A; r_{a} = R_{A}/A; z_{b} = Z_{B}/B; r_{b} = R_{B}/B). The integral over the azimuthal angle, φ_{B}, forthe upper hemisphere was simply replaced by 2π, and the integral over r_{A} was doable analytically. Thus, as stated above, the 6-fold integral was reduced to a 4-fold integral. Since the 4-fold integral did not appear to be analytically evaluable, it was numerically evaluated by using the Wolfram Mathematica command “NIntegrate” as applied to multiple integration. The results of the direct integration are shown in Figure 6 for b = 1 and b = 1/4. Configuration 1 refers to the first panel in Figure 1. The basic behavior is the same for both values of b: a rapid drop off for small values of S_{A}, followed by a transition to very small values of (1/R_{eff}_{,A})^{2} at large values of

Figure 6. The inverse square of the scaled separation, (1/R_{eff}_{,A})^{2}, for the 1:1 and the 0.25:1 of configuration 1 over separations from 0 to 50 A. The values for b = 1/4 are larger than for b = 1.

S_{A}. The values for b = 1/4 are consistently larger than those for b = 1. At S = 0, for example, the two are 1.31 and 0.75, respectively. The force for b = 1 is considerably larger, however, because the product of the masses is 64 times as large. Given the value of 0.75 for (1/R_{eff}_{,A})^{2} with b = 1, an expression for the magnitude of the force between the two hemispheres in contact with one another is (1/3) [G(πρA^{2})^{2}], where ρ is the mass density of both hemispheres. This is the same result obtained in [1] by different methods from the one used here (see Equations (5) and 19 of that reference).

3.1. Gravitational Force between Two Halves of a Sphere

Since the gravitational force between the two halves of a solid sphere, such as the earth, is of some interest to inquiring minds [2] , we shall consider two applications of the result immediately above. The first concerns how the gravitational force between two contacting spheres, each of the same mass as the individual hemispheres, compares to the force between the two hemispheres. Since the ratio of the radii, sphere to hemisphere, is 1/2^{1/3}, the ratio of the two forces, spheres to hemispheres, can easily be deduced as (2^{2/3}/3) ≈ 0.53. Thus, the force between the two hemispheres is almost twice that of the two spheres. The second application is concerned with the rotation rate of the earth required to separate its two halves under only the combined action of the fictitious centrifugal force and gravity. The calculation is done with reference to Figure 7, in which is shown a sphere of radius A split along a plane perpendicular to the plane of the paper and passing through the central axis of rotation. The split sphere is rotating with angular velocity ω. The coordinates (r, φ, z) define the position of a rotating point. The total “centrifugal force”, F_{c}, is obtained by integrating the component of the force perpendicular to the interface plane over the hemisphere.

Figure 7. The geometry of gravitational force between the two halves of a solid sphere, such as the earth.

${F}_{c}=2{\omega}^{2}\rho {\displaystyle {\int}_{0}^{A}\text{d}z}{\displaystyle {\int}_{0}^{{\left({A}^{2}-{z}^{2}\right)}^{1/2}}{r}^{2}\text{d}r}{\displaystyle {\int}_{-\frac{\text{\pi}}{2}}^{\frac{\text{\pi}}{2}}\mathrm{cos}\phi \text{d}\phi}=\frac{\text{\pi}}{4}{\omega}^{2}\rho {A}^{4}$ (6)

Equating F_{c} to [G(πρA^{2})^{2}]/3, solving for ω, and expressing ω as 2πN_{r}, where N_{r} is the number of rotations per second, we obtain N_{r} = [Gρ/(3π)]^{1/2}. This results in ≈ 20 × 10^{−5} rotations/sec. required to just separate the two halves of the sphere. The earth’s rotational speed ≈ 1.16 × 10^{−5} rotations/sec., which ≈ (1/17) of the required speed. Of course, the separation speed would be much greater if the cohesion within the earth, which serves a kind of geological corpus callosum, were taken into account. Furthermore, tidal friction is causing the earth’s rotational speed to slowly diminish. Thus, we have nothing to worry about.

3.2. An Array of Results

In this section, we present a graphical array of results, in addition to the one already presented in Figure 6. The graphical array will be of (1/R_{eff}_{,A})^{2} vs. S_{A} for the various configurations in Figure 1 and of Alpha (the spelled out version of the Greek letters in Figure 2) for those same configurations. Figure 8 is a more detailed view of the behavior of (1/R_{eff}_{,A})^{2} for S_{A} ≤ 3, b = 1, and for its three distinguishable configurations. The dotted line is added to show the linear behavior in the immediate vicinity of S_{A} = 0. The ordering in size is consistent with the amount of mass in one hemisphere that is close to the other hemisphere. That order is flat-to-flat, flat-to-curved, curved-to-curved. Figure 9 is similar to Figure 8, but with b = 0.25. For these graphs, (1/R_{eff}_{,A})^{2} is larger than in Figure 8 because one of the hemispheres is so much smaller. The ordering in size in Figure 9, in general, is also consistent with the amount of mass in one hemisphere that is close to the other. Curiously, however, there is a cross over between con-

Figure 8. A view of (1/R_{eff,A})^{2} for the three distinguishable configurations for b = 1 over the range over separations from 0 to 3.

Figure 9. Similar to Figure 8, but for b = 0.25. Configurations 2A and 2B are now distinguishable.

figurations 2A and 2B at S_{A} = 0.25, which is maintained down to S_{A} = 0. Figure 10 and Figure 11 are the same as Figure 8 and Figure 9, except that the entire range of values are shown on a log-log scale. The dashed line in Figure 10 shows how the behavior at large values of S_{A} varies as
$1/{S}_{A}^{2}$
. Figure 12 and Figure 13 illustrate the behavior of Alpha for various configurations and the two values of b (1 and 1/4). Alpha 1 (or Alpha 1/4) is the value of Alpha for a hemispherical radius = 1 (or 1/4). As expected from the comment in the introduction, this parameter approaches (3/8) or (5/8) for large values of S_{A}, depending on whether

Figure 10. The same as Figure 8, but on a log-log scale for separations between 0 and 50 A. The dashed line behaves as [(1/(S/A)]^{2}.

Figure 11. The same as Figure 9, but on a log-log scale for separations between 0 to 50 A. Configurations 2B and 3 are indistinguishable on this log-log graph.

Figure 12. The top two curves, configuration 2A curved and configuration 3 curved, are indistinguishable on this graph. However, Alpha 1 for configuration 3 is slightly greater than that for configuration 2A because the force between the two hemispheres is slightly greater in configuration 2A, as is made clear in Figure 1. A similar comment applies to the other pair of curves, where Alpha 1 for configuration 2 flat is slightly greater than that for configuration 1 flat.

Figure 13. Alpha 1/4 vs. S/A from 0 to 50 A for b = 1/4. The top two curves, which are almost indistinguishable, represent the curved side facing forward, and approach 5/8 at large S/A. The bottom curve represents the flat side facing forward and approaches 3/8 at large S/A.

the flat or curved side faces the other hemisphere. That condition is denoted by “flat” or “curved” in Figure 12. The curves in Figure 12 and Figure 13 can be interpreted as indicating that the force between the two hemispheres decays away more slowly than 1/S^{2} and only approaches that behavior at very large vales of S (as shown in Figure 10 and Figure 11). As Alpha diminishes, the increasing separation between the two fictitious point masses with S is partially offset by their move toward each other within their respective hemisphere.

3.3. The Use of Power Expansions and Legendre Polynomials

Without the use of high-speed computers and sophisticated computational software, such as Wolfram Mathematica or, perhaps, MATLAB, it would be extremely difficult to perform the four-fold integration that makes up Equation (5). This is particularly true because of the coupling among the variables of integration. However, the use of a power expansion of the on-axis gravitational potential of hemisphere A (or B), followed by a term-by-term multiplication by Legendre polynomials to produce off-axis values of the potential simplifies the calculation. This is commonly done when dealing with electrostatics problems containing azimuthal symmetry, and it is discussed in many textbooks [3] . This method can be applied to the current problem because, as in electrostatics, it requires a solution of Laplace’s equation. The only possible drawback is that the result is an infinite series, which may require a largenumber of terms to produce the desired accuracy. Fortunately, this will not be true of the number of terms required for an accurate solution to the current problem.

This method of solution for configuration 1 is presented in conjunction with Figure 14, where a hemisphere of radius A, assumed to be the source of the potential, Φ(R, θ), is shown with the hemisphere of radius B, assumed to be in the gravitational field of hemisphere A.

Once the potential, Φ(r, θ), is determined, its negative z-derivative (the z- component of its gravitational field) is calculated and integrated over hemisphere B to determine the total force between the two hemispheres. In the case of Figure 14, when R ≤ A, the on-axis potential can be written as a power series in R^{l}, whereas for R ≥ A, it can be written as a power series in R^{−(}^{l}^{+1)}. Thus, the potential at angle θ is of the form:

Figure 14. The geometry for the use of power series expansions and Legendre Polynomials for configuration 1.

$\Phi \left(R,\theta \right)={\displaystyle \sum}_{l=0}^{\infty}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{C}_{l}{R}^{l}{P}_{l}\left(\mathrm{cos}\theta \right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}R\le A$ (7a)

$\Phi \left(R,\theta \right)={\displaystyle \sum}_{l=0}^{\infty}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{D}_{l}{R}^{-\left(l+1\right)}{P}_{l}\left(\mathrm{cos}\theta \right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}R\ge A$ (7b)

C_{l} and D_{l} are parameters to be determined and
${P}_{l}\left(\mathrm{cos}\theta \right)$
is the Legendre Polynomial of order l, all of which are unity at θ = 0. In order to proceed, we must first determine the on-axis potential. After a simple integration it can be shown to be:

$\Phi \left(R,0\right)=-G\left(2\text{\pi}{\rho}_{a}\right)\left[\frac{{\left(A+R\right)}^{3}}{3R}-\frac{{\left({A}^{2}+{R}^{2}\right)}^{\frac{3}{2}}}{3R}-AR-\frac{{A}^{2}}{2}\right]$ (8)

Ignoring constants in the above expression, which contribute nothing to the force, results in the following expansions for R ≤ A and R ≥ A, respectively.

and

From these two equations, the constants C_{l} and D_{l} can be directly determined. To obtain Φ(R, θ), we multiply each term in the above two expansions by the appropriate Legendre polynomial, express the expansion in terms of the hemispheric coordinates z, r, φ and, to obtain the force, integrate its negative z-deriv- ative over z from s to
$\left(s+{\left({B}^{2}\u2013{r}^{2}\right)}^{1/2}\right)$
and over r from 0 to B. The azimuthal integration over φ is simply replaced by 2π. In addition, we note that
$\mathrm{cos}\theta =z/{\left({z}^{2}+{r}^{2}\right)}^{1/2}$
. When all of this and a few other manipulations are completed, we identify (1/R_{eff,A})^{2} for R ≤ A as:

We note the use of scaled coordinates and the fact that the z_{a}-integral is particularly simple because it is the integral of a derivative. The integral over r_{a} was done numerically. For R ≥ A, the results are similar:

For simplicity, only the first two terms of each expansion are shown in Equations (10a) and (10b). Figure 15 and Figure 16 show the results of this method for b = 1 and b = 0.25, respectively, along with those using Wolfram Mathematica. With only a very few terms of the expansion, the two sets are virtually indistinguishable on the graphs.

Figure 15. The results of the two methods for b = 1, showing how well the two methods agree after the use of only a very few terms in the expansion. The results from Wolfram Mathematica, the use of terms 1 - 2 of the series, and the use of terms 1 - 3 of the series are essentially indistinguishable.

Figure 16. Similar to Figure 15, but for b = 1/4. The results from Wolfram Mathematica, the use of terms 1 - 4 of the series, and terms 1 - 6 of the series are essentially indistinguishable.

4. Summary

This paper consists of a study of the gravitational attraction between two inline hemispheres over a range of separations from zero to 50 times the radius of the larger of the two hemispheres. Two relative radii were chosen, 1:1 and 0.25:1 as sufficiently representative, along with three relative orientations: flat side facing flat side, curved side facing curved side, and flat side facing curved side. Initially, the first was calculated through a “brute force”, four-fold numerical integration with the help of Wolfram Mathematica. The second and third orientations were derived from the first by a superposition of the first with a simpler, easily calculated, intermediate configuration. The results were expressed in terms of the canonical form of two point masses separated by an effective distance. That “effective distance” was calculated, as was the location of the “point masses” within the hemispheres. An application of the 1:1, flat-to-flat configuration was made to a split, rotating earth. Finally, the flat-to-flat configuration was analyzed using a power series expansion for the on-axis potential and Legendre polynomials. With the use of very few terms in the expansion, the results agreed extremely well with the four-fold integration.

Based on the methods of calculation presented here, future work could involve other pairs of masses that form a cylindrically symmetric configuration. Elongating or compressing one or both hemispheres along their individual axes immediately comes to mind.

Acknowledgements

The author is grateful to Joseph S. Accetta, CEO of JSA Photonics, and William F. Filter, retired Sandia National Laboratories physicist, for reviewing the manuscript and making valuable suggestions for its improvement.

Cite this paper

Weiss, J.D. (2017) The Gravitational Attraction between Hemispherical Masses. Applied Mathematics, 8, 820-834. https://doi.org/10.4236/am.2017.86064

References

- 1. Gray, A. (1906) Solutions of Physical Problems. Proceedings of the Royal Philosophical Society of Glasgow, 37, 192.
- 2. If the earth is cut into two halves, what will be the gravitational force between them? https://www.quora.com/If-the-earth-is-cut-into-two-halves-what-will-be-the-gravitational-force-between-them
- 3. Jackson, J.D. (1999) Classical Electrodynamics. 3rd Edition, John Wiley & Sons, New York.