Ygtjki

I hate to be that guy to just post an integration problem and ask how to solve it so I'll give a little relevant info

Okay, so I'm working on a physics project and my professor proposed that the following double integral could potentially solve a problem that I've used an alternative method to solve:

$$I=int_0^piint_0^rhofract^2sinphileft(tcosphi-dright)left[t^2sin^2phi+left(tcosphi-dright)^2right]^3/2;dt dphi$$

1. 
   $rho$ is an arbitrary, strictly positive real constant


2. 
   $d$ is a real constant that satisfies $d>rho$
   

This integral's value could provide immense insight into fields of uniform, solid spherical objects, so it's actually pretty important for my work.

After some quick attempts to simplify, I decided to try some integral calculators with set values. Needless to say, the result after the first integral seemed so hopeless that I couldn't imagine simplifying and integrating again--not to mention then generalising constant inputs to their original variable form.

However, there is a strong likelihood that $I$ simplifies to one of the following two solutions:

$$text1.This solution comes from inverse square laws$$

$$I=frac1d^2$$

$$text2. This solution comes from a separate computation that I did (integrals below)$$

$$I=left(1-fracrho^25d^2right)left[frac32rho^2+frac3(rho^2-d^2)4drho^3lnleft(fracd+rhod-rhoright)right]$$

Although it looks like these are vastly different answers, given $rho=1$ and $d=10$, you get the following outputs from $(1)$ and $(2)$:

$$1.; I=0.01$$
$$2.; Iapprox 0.01000046$$

Here's the ratio of solution (2) over (1) for $rhoin(0,1),;din(0,50)$

I tried to tackle this problem differently than my professor, and set up the following integrals to solve the problem that lead to solution $(2)$:

$$frac94rho^6left[;intlimits_d-rho^d+rhoxleft[x-fracx^2+d^2-rho^22dright]left[frac(x+d)^2-rho^24dcdot xright];dxright]cdotleft[;intlimits_d-rho^d+rhofracrho^2-(x-d)^22dcdot x;dxright]$$

Where you come in

If the double integral is correctly composed (which my professor felt confident with), I need someone skilled in integration to solve said double integral. I've given two possible solutions and it's probable that the answer will be one of those. If it's solution $(1)$, I know that mine will have an error and you will essentially have proved the inverse square law for gravitational and electric fields. If it's solution $(2)$, then this will be far more exciting to me but less likely. If it's neither, then there are several possible implications

BOUNTY

I'm willing to award the following bounties for solving the double integral at the beginning. Since certain solutions have stronger implications (as explained above), I'm rewarding the following bounties:

1. +200 rep if you verify solution $(1)$
   


2. +500 rep if you verify solution $(2)$
   


