clc
clear
area_cons = 10000; % Area under consideration
L = 50; % Consider the square with boundaries [-L L] x [-L L] with are 10000
lambda = 10^(-5); % Intensity of the Poisson point process
lambda_dash = lambda.*area_cons; % Average number of users in the considered area
n_point = poissrnd(lambda_dash); % Number of points, Poisson distributed
x = 2*L*rand(1,n_point)-L; % Uniformly distributing x-coordinates of points
y = 2*L*rand(1,n_point)-L; % Uniformly distributing y-coordinates of points
scatter(x,y)
A word of caution: The intensity considered is very small. For the intensity of 10^(-5) and area of 10000, you will get only 0.1 point, so almost no point. Try changing the intensity to somewhat higher value, e.g., 10^(-3), to see a realization.
