esti = zeros(10000,1);

for j=1:1:10000

    y1 =  normrnd(0,sqrt(2),10000,1);
    y2 =  normrnd(0,sqrt(2),10000,1);


Z = y1./y2;

esti(j,1) = mean(Z);
end;

mean(esti)
var(esti)

y1so = sort(abs(y1));
y2so = sort(abs(y2));
Zso = sort(abs(Z));

prob = normpdf(1,0,sqrt(2)) - normpdf(-1,0,sqrt(2));