Select points with boundary

I have a matrix of 2D points like trajectory = [1,1; 1,2; 2,3; 3,4; 4,5; 5,5]

I want to limit it in a square of (0,0), (4,4), (0,4), (4,0).

Is there a way I can delete some rows out of the boundary and get the new trajectory = [1,1; 1,2; 2,3; 3,4]?

I know it is possible to do it in a lot of if statements. Is there a easy way to do that?

Answer 1

[Walter Roberson]

traj= [1,1; 1,2; 2,3; 3,4; 4,5; 5,5];
mask = traj(:,1) >= 0 & traj(:,1) <= 4 & traj(:,2) >= 0 & traj(:,2) <= 4;
new_traj = traj(mask,:);