[latlim,lonlim]
= geoquadpt(lat,lon) returns
the limits of the tightest possible geographic quadrangle that bounds
a set of points with geographic coordinates lat and lon.

In most cases, tf = ingeoquad(lat,lon,latlim,lonlim)
will return true, but tf may be false for points
on the edges of the quadrangle, due to round off. tf will
also be false for elements of lat that fall outside
the interval [-90 90] and elements of lon that
are not finite.

In this case the output quadrangle straddles
the 180-degree meridian, hence the elements of lonlim are
in descending numerical order, although they are ordered from west
to east.

Read a set of points and then create a bounding box around
the points.

Latitude limits of a geographic quadrangle, returned as a 1-by-2
vector of the form [southern_limit northern_limit],
in units of degrees. The elements are in ascending order, and both
lie in the closed interval [-90 90].

Longitude limits of a geographic quadrangle, returned as a 1-by-2
vector of the form [western_limit eastern_limit],
in units of degrees. The limits are wrapped to the interval [-180
180]. They are not necessarily in numerical ascending order.