powered by
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))