Do you have to implement the gaussian copula yourself? Otherwise there is an open source statistical software package called R that has an implementation for the Gaussian, t and certain Archimedean copulas. This library is relatively easy to use. If R isnt installed, download R from
http://www.r-project.org/. When R is installed - install the copula library (Package -> install packages -> copula.zip) and then load the library. Often you are not interested in the copula itself but rather a probability distribution which implicitly has a certain copula dependence structure. When you have specified a dependence structure along with marginal distributions you have actual constructed a multivariate distribution function. In R you can simulate from this distribution but you also have access to the probability and density function.#R-codelibrary(copula)#Construction of the multivariate probability function with a gaussian copula (correlation = 0.3). The marginal distributions are a beta(5; 200) and a log-normal(-4; 0.5).thisCopula <- normalCopula(0.3, dim = 2, corstr = "un")this.mvdc <- mvdc(thisCopula, c("beta", "lnorm"), list(list(shape1 = 5, shape2 = 200), list(meanlog = -4, sdlog = 0.50)))#10 Random draws from the probability distribution with a Gaussian dependence structure. rmvdc(this.mvdc, 10) #Probability: P(marginal.1 < 0.05; marginal.2 < 0.09) pmvdc(this.mvdc, c(0.05,0.09)) #Density: dmvdc(this.mvdc, c(0.05,0.09)) #End R-codeHope this helps/Jonas