3. +75 rep for any other solutions (note they'll have to be verified by a second user first)

QUESTIONS

If you have any additional questions feel free to ask, and thanks for reading all this!

As a mathematician, I would divide by force in the numerator and denominator by $d^3$, substitute $t/d$ by something, thus reducing to the case $d=1$. But here, let it be, we conserve the homogeneous setting as a control of the computations.

We split the numerator, compute first
$$
beginaligned
J_1
&=
int_0^rho
dt
int_0^pi
frac
t^2sinphicdot tcosphi
left[t^2sin^2phi+left(tcosphi-Dright)^2right]^3/2; dphi
\
&=
int_0^rho
dt
int_0^pi
frac
t^2(-cosphi)'cdot tcosphi
left[t^2-2Dtcosphi+D^2right]^3/2; dphi
\
&qquadtext Substitution: u=cos phi ,
\
&=
int_0^rho
dt
int_-1^1
frac
t^3; u
left[t^2-2Dt;u+D^2right]^3/2; du
\
&qquadtext Substitution (for $u$, fixed $t$) of the radical v=sqrtt^2-2Dt;u+D^2 ,
\
&qquad u=frac 12Dt(t^2+D^2-v^2) , du=-frac vDt; dv\ ,
\
&=
-
int_0^rho
dt
int_sqrtt^2+2Dt+D^2^sqrtt^2-2Dt+D^2
frac
t^3; frac 12Dt(t^2+D^2-v^2)
v^3; frac vDt; dv
\
&=
int_0^rho
t;dt
int_D-t^D+t
frac 12D^2
cdot
frac t^2+D^2-v^2
v^2; dv
\
&=
int_0^rho
t;dt
;frac 12D^2
left[
-(t^2+D^2)frac 1v
-1
right]_v=D-t^v=D+t
\
&=
int_0^rho
dt
;frac t2D^2
left[
(t^2+D^2)left(frac 1D-t-frac 1D+tright)
-
2t
right]
\
&=
int_0^rho
dt
left[
frac DD+t
+frac DD-t
-2fracD^2+t^2D^2
right]
\
&=
Dlnfrac D+tD-t
-
2rholeft(1+frac rho^23D^2right)
.
endaligned
$$
Computer check for $D=2$, $rho=1$ (pari/gp code):

? D=2; r=1;
? intnum(t=0,r, intnum(s=0, Pi, t^2*sin(s)*t*cos(s) / (t^2-2*t*D*cos(s)+D^2)^(3/2) ) )
%19 = 0.030557910669552716123823807178384744388
? D*log( (D+r)/(D-r) ) - 2*r*(1+r^2/3/D^2)
%20 = 0.030557910669552716123823807178384742634
? 
? D=223; r=101;
? intnum(t=0,r, intnum(s=0, Pi, t^2*sin(s)*t*cos(s) / (t^2-2*t*D*cos(s)+D^2)^(3/2) ) )
%22 = 1.9969022076015148346071622544965636670
? D*log( (D+r)/(D-r) ) - 2*r*(1+r^2/3/D^2)
%23 = 1.9969022076015148346071622544965636629

The other integral. I will integrate here first w.r.t. $t$.

$$
beginaligned
J_2
&=
-D
int_0^pi
dphi
int_0^rho
frac
t^2
left[(t-Dcosphi)^2+Dsin^2phiright]^3/2
; dt
\
&qquadtext and we consider separately (without the factor $-D$)
\
J_2(phi)
&=
int_0^rho
frac
t^2
left[(t-Dcosphi)^2+Dsin^2phiright]^3/2
; dt
\
&=
int_0-Dcosphi^rho-Dcosphi
frac
(u+Dcosphi)^2
(u^2+a^2)^3/2
; du ,qquad a:= Dsinphi
.
\
&qquad
text Now the integrals can be computed
\
int fracu^2
(u^2+a^2)^3/2
; dt
&=
-frac t(u^2+a^2)^1/2+operatornamearcsinh frac ta+C ,
\
int fracu
(u^2+a^2)^3/2
; dt
&=
-frac 1(u^2+a^2)^1/2+C ,
\
int frac1
(u^2+a^2)^3/2
; dt
&=
-frac a^2;u(u^2+a^2)^1/2+C ,
endaligned
$$
and the computation goes on.
If my calculus is ok, then
$$
beginaligned
J_2(phi)
&=
int_0^pi
dphi;
Bigg[
operatornamearcsinh fract-Dcos phiDsinphi
\&qquadqquadqquad+
fract-Dcosphi(t^2-2Dtcosphi+D^2)^1/2sin^2phi
\&qquadqquadqquadqquadqquadqquad
+frac2(t^2-2Dtcosphi+D^2)^1/2
Bigg]_0^rho .
endaligned
$$
I have to submit, hope this is helpful to check with the own computations.
I'll be back, but typing kills a lot of time.

Hint:

With the change of variable $u=cosphi$, the integral on $phi$ becomes

$$int_-1^1fract^2(tu-d)sqrt(u-dt)^2+d^2(1-t^2)du.$$

By decomposition of the numerator, you will get a term

$$c(t)log((u-dt)^2+d^2(1-t^2))$$

and another

$$c'(t)arctanfracu-dtdsqrt1-t^2.$$

These terms do not simplify at the bounds of the integration interval.

The integral on $t$ (cubic in $t$ at the denominator) is worse. I am not optimisitc about existence of a closed-form.

beginalign*
&iint fract^2 sin(phi) (t cos(phi) - d)(t^2 sin^2(phi) + (t cos(phi) - d)^2)^3/2 ,mathrmdt,mathrmdphi \
&= fracsqrtd^2 + t^2 - 2 d t cos(phi)(d^2 - 2 t^2 - 2 d t cos(phi) - 3 d^2 cos(2 phi))6d^2 \
&+ d cos(phi) lnleft(t - d cos(phi) + sqrtd^2 + t^2 - 2 d t cos(phi)right) sin^2(phi) text, endalign*
as one can readily verify. Then $I = frac-2 rho^33 d^2$.

I think for your case $1$, you mean $I propto frac1d^2$. The integral can't be positive because:

• 
  $t^2 geq 0$ and


• 
  $sin(phi) geq 0$ since $phi in [0,pi]$, but


• 
  $t cos(phi) - d < 0$ because $0 < t < rho < d$, while


• the denominator is $geq 0$, so


• the integrand is (zero or) negative everywhere.