[Date Index]
[Thread Index]
[Author Index]
Re: Beta function, Integral
*To*: mathgroup at smc.vnet.net
*Subject*: [mg79945] Re: Beta function, Integral
*From*: Asim <maa48 at columbia.edu>
*Date*: Thu, 9 Aug 2007 05:23:08 -0400 (EDT)
*References*: <f990ds$btq$1@smc.vnet.net><f9c1q6$75a$1@smc.vnet.net>
On Aug 8, 2:19 pm, dimitris <dimmec... at yahoo.com> wrote:
> On 7 , 08:37, Asim <ma... at columbia.edu> wrote:
>
> > Hi
>
> > The following integral does not seem to give the correct answer. The
> > answer should be the Euler Beta function, Beta[p,q]. Can anybody let
> > me know what I am doing wrong? Or is this a bug?
>
> > In[12]:= Integrate[t^{p - 1}*(1 - t)^(q - 1), {t, 0, 1}, Assumptions -
>
> > > {p > 0, q > 0}]
>
> > Out[12]= {(\[Pi] Csc[\[Pi] q] Gamma[p])/(Gamma[1 - q] Gamma[p + q])}
>
> > Thanks
>
> > Asim Ansari
>
> First note that you used List where you should have used parentheses!
> A common mistake.
> It must be t^(p-1); not {p-1}.
>
> In the Mathematica I work, I took:
>
> In[18]:=
> $Version
>
> Out[18]=
> "5.2 for Microsoft Windows (June 20, 2005)"
>
> In[20]:=
> Integrate[t^(p - 1)*(1 - t)^(q - 1), {t, 0, 1}, Assumptions -> {p > 0,
> q > 0}]
> FunctionExpand[Beta[p, q] - %]
>
> Out[20]=
> (Gamma[p]*Gamma[q])/Gamma[p + q]
>
> Out[21]=
> 0
>
> Ommiting { } from your output, we have also
>
> In[28]:=
> FullSimplify[Beta[p, q] - (Pi*Csc[Pi*q]*Gamma[p])/(Gamma[1 -
> q]*Gamma[p + q])]
>
> Out[28]=
> 0
>
> as it must be.
>
> Note also that
>
> In[30]:= Beta[p, q] // FunctionExpand
> Out[30]= (Gamma[p] Gamma[q])/Gamma[p + q]
>
> and
>
> In[46]:=
> FullSimplify[(Gamma[p]*Gamma[q])/Gamma[p + q] == Beta[p, q]]
>
> Out[46]=
> True
>
> but
>
> In[47]:=
> FullSimplify[c
>
> Out[47]=
> c
>
> Regards
> Dimitris
Thanks for all who responded to by Beta function integral post.
I used Mathematica 6.0, and by mistake used t^{p-1} in the integral,
instead of t^(p-1).
I was expecting (Gamma[p] Gamma[q])/Gamma[p + q] which can be
obtained from Mathematica 5.2 as the answer, but Mathematica 6.0
returns (\[Pi] Csc[\[Pi] q] Gamma[p])/(Gamma[1 - q] Gamma[p + q]). The
two are indeed the same, but the former is the more "standard" answer
given the definition of the Beta function.
Regards
Asim Ansari
Prev by Date:
**Re: FindRoot can NOT handle mixed real and complex variables**
Next by Date:
**Re: How show Locator point and its image under a function?**
Previous by thread:
**Re: Beta function, Integral**
Next by thread:
**Re: Beta function, Integral**
| |