function axis_y(ylims)
xlims=get(gca,'Xlim');
axis([xlims ylims]);