function [ X1,Y1,X2,Y2,X3,Y3 ] = get_dFdz_arrays( rez, t )
%GET_DADZ_ARRAYS Summary of this function goes here
%   Detailed explanation goes here
[ 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/100;
di = i/100;
dnu = nu/100;

w_max = 1.6;
w_step=0.1;
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_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);

    [mag1,phase1] = bode(G1_mu,w);
    %[mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3_mu,w);
    X1(step)=w;
    Y1(step)=(phase3-phase1)/(2*dmu);

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

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

    [mag1,phase1] = bode(G1_nu,w);
    %[mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3_nu,w);
    X3(step)=w;
    Y3(step)=(phase3-phase1)/(2*dnu);
    
    % градусы преобразуем в радианы
    Y1(step) = Y1(step)/57.3;
    Y2(step) = Y2(step)/57.3;
    Y3(step) = Y3(step)/57.3;
end

end