clear
N = 100;
Mrep = 1000;
std( 4*mean( sum( rand(N,2,Mrep).^2, 2 )<1, 1 ) )