function val=rbivnorm(m,V)
% simulates a single bivariate normal observation
% m is row vector [m1 m2], V is 2 x 2 matrix

a=chol(V);
val=randn(1,2)*a+m;