"r2norm" <-
function (N,m1=0,m2=0,sd1=1,sd2=1,rho=0) 
{
# Erzeugt N Zufallszahlenpaare mit zweidimensionaler Normalverteilung
# mit Parametern m1, m2, sd1, sd2, rho
Sigma<-matrix(c(sd1^2,rho*sd1*sd2,rho*sd1*sd2,sd2^2),ncol=2)
m<-c(m1,m2)
zmatrix<-matrix(rnorm(2*N),nrow=2)
A<-chol(Sigma)
r2norm<-t(A)%*%zmatrix +m
list(mu=m,Sigma=Sigma, r2norm=t(r2norm))
}