function [ awo, fwo ] = getMaxAFwo( G, w1,w2 )
% Возвращает максимальные по модулю значения производных dA/dw, dF/dw 
% у функции G на интервале w1,w2
% знак в конце концов учитывается

w_step=(w2-w1)/30;
steps = uint32(w2/w_step+1);
X1=zeros(steps,1);
Y1=zeros(steps,1);
Y2=zeros(steps,1);

for w = w1:w_step:w2
    step = uint32(w/w_step+1);
    [mag1,phase1] = bode(G,w);
    [mag2,phase2] = bode(G,w+w_step);
    X1(step)=w;
    Y1(step)=(mag2-mag1)/(w_step);
    Y2(step)=(phase2-phase1)/57.3/(w_step);
end

[x1,y1]=get_max_val( X1, Y1, w1, w2 );
[x2,y2]=get_max_val( X1, Y2, w1, w2 );

awo = y1;
fwo = y2;

end