testDensityMultiNormal function

3d Mutivariate Normal likelihood