программирование в матлабе

[илья123456]

Всем привет. Есть такой код на матлабе, но у меня не выходит сделать так чтобы после нажатия кнопку Fractal, а потом Zoom и потом снова Fractal увеличивалась та область которую выделили. мое предположение что нужно исправить

 case 3
            set(gcf,'Units','normalized')
            k = waitforbuttonpress;
            rect_post = rbbox;
            %annotation('rectangle',rect_pos,'Color','red') /(16*xmax)
         
           xmax=rect_post(1)/xmax*xmin
           ymax=rect_post(2)/ymax*ymin
           xmin=-xmax
           ymin=-ymax

этот блок а именно

   xmax=rect_post(1)/xmax*xmin
           ymax=rect_post(2)/ymax*ymin
           xmin=-xmax-rect_post(3)
           ymin=-ymax-rect_post(4)

тут мои попытки

clc
clear
n=400; %размер окна
m=50; %количество итераций
colormap(hsv(m));
close all
xmin=-2;
xmax=2;
ymin=-2;
ymax=2;
while true
    switch menu ('Menu', 'Test','Fractal', 'Zoom','Reset','ZOOM','Quit')
        case 1
            c=4*(rand-.5)+2*1i*(rand-.5);
            n=50;
            y=rand(3);
            z=zeros(n,1);
            for i=2:n
                z(i)=z(i-1)^2+c;
            end
            i=0;
            plot(c)
            grid
            shg
            axis equal
            title(num2str(c))
            
        case 2
            if abs(xmax-xmin)>abs(ymax-ymin)
                d=(xmax-xmin)/n;
                %nx=n;
                %ny=(ymax-ymin)*nx/(xmax-xmin);
            else
                d=(ymax-ymin)/n;
               % ny=n;
              % nx=(xmax-xmin)*ny/(ymax-ymin);
            end
           
            [X,Y] = meshgrid(xmin:d:xmax,ymin:d:ymax);
            C=X+1i*Y;
            %zeros(size(C));
            Z=zeros(size(C));
            Q=Z;
            for i=1:m
                % Z[i]=Z(i-1)*Z(i-1)+C;
                Z=Z.^2+C;
                %isinf(Z);
                Q=Q+isinf(Z);
            end
            L=abs(Z)<inf;
            image(Q); shg;
            
            % plot(C(L),'.'), grid,
        case 3
            set(gcf,'Units','normalized')
            k = waitforbuttonpress;
            rect_post = rbbox;
            %annotation('rectangle',rect_pos,'Color','red') /(16*xmax)
         
           xmax=rect_post(1)/xmax*xmin
           ymax=rect_post(2)/ymax*ymin
           xmin=-xmax
           ymin=-ymax
               
        case 4
            xmin=-2;
            xmax=2;
            ymin=-2;
            ymax=2;
       
        case 5
              zoom 
        otherwise
            break
    end
end

кто знает или сталкивался с этим помогите)