Originalmente fue preguntado aquí. Este también fue preguntado aquí.
Me he encontrado con algunas dificultades a la hora de la siguiente integral:
I=∫2π0dϕ∫π0dθ sinθ∫∞0dr r23x2y2cos(ursinθcosϕ)cos2θ(y2cosϕ+x2sin2ϕ)sin2θ+x2y2cos2θe−r22,
donde x, y, y u son reales constantes positivas. Traté de al menos dos maneras de resolver esta integral:
- Primer intento:
Comencé a resolver el r integrales primeras. Mediante el uso de Mathematica, entonces
I=∫2π0dϕ∫π0dθ sinθ3x2y2(1−u2sin2θcos2ϕ)cos2θ(y2cosϕ+x2sin2ϕ)sin2θ+x2y2cos2θe−u22sin2θcos2ϕ.
After that, I looked for a solution for ϕ integral. My best attempt was:
Iϕ(x,y,u,θ)=2B[B(12)F1(12,1,−;1;ν,−a2)−aB(12)F1(32,1,−;2;ν,−a2)],
where B=x2sin2θ+x2y2cos2θ, a=u2sin2θ, and ν=x2−y2x2+x2y2cot2θ. In this way, the final results it's something like that:
I=∫π0dθ 3x2y2sinθcos2θ Iϕ(x,y,u,θ)..
Eq. (3) no puede ser más simplicada, en general, y es el resultado nal.
- Segundo intento:
Para evitar que la función hipergeométrica F1, he intentado iniciar con el ϕ integral. En este caso, mi problema inicial es una parte integral algo así como que:
∫2π0dϕcos(Acosϕ)a2cos2ϕ+b2sin2ϕ.
This integral (4) can be solved by series (see Vincent's answer and Jack's answer). However those solutions, at least for me, has not a closed form. This is my final step on this second attempt :(
What is the point? It turns out that someone has managed to solve the integral (1), at least the integral in r and ϕ. El final resuls encontrado por esta persona fue:
IG=12πx y(1−x2)3/2∫√1−x20dkk2exp(−u22x2k2(1−x2)(1−k2))√1−k2√1−k21−y21−x2,
donde, creo, k=√1−x2cosθ. Como puedes ver en el siguiente código realiza en Mathematica
IG[x_, y_, u_] :=
Sqrt[Pi/2] NIntegrate[(12 Pi x y)/(1 - x^2)^(3/2)
(v^2 Exp[-(u^2 x^2 v^2)/(2 (1 - x^2) (1 - v^2))])/(Sqrt[1 - v^2] Sqrt[1 - v^2 (1 - y^2)/(1 - x^2)]), {v, 0, Sqrt[1 - x^2]}]
IG[.3, .4, 1]
** 4.53251 **
I[x_, y_, u_] :=
NIntegrate[(r^2 Sin[a] Cos[
u r Sin[a] Cos[b]] 3 x^2 y^2 Cos[a]^2 Exp[-r^2/
2])/((y^2 Cos[b]^2 + x^2 Sin[b]^2) Sin[a]^2 +
x^2 y^2 Cos[a]^2), {r, 0, Infinity}, {a, 0, Pi}, {b, 0, 2 Pi}]
I[.3, .4, 1]
** 4.53251 **
las integrales de I IG son iguales. De hecho, desde la que emergen desde el mismo problema físico.
Entonces, mi pregunta es: ¿cuáles son los pasos que se aplica para que la integral de la I da la integral de la IG?
Editar
Desde que mi pregunta no fue resuelto todavía, yo creo que es porque es una pregunta difícil, voy a mostrar un caso particular de la integral de la I, dejando u=0. Espero con esto ayudar a que me ayude.
En este caso, el r integral en (1) es trivial y el integral tiene la forma:
IP=∫2π0dϕ∫π0dθ sinθ3x2y2cos2θ(y2cosϕ+x2sin2ϕ)sin2θ+x2y2cos2θ.
The ϕ integral can be integrated with the help of Eq. 3.642.1 in Gradstein and Ryzhik's tables of integrals. Thereby, the IP takes the for:
IP=3xy∫π0dθsinθcos2θ√1+(x2−1)cos2θ√1+(y2−1)cos2θ.
Now the change of variable k=√1−x2cosθ bring expression (6) to the form
IP=(const)x y(1−x2)3/2∫√1−x20dkk2√1−k2√1−k21−y21−x2.
Did you notice how IG and IP are similar? Do you think a similar approach can be applied to my original problem? Please, let me know.
Edit 2
The integral (1) también es evaluado en el Apéndice A. 4 de esta tesis. Sin embargo, no se utiliza simetría cilíndrica.