r Sample correlated bi-variate normal using R. Need expm package to compute sqrtm. Can compute derivatives automatically using pracma package.