Problems with "integrate"

[yiwu ye <yiwu21111958@yahoo.com>]

Dear list,

I have an integration to do. The function is a gamma function as follows,

g_function(x,alpha.L,beta.L){1/(beta.L^alpha.L*gamma(alpha.L))*x^(alpha.L-1)*exp(-x/beta.L)}

No problems at all with the codes below:
beta_2
integrate(g,lower=0,upper=165,alpha.L=163/beta,beta.L=beta)
integrate(g,lower=0,upper=170,alpha.L=167.8/beta,beta.L=beta)

However, when beta=1.2, an error message says "singularity encountered..."
beta_1.2
integrate(g,lower=0,upper=165,alpha.L=163/beta,beta.L=beta)
integrate(g,lower=0,upper=170,alpha.L=167.8/beta,beta.L=beta)

Can anyone kindly help and explain why?

Thanks so much in advance,

Yiwu

---

Shape Yahoo! in your own image.  
Join our Network Research Panel today!