ch = get (gca, "children"); lims = NaN (2, 2); for ii=1:numel (ch) lims(1, 1) = min (lims(1, 1), min (get (ch(ii), "xdata"))); lims(1, 2) = max (lims(1, 2), max (get (ch(ii), "xdata"))); lims(2, 1) = min (lims(2, 1), min (get (ch(ii), "ydata"))); lims(2, 2) = max (lims(2, 2), max (get (ch(ii), "ydata"))); endfor ## Print [xmin xmax; ymin ymax] lims