Abox = A(A(:,1)<X1 & A(:,1)>X2 & A(:,2)<Y1 & A(:,2)>Y2,:));

indices = find(A(:,1)<X1 & A(:,1)>X2 & A(:,2)<Y1 & A(:,2)>Y2);
Abox = A(indices,:);

Axrange = A(A(:,1)<X1 & A(:,1)>X2,:));
Abox = Axrange(Axrange(:,2)<Y1 & Axrange(:,2)>Y2,:);

find(A(:,1)<X1 & A(:,1)>X2 & A(:,2)<Y1 & A(:,2)>Y2)

indices = find(A(:,1)<X1 & A(:,1)>X2 & A(:,2)<Y1 & A(:,2)>Y2)
A2 = A(indices,:)