Add orthogonal range searching

I've extended this library for a project I'm doing and figured I'd upstream the changes.

It's nothing special - just some basic orthogonal range searching. Given two points, the function iteratively searches the tree and returns all nodes that fall in or on the hyperrectangle defined by the points.