function [ X1,Y1,X2,Y2,X3,Y3 ] = get_dAdz_arrays( rez, t )
% Возвращает численные значения графиков функций dA/dmu, dA/di, dA/dnu
% 1 - mu
% 2 - i
% 3 - nu

[ a1, a2, a3, a4, a5, a6, a7, b1, b2, b3, b4, b5, b6, b7 ]=make_koefs(rez);

mu = (18-b1*t)/(b3*t);
i=108/(b3*t^2);
nu=216/(b3*t^3);

dmu = mu/50;
di = i/50;
dnu = nu/50;

w_max = 1.6;
w_step=0.01;
steps = uint32(w_max/w_step+1);
X1=zeros(steps,1);
Y1=zeros(steps,1);
X2=zeros(steps,1);
Y2=zeros(steps,1);
X3=zeros(steps,1);
Y3=zeros(steps,1);

G1_mu = get_tf_kren_zamkn(mu-dmu, i, nu, b1, b3);
%G2 = get_tf_kren_zamkn(mu, i, nu, b1, b3);
G3_mu = get_tf_kren_zamkn(mu+dmu, i, nu, b1, b3);

G1_i = get_tf_kren_zamkn(mu, i-di, nu, b1, b3);
%G2 = get_tf_zamkn(mu, i, nu, b1, b3);
G3_i = get_tf_kren_zamkn(mu, i+di, nu, b1, b3);

G1_nu = get_tf_kren_zamkn(mu, i, nu-dnu, b1, b3);
%G2 = get_tf_zamkn(mu, i, nu, b1, b3);
G3_nu = get_tf_kren_zamkn(mu, i, nu+dnu, b1, b3);

for w = 0:w_step:w_max
    step = uint32(w/w_step+1);
    %fprintf('%f\n', step);
    [mag1,phase1] = bode(G1_mu,w);
    %[mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3_mu,w);
    X1(step)=w;
    Y1(step)=(mag3-mag1)/(2*dmu);

    %%%%%%%%%%%%%
    [mag1,phase1] = bode(G1_i,w);
    %[mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3_i,w);
    X2(step)=w;
    Y2(step)=(mag3-mag1)/(2*di);

    %%%%%%%%%%%%%%

    [mag1,phase1] = bode(G1_nu,w);
    %[mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3_nu,w);
    X3(step)=w;
    Y3(step)=(mag3-mag1)/(2*dnu);
    
    %Y1(step) = Y1(step)/15;
    %Y2(step) = Y2(step)/15;
    %Y3(step) = Y3(step)/15;
end

end