clc;
clear all;
close all;

x=unifrnd(0,100,50,1);
y=unifrnd(0,100,50,1);
for i=1 : 50
if x(i)>40 && x(i)<60
plot(x(i),y(i),'ro')
hold on;
else
plot(x(i),y(i),'o')
end
end