Re: a function containg MeijerG (limit behavior)
- To: mathgroup at smc.vnet.net
- Subject: [mg73329] Re: a function containg MeijerG (limit behavior)
- From: "dimitris" <dimmechan at yahoo.com>
- Date: Mon, 12 Feb 2007 05:06:23 -0500 (EST)
- References: <eqmche$bf4$1@smc.vnet.net>
> Sorry for the inconvenience.
No problem at all!
> The expression 'expr' above is what FunctionExpand[f[x1]] should have returned.
Thanks a lot firstly for your timing and secondly for reply!
"=2E..Enemy of the good is the better..."
I am sure next edition of Mathematica will not face this problem.
Cheers,
Dimitris Anagnostou
Research Associate
National University of Athens, Greece
=CF/=C7 Bhuvanesh =DD=E3=F1=E1=F8=E5:
> As I recall, there was a Sum issue in 5.2 that caused unexpected Interval results from MeijerG. So the Interval returned by the following input, for example, is incorrect.
>
> FunctionExpand[MeijerG[{{0}, {}}, {{0, 0}, {1/2}}, x]]
>
> For your example, after accounting for the FunctionExpand behavior, I get the following in version 5.2:
>
> In[1]:= expr = Cosh[2*x1]*(CoshIntegral[2*x1] + (Log[(-I)*x1] + Log[I*x1])/2 - Log[x1]) - Log[Abs[x1]] -
> Sinh[2*x1]*SinhIntegral[2*x1];
>
> In[2]:= Limit[expr, x1->0]
>
> Out[2]= EulerGamma + Log[2]
>
> In[3]:= Series[expr, {x1,0,2}, Assumptions -> 0<x1<1/100] //InputForm
>
> Out[3]//InputForm= SeriesData[x1, 0, {EulerGamma + Log[2], 0, -3 + 2*EulerGamma + 2*Log[2] + 2*Log[x1]}, 0, 3, 1]
>
> In[4]:= Limit[Normal@%, x1->0]
>
> Out[4]= EulerGamma + Log[2]
>
> In[5]:= Limit[expr, x1->Infinity]
>
> Out[5]= -Infinity
>
> The expression 'expr' above is what FunctionExpand[f[x1]] should have returned.
>
> Sorry for the inconvenience.
>
> Bhuvanesh,
> Wolfram Research.