matlab蒙特卡罗

%%求圆周率
sum = 0;
j = 0;
x = rand(1);
y = rand(1);
while j < 300000
    if x^2 + y^2 <= 1
       sum = sum + 1 ;
    end
    j = j + 1;
    x = rand(1);
    y = rand(1);
end
plot(x,y)
%%d单位圆面积 pi * 1 * 1  正方形面积 2*2
%%可得pi
p = 4 * sum / 300000