InRect takes four subscript indices corresponding to rectangle vertices and the number of rows of the sampled image/matrix, and returns a row vector of linear indices corresponding to pixels that lie within the rectangle.

This simple algorithm is very useful for performing a fast convolution of a rectangle over an image.