[R] help me to debug this part of code?

[francogrex]

jinjin wrote:
> 
> I am trying to solve the integration equation, for different values of K
> from 4 to 25, the integration is with respect to u,
> Here is the equation:    gamma(k/2) / ( sqrt(k-1)*gamma((k-1)/2) ) *
> integrate(f= (1+u^2/k-1)^(-k/2), lower=0, upper= sqrt(a^2*k/(k+1-a^2)) ) = 
> the similar expression as te left hand except k becomes k+1
> my code is below...
> 

I don't see that your expression at least as you wrote it above can be
integrated between the k values of 4 to 25, however if your k value is
negative, the expression (and the integral) has an analytical solution:
-((((a^2*k)/(1 - a^2 + k))^((1 - k)/2)*gamma(k/2))/-((-1 +
k)^1.5*(1/k)^(k/2)*gamma((-1 + k)/2)))
That's all I can help with, sorry. maybe someone with more mathematical
skills can help you out better.


-- 
View this message in context: http://www.nabble.com/help-me-to-debug-this-part-of-code--tp16763482p16799564.html
Sent from the R help mailing list archive at Nabble.com.