       Re: N function problematic performance (V. 5.2)

• To: mathgroup at smc.vnet.net
• Subject: [mg74809] Re: N function problematic performance (V. 5.2)
• From: "dimitris" <dimmechan at yahoo.com>
• Date: Fri, 6 Apr 2007 04:19:11 -0400 (EDT)
• References: <euvm6e\$d7s\$1@smc.vnet.net><ev2av0\$kcp\$1@smc.vnet.net>

```But the question about the problematic behavior of the N function

(Warning messages

Infinity::"indet"...
NSum::nsnum :...

and failure to provide a numerical approximation in reasonable timing)

still holds of course!

Dimitris

=CF/=C7 dimitris =DD=E3=F1=E1=F8=E5:
> As I said somewhere else we have Mathematica on one side
> and on the other our mind. We can get the desirable results
> if we bring mind&Mathematica together...
>
> As we see using Integrate on Sin[x - 1]/Sqrt[x*(x - 2)] gives a result in=
terms
> of MeijerG functions. But then N faces problems to get a numerical
> estimation.
> On the hand applying TrigToExp first gives a result again in terms of Mei=
jerG
> functions which is opposite of the correct result.
>
> However a trivial subsitution results in a very simple result very
> quickly (bear in mind that
> in previous result containg MeijerG functions neither FullSimplify nor
> FunctionExpand
> provides any simplifications).
>
> So...
>
> (*Ins*)
> int = Simplify[(Sin[x - 1]/Sqrt[x*(x - 2)])*Dt[x] /. x -> t + 1] /.
> Dt[t] -> 1
> Plot[int, {t, 1, 10}]
> Timing[Integrate[int, {t, 1, Infinity}]]
> {N[%[],11], NIntegrate[int, {t, 1, Infinity}, Method ->
> Oscillatory,WorkingPrecision->20,PrecisionGoal->10]}
>
> (*Outs*)
> Sin[t]/Sqrt[-1 + t^2]
> Graphics
> {0.672*Second, (1/2)*Pi*BesselJ[0, 1]}
> {1.2019697153, 1.2019697153}
>
> Dimitris
>
> =CF/=C7 dimitris =DD=E3=F1=E1=F8=E5:
> > Consider the following the definite integral
> >
> > h = HoldForm[Integrate[Sin[x - 1]/Sqrt[x*(x - 2)], {x, 2, Infinity}]]
> >
> > Here is a numerical estimation
> >
> >
> > ReleaseHold[h /. Integrate[x___] :> NIntegrate[x, Method ->
> > Oscillatory]]
> > 1.2019697137297456
> >
> > Here is the symbolic result by Mathematica
> >
> > Timing[hh = ReleaseHold[h]]
> >
> > {30.031000000000002*Second, (Pi*(Cos*MeijerG[{{}, {1/4, 3/4}},
> {{0,
> > 1/2, 1/2}, {0}}, 1] -
> >      MeijerG[{{}, {1/4, 3/4}}, {{0, 0, 1/2}, {1/2}}, 1]*Sin))/
> > Sqrt}
> >
> > Application of N function results in unexpected warnings and not in
> > any estimation in reasonable timing...
> >
> > TimeConstrained[N[hh], 120]
> > Infinity::"indet"...
> > NSum::nsnum :...
> > \$Aborted
> >
> > What is going here?
> > Did I encouter a known (problematic) situation?
> >
> > BTW,
> >
> > Integrate[TrigToExp[Sin[x - 1]/Sqrt[x*(x - 2)]], {x, 2, Infinity}]
> > ((1/2)*I*Pi^(3/2)*(E^(2*I)*MeijerG[{{}, {1/2, 1/2}}, {{0, 0}, {1/2}},
> > -2*I] - MeijerG[{{}, {1/2, 1/2}}, {{0, 0}, {1/2}}, 2*I]))/
> >   E^I
> >
> > Chop[N[%]]
> > -1.2019697153172064
> >
> > which is opposite of the correct result.

```

• Prev by Date: Re: ftp from mathematica code
• Next by Date: Re: Enquirey
• Previous by thread: Re: N function problematic performance (V. 5.2)
• Next by thread: "cannot find the required native library named JLinkNativeLibrary" with 64-bit Linux and Java