function [M, xy, id] = PuntosDentroDelCirculo(N)

xy = GenerarNumerosRand(N);

id = sum(xy.^2) < 1;
M = sum(id);