       Re: GenerateConditions->False gives fine result!

• To: mathgroup at smc.vnet.net
• Subject: [mg73976] Re: GenerateConditions->False gives fine result!
• From: "dimitris" <dimmechan at yahoo.com>
• Date: Sat, 3 Mar 2007 23:54:55 -0500 (EST)
• References: <es9238\$3ig\$1@smc.vnet.net>

```Hi and thanks a lot regarding answering to my messages/questions
apperared recently

> This is because the integral does in fact diverge in the Riemann sense. T=
aking the simpler example with {t->3, x->2}:

Now I did plots of the integrand I understand where was my mistake.
Thanks again!

In fact what confused me (which should have confused!) was the fact
that

In:=
(1/(2*Pi*I))*Integrate[q/(E^(2*I*q)*Sqrt[q^2 + 9]), {q, -Infinity,
Infinity}, GenerateConditions -> False]
N[%]
Out=
-((3*BesselK[1, 6])/Pi)
Out=
-0.001283348797177626

In:=
Chop[(1/(2*Pi*I))*NIntegrate[q*(Exp[(-I)*q*x]/Sqrt[q^2 + t^2]) /. {t -
> 3, x -> 2}, {q, -Infinity, Infinity},
Method -> Oscillatory]]
SequenceLimit::seqlim: The general form of the sequence could not be \
determined, and the result may be incorrect.
Out=
-0.0012833487971660698

I=2Ee. the two results are almost the same; this lead me to ignore the
messages from SequenceLimit!
Very amateurish mistake!

Anyway...Integrate rules! I need more reading!

Kind Regards
Dimitris

=CF/=C7 Bhuvanesh =DD=E3=F1=E1=F8=E5:
> This is because the integral does in fact diverge in the Riemann sense. T=
aking the simpler example with {t->3, x->2}:
>
> In:= integrand = q/(E^(2*I*q)*Sqrt[q^2 + 9]);
>
> In:= Limit[integrand, q->Infinity]
>
> Out= (1 + I) Interval[{-1, 1}]
>
> GenerateConditions->False, in addition to checking convergence and lookin=
g for singularities, also does Hadamard-type integrals. Here's another exam=
ple:
>
> In:= Integrate[1/(x^(1 + I/2)*(1 + x)), {x, 0, 1}]
>
>                                     1
> Integrate::idiv: Integral of ---------------- does not converge on {0, 1}.
>                               1 + I/2
>                              x        (1 + x)
>
>                          1
> Out= Integrate[----------------, {x, 0, 1}]
>                    1 + I/2
>                   x        (1 + x)
>
> In:= Integrate[1/(x^(1 + I/2)*(1 + x)), {x, 0, 1}, GenerateCondition=
s->False] //InputForm
>
> Out//InputForm= (-PolyGamma[0, -I/4] + PolyGamma[0, 1/2 - I/4])/2
>
> Bhuvanesh,
> Wolfram Research.

```

• Prev by Date: Re: nested list creation with index skipping
• Next by Date: Re: Problems with simplify and Sqrt[a^2]
• Previous by thread: Re: GenerateConditions->False gives fine result!