clear;

% --- Number of sampled points:

M = 10000;

% --- Main loop:

m = 0;

for i = 1:M;

% --- Sample point:

x = rand;
y = rand;

% --- Test and add to count:

if (x^2 + y^2 < 1.0);
m = m + 1;
endif;

endfor;

% --- Output:

printf("pi = %1.5f\n", 4*m/M);