"find" is not finding what I want it to?

I have a data set of 8341 numbers and I want to create two vector of the numbers between 3 and 200, and another of the numbers outside of that range. AI is the area included (between 3 and 200), and AE is the area excluded (below 3 or above 200) but when I run them they keep giving me the same amount (8341). Any idea what's going on?

I'm trying to plot the area included as a yellow ellipse, and the area excluded as a green ellipse.

 s = regionprops(BW, 'Area','Orientation', 'MajorAxisLength', ...
                    'MinorAxisLength','Area', 'Eccentricity', 'Centroid'); 
            hold on
            phi = linspace(0,2*pi,50);
            cosphi = cos(phi);
            sinphi = sin(phi);
            
            areas=[s.Area]
            AI=find(3<areas<200);
            AE=find(3.0>areas<200);
            for k = AI
                            xbar = s(k).Centroid(1);
                            ybar = s(k).Centroid(2);
                            a = s(k).MajorAxisLength/2;
                            b = s(k).MinorAxisLength/2;
                            theta = pi*s(k).Orientation/180;
                            R = [ cos(theta)   sin(theta)
                                 -sin(theta)   cos(theta)];
                            xy = [a*cosphi; b*sinphi];
                            xy = R*xy;
                            x = xy(1,:) + xbar;
                            y = xy(2,:) + ybar;
                            plot(x,y,'y','LineWidth',2);     
                            hold on;
            end
            hold on
                        for k = AE
                            xbar = s(k).Centroid(1);
                            ybar = s(k).Centroid(2);
                            a = s(k).MajorAxisLength/2;
                            b = s(k).MinorAxisLength/2;
                            theta = pi*s(k).Orientation/180;
                            R = [ cos(theta)   sin(theta)
                                 -sin(theta)   cos(theta)];
                            xy = [a*cosphi; b*sinphi];
                            xy = R*xy;
                            x = xy(1,:) + xbar;
                            y = xy(2,:) + ybar;
                            plot(x,y,'g','LineWidth',2);     
                            hold on;
            end 
           

Answer 1

[Image Analyst]

Try this:

logicalIndexes = areas > 3.0 & areas < 200;
AI=find(logicalIndexes); % Indexes within the area range.
AE= find(~logicalIndexes); % Indexes outside the area range.

Answer 2

[Chunru]

Try the following:

 AI=find(3.0<areas & areas<200);
 % AE=find(3.0>areas & areas<200);
 AE=find(3.0>areas | areas<200);

[Image Analyst]

AI=find(3.0<areas & areas<200);
AE=find(3.0>areas & areas<200);

I don't think your AE is right. See my answer below. Or use

AE=find(3.0>areas | areas>200);

[Chunru]

@Image Analyst You are right. I did not test the code as the data is not available. Edited.

[Kristin Aldridge]

@Image Analyst thank you for being so knowledgable, I feel like you answer all of my questions!

[Image Analyst]

@Kristin Aldridge, yes I did. Perhaps you overlooked it below. Maybe you could "Accept" it, or at least "Vote" for it.