Dbvn(p,design=function(p,...) {
return(list(mu=cbind(p[1],p[1]),
dmu=cbind(1,1),
S=matrix(c(p[2],p[3],p[3],p[4]),ncol=2),
dS=rbind(c(1,0,0,0),c(0,1,1,0),c(0,0,0,1))) )},
Y=cbind(0,0))