MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Help in analytical integration 3

  • To: mathgroup at smc.vnet.net
  • Subject: [mg75202] Re: Help in analytical integration 3
  • From: dimitris <dimmechan at yahoo.com>
  • Date: Fri, 20 Apr 2007 00:40:02 -0400 (EDT)
  • References: <f04nr2$8nu$1@smc.vnet.net>

Searching more the G-R tables of Integrals I could not find any
relevant formula.
So I should have been more helpful as regards my previous replies.

I tried various methods in order to get a closed form result by myself
but I
failed. So I consulted another forum and I posted your integral; in
order to make the things
easier I replaced specific values for {A1,B1,C1,D1}={1,2,3,4}

That is, I took the following integral

In[24]:=
Integrate[1/Sqrt[(1 - w^2)*(4 - w^2)*(9 - w^2)*(16 - w^2)], {w, 0, 1}]

Mathematica can't return a result. Trying now the substitution w^2=v.
So that

In[32]:=
(FullSimplify[#1, v > 0] & )[Dt[w]/Sqrt[(1 - w^2)*(4 - w^2)*(9 -
w^2)*(16 - w^2)] /. w -> Sqrt[v] /. Dt[v] -> 1]
Integrate[%, {v, 0, 1}]

Out[32]=
1/(2*Sqrt[(-16 + v)*(-9 + v)*(-4 + v)*(-1 + v)*v])

Out[33]=
Integrate[1/(2*Sqrt[(-16 + v)*(-9 + v)*(-4 + v)*(-1 + v)*v]), {v, 0,
1}]

Again Mathematica fails.

Note that also the other CAS I frequently used does not exhibit better
performance in this case.
Again the integrals are remained unevaluated.

Noting also on the timings spent by Mathematica for this integral I am
quite
convinced that it does not how to deal with them. For integrals like
the above one
(I think that) table look up is taken place. Mathematica through
pattern matching
searches for these integrals but it does not find anything relevant
and it returns quickly the
integral back.

In the other forum in which I consulted to, professor G. A. Edgar
replied me (and I have many
reasons to trust him!) that may be, there is no answer more elementary
than Schwarz-Christoffel
transformation.

See

http://planetmath.org/?op=getobj&from=objects&id=6289
http://en.wikipedia.org/wiki/Schwarz-Christoffel_mapping
http://mathworld.wolfram.com/ConformalMapping.html

In our case (after the change of the integration variable!) we have
to find the coordinate of the vertex of the appropriate
hexagon on the positive real axis.

Enough for now!
I hope you found my current more helpful (and more...hopeful!)
than previous ones.

I will be in touch for anything new!

Regards
Dimitris

PS

1)

Don't wait something simple. Quite the opposite!
For example the much more simpler definite integral

In[41]:=
Integrate[1/Sqrt[(1 - x^2)*(2 - x^2)], {x, 0, 1}]

returns

Out[41]=
EllipticK[1/2]/Sqrt[2]

whereas of course

In[42]:=
(Information[#1, LongForm -> False] & )[EllipticK]

>From In[42]:=
"EllipticK[m] gives the complete elliptic integral of the first kind
K(m)."

2)

Think also of contour integration. Take this advice as a further hint.
I don't try it by myself.

3)

Don't forget that if the things don't go well, you have NIntegrate in
your hands!
Symbolic results is better from numerical ones (especially working in
applied mathematics
field). But also numerical results is better than no...results at all!






=CF/=C7 ashesh =DD=E3=F1=E1=F8=E5:
> Hi all,
>
> I need to solve the following integration analytically:
>
> Integrate[1/Sqrt[(A1^2 - w^2)*(B1^2 - w^2)*(C1^2 - w^2)*(D1^2 - w^2)],
> {w, 0, A1}]
>
>  0 < A1 < B1 < C1 < D1
>
> Mathematica is unable to perform the same.
>
> Any help / tips in solving the same will be highly appreciated.
>
> Regards.



  • Prev by Date: Re: Dirac Delta Function: Basics
  • Next by Date: Re: Importing and retaining graphics
  • Previous by thread: Updated BiokmodWeb, the web version of BIOKMOD
  • Next by thread: sorry for the previous massage