def get_bound(pts):
(x0, y0, x1, y1) = (INF, INF, -INF, -INF)<EOL>for (x, y) in pts:<EOL><INDENT>x0 = min(x0, x)<EOL>y0 = min(y0, y)<EOL>x1 = max(x1, x)<EOL>y1 = max(y1, y)<EOL><DEDENT>return (x0, y0, x1, y1)<EOL>
Compute a minimal rectangle that covers all the points.