function [ X1,Y1,X2,Y2,X3,Y3 ] = get_dAdz_arrays( rez )
%GET_DADZ_ARRAYS Summary of this function goes here
%   Detailed explanation goes here

mu = (18-b1*t_reg_min)/(b3*t_reg_min);
i=108/(b3*t_reg_min^2);
nu=216/(b3*t_reg_min^3);
dmu = mu/40;
di = i/40;
dnu = nu/40;

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);

for w = 0:w_step:1.6
    step = uint32(w/w_step+1);

    G1 = get_tf_zamkn(mu-dmu, i, nu, b1, b3);
    %G2 = get_tf_zamkn(mu, i, nu, b1, b3);
    G3 = get_tf_zamkn(mu+dmu, i, nu, b1, b3);
    [mag1,phase1] = bode(G1,w);
    [mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3,w);
    X1(step)=w;
    Y1(step)=(mag3-mag1)/(2*dmu);

    %%%%%%%%%%%%%
    G1 = get_tf_zamkn(mu, i-di, nu, b1, b3);
    %G2 = get_tf_zamkn(mu, i, nu, b1, b3);
    G3 = get_tf_zamkn(mu, i+di, nu, b1, b3);
    [mag1,phase1] = bode(G1,w);
    [mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3,w);
    X2(step)=w;
    Y2(step)=(mag3-mag1)/(2*di);

    %%%%%%%%%%%%%%
    G1 = get_tf_zamkn(mu, i, nu-dnu, b1, b3);
    %G2 = get_tf_zamkn(mu, i, nu, b1, b3);
    G3 = get_tf_zamkn(mu, i, nu+dnu, b1, b3);
    [mag1,phase1] = bode(G1,w);
    [mag2,phase2] = bode(G2,w);
    [mag3,phase3] = bode(G3,w);
    X3(step)=w;
    Y3(step)=(mag3-mag1)/(2*dnu);
end

end