clear; clc;

n = 1e6;
x = rand(1, n);
y = rand(1, n);
m = sum(x .^ 2 + y .^ 2 < 1);
disp(4 * m / n);

% 
% mask = x .^ 2 + y .^ 2 >= 1;
% plot(x(~mask), y(~mask), ".");
% hold on;
% plot(x(mask), y(mask), "r.");
    