Integrating a complicated expression involving Sign[...] etc.
- To: mathgroup at smc.vnet.net
- Subject: [mg56161] Integrating a complicated expression involving Sign[...] etc.
- From: Christian Mikkelsen <s010132 at student.dtu.dk>
- Date: Sun, 17 Apr 2005 03:07:07 -0400 (EDT)
- Organization: UNI-C
- Sender: owner-wri-mathgroup at wolfram.com
I want to ask some general questions. I get some very complicated
intermediate results (see expressions below) that Mathematica does not
proces (Integrate) within a day or two.
1) Does it normally pay off just leaving it to it? If Mathematica is still
running does it mean that it is doing meaningful transformations and
calculations?
2) I would like to understand the structure of some of the intermediate
results ("Thetaintegrant" in particular) to see if Mathematica might need
a little help (see my earlier post and the kludge below). Is there any
good way to do that? Using TreeForm, for example if I could just look at
the upper most levels of the expression tree?
3) My problem seems to have a lot of structure so I am hopeful that
something can be done to help Mathematica a little.
Oh, and a easy one... :-)
4) How do I use $Assumptions to tell Mathematica that A is real etc.
globally?
-----
(* This is actually a spherical Bessel-function *)
j[k_, n_] = ((-I)^n*Sqrt[2*Pi]*BesselJ[1/2 + n, k])/Sqrt[k]
(* A, B, C are to substituted by sines and cosines to make up polar
coordinates *)
StorMatrix = {{B^2+C^2, A B, A C},
{B A, A^2+C^2, B C},
{C A, C B, A^2+B^2}};
(* My integrant in the case (0, 0) *)
Integrant = Simplify[StorMatrix j[k A, n] j[k B, m] j[k a, 0] Exp[I k (A
x + B y + C z)]/.{n->0, m->0},
Assumptions -> {A \[Element] Reals, B \[Element] Reals,
C \[Element] Reals, a > 0,
x \[Element] Reals, y \[Element] Reals,
z \[Element] Reals}];
(* Mathematica refuses to do the definite integral but happily does the
indefinite one *)
KIntegral = Integrate[Integrant, k,
Assumptions -> {A \[Element] Reals, B \[Element] Reals,
C \[Element] Reals, a > 0,
x \[Element] Reals, y \[Element] Reals,
z \[Element] Reals}];
(* I calculate the lower limit, no problems here *)
Klowlimit = Limit[KIntegral, k -> 0,
Assumptions -> {A \[Element] Reals, B \[Element] Reals,
C \[Element] Reals, a > 0,
x \[Element] Reals, y \[Element] Reals,
z \[Element] Reals}]//Simplify;
(* Upper limit, with ComplexExpand and cheating. True? *)
KinfIntegral = ComplexExpand[KIntegral] /.
{ExpIntegralEi[ I A_ k] -> Sign[ I A] Pi,
ExpIntegralEi[-I A_ k] -> Sign[-I A] Pi};
Kinflimit = Limit[KinfIntegral, k -> \[Infinity],
Assumptions -> {A \[Element] Reals, B \[Element] Reals,
C \[Element] Reals, a > 0,
x \[Element] Reals, y \[Element] Reals,
z \[Element] Reals}];
Thetaintegrant = (Kinflimit-Klowlimit)/.{A->Sqrt[1-w^2] P, B->Sqrt[1-w^2]
Q, C-> w};
(* Warning: Mathematica 5.1 gets stuck on this integral : *)
Thetaintegral =
Integrate[Thetaintegrant, {w, -1, 1},
Assumptions -> {P >= -1, P <= 1, Q >= -1, Q <= 1, a > 0}]