function [prob] = guessing()
j = 0;
for i = 1 : 10000
 x = ceil(3 * rand);
 y = ceil(3 * rand);
 j += (x-y) / x;
end
prob = j / 10000;