function Q2 = zoomFFT_plot(fp,zp,nz,fmin,fmax,ymax); 
    fmax=fmin*(2^(1/3));
    choice=1;
    while(choice==1)
%   
        plot(fp,zp);
        xlabel('Frequency (Hz)');
        ylabel('Magnitude');
        out5 = sprintf(' ZOOM FFT  %d:1 ',nz);
        title(out5);
        grid;
        ymin=0.;
        ymax=max(zp)*1.2;
        iflag=0;
        for(kj=1:100)
            ne=kj-50;
            for(jq=1:100)
                limit=(jq/10)*(10^ne);
                if(ymax<limit)
                    ymax=limit;
                    iflag=1;
                    break;
                end
            end    
            if(iflag==1)
                break;
            end
        end
        axis([fmin,fmax,ymin,ymax]);
%
        disp(' ');
        disp('** See Figure 1 **');
        disp(' ');
        disp(' Change frequency limits on plot?  1=yes 2=no' );
        choice=input(' ');
        if(choice==1)
            disp(' ');
            disp(' Enter lower frequency ');
            fmin=input(' ');
            disp(' Enter upper frequency ');
            fmax=input(' ');   
        else
            break;
        end      
    end
Q2=ymax;  