Freqz точно в указанном axes

function Calculate
global hAxes hEd1 hEd2 hEd3 hEd4 hEd5
cla;
str1 = get(hEd1,'string'); 
str2 = get(hEd2,'string'); 
str3 = get(hEd3,'string'); 
str4 = get(hEd4,'string'); 
str5 = get(hEd5,'string'); 
Wn1 = str2num(str1); 
Wn2 = str2num(str2);
Rp = str2num(str3); 
Rs = str2num(str4); 
n  = str2num(str5);
Fs = Wn2*2.5;
Ws1 = Wn1-50;
Ws2 = Wn2+50;
Wn = [Wn1 Wn2]/Fs*2; 
Ws = [Ws1 Ws2]/Fs*2;
nfft = 512;                      % Число точек характеристик
[n,Wn] = ellipord(Wn,Ws,Rp,Rs);
%[b1,a1]  = butter(n,Wn);      
%[b2,a2]  = cheby1(n,Rp,Wn);
%[b3,a3]  = cheby2(n,Rs,Wn);
[b4,a4]  = ellip(n,Rp,Rs,Wn);    
[h,f,s] = freqz(b4,a4,nfft,Fs); % Частотная характеристика первого фильтра 
                                 % Далее используются те же значения nfft и Fs, поэтому нам 
                                 % не нужно второй раз получать выходные параметры f и s
%h2       = freqz(b2,a2,nfft,Fs); % Частотная характеристика второго фильтра 
%h3       = freqz(b3,a3,nfft,Fs);
%h       = freqz(b4,a4,nfft,Fs);
%h = [h1 h2 h3 h4];                     % Матрица из частотных характеристик двух фильтров
freqzplot(h,f,'s')           % Выводим график квадрата амплитуды
s.plot   = 'mag';                % Только график АЧХ
s.xunits = 'Hz';                % Частотная ось проградцирована в килогерцах
s.yunits = 'dB';            % Выводится АЧХ, в дБ
freqzplot(h,f,s);                % Сравниваем характеристики 4х фильтров
hAxes(axis);