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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
