ylim returns the present axis limits as an array, so just use one return variable if you want to retrieve the existing limits (and, optionally, set new ones).
Looks like what would help more in this case would be to expand the x axis rather than y, though (altho both may help some)
xlim([0 10])
looks like would cover most of the initial transient period.
You might also want to investigate
doc zoom
