複素数とガンマ関数・ベータ分布

library(fAsianOptions)

cgamma(1i)

my.cbeta <- function(a,b){
	ret <- cgamma(a)*cgamma(b)/cgamma(a+b)
	return(ret)
}

my.dcbeta <- function(x,a,b){
	x^(a-1) * (1-x)^(b-1) / my.cbeta(a,b)
}

x <- seq(from=-1,to=2,length=100)
a <- -0.5
b <- -0.5
y <- my.dcbeta(x,a,b)

plot(x,y)