Re: Beta function, Integral
- To: mathgroup at smc.vnet.net
- Subject: [mg79917] Re: Beta function, Integral
- From: dimitris <dimmechan at yahoo.com>
- Date: Thu, 9 Aug 2007 05:08:32 -0400 (EDT)
- References: <f990ds$btq$1@smc.vnet.net><f9c0i1$5rq$1@smc.vnet.net>
On 8 , 11:57, Jean-Marc Gulliet <jeanmarc.gull... at gmail.com> wrote: > Asim wrote: > > 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])} > > The /mathematical/ expression that is returned by Mathematica 6 is > correct, indeed. However, the answer has funny extra curly brackets and > if one try to use the expression to check its validity against Beta[p,q] > the expression returned unevaluated. On the other hand, if one type in > the expression by hand, the simplification occurs and the identity is > checked positively. > > In[1]:= $Version > > Out[1]= "6.0 for Microsoft Windows (32-bit) (June 19, 2007)" > > In[2]:= sol = > Integrate[t^{p - 1}*(1 - t)^(q - 1), {t, 0, 1}, > Assumptions -> {p > 0, q > 0}] > > Out[2]= {(\[Pi] Csc[\[Pi] q] Gamma[p])/(Gamma[1 - q] Gamma[p + q])} > --------^---------------------------------------------------------^ > Note the spurious curly brackets. > > In[3]:= FullSinplify[Beta[p, q] == sol[[1]]] > > Out[3]= FullSinplify[ > Beta[p, q] == (\[Pi] Csc[\[Pi] q] Gamma[p])/( > Gamma[1 - q] Gamma[p + q])] > > Even though we took out the contains of the list, Mathematica returns > the FullSimplify unevaluated. > > In[4]:= FullSimplify[ > Beta[p, q] == (Pi Csc[Pi q] Gamma[p])/(Gamma[1 - q] Gamma[p + q])] > > Out[4]= True > > Now, having entered the expression by hand, Mathematica is able to check > the identity. > > A similar behavior can be seen with Mathematica 5.2, though the > expression returned is different. > > In[1]:= > $Version > > Out[1]= > "5.2 for Microsoft Windows (June 20, 2005)" > > In[2]:= > sol = Integrate[t^{p - 1}*(1 - t)^(q - 1), {t, 0, 1}, > Assumptions -> {p > 0, q > 0}] > > Out[2]= > {(Gamma[p]*Gamma[q])/Gamma[p + q]} > > In[3]:= > FullSinplify[Beta[p, q] == sol[[1]]] > > Out[3]= > FullSinplify[Beta[p, q] == (Gamma[p]*Gamma[q])/ > Gamma[p + q]] > > In[4]:= > FullSimplify[Beta[p, q] == (Gamma[p]*Gamma[q])/ > Gamma[p + q]] > > Out[4]= > True > > -- > Jean-Marc Hi Jean Marc. I think the person who made the question should wanted to write Integrate[t^(p - 1)*(1 - t)^(q - 1), {t, 0, 1},Assumptions -> {p > 0, q > 0}] instead of Integrate[t^{p - 1}*(1 - t)^(q - 1), {t, 0, 1},Assumptions -> {p > 0, q > 0}] So that's explain the presence of the "spurious" curl brackets. Regards Dimitris