Dear all,

I am trying to make a numerical integral of a Hypergeometrical function 2F1
but it is not working. I dont know how to pass the arguments to my
function. Thank you beforehand.
Chuse.

Here is my code:

H <- 0.7;

t <- 1;
s <- seq(0,t,0.1);
x <- 1-t/s;

library(gsl)
Gauss2F1 <- function(H,x){
  a <- H-1/2;
  b <- 1/2-H;
  c <- H+1/2;
#  if(x>=0 & x<1){
#    hyperg_2F1(a,b,c,x)
#  }else{
    hyperg_2F1(c-a,b,c,1-1/(1-x))/(1-x)^b
#  }
}

F <- Gauss2F1(a,b,c,x)

integrate(F,lower=-1,upper=1,H=0.7)

	[[alternative HTML version deleted]]

