Re: Beta function, Integral
- To: mathgroup at smc.vnet.net
- Subject: [mg79914] Re: Beta function, Integral
- From: dimitris <dimmechan at yahoo.com>
- Date: Wed, 8 Aug 2007 05:00:54 -0400 (EDT)
- References: <f990ds$btq$1@smc.vnet.net>
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[(Gamma[p]*Gamma[q])/Gamma[p + q]] Out[47]= (Gamma[p]*Gamma[q])/Gamma[p + q] Regards Dimitris