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

% что получается:
[ X1,Y1,X2,Y2,X3,Y3 ] = get_dAdz_arrays( 2, 2 );
plot(X1,Y1,X2,Y2,X3,Y3)
%grid on


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

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

G1 = get_tf_kren_zamkn(mu, i, nu, b1, b3);
%G2 = get_tf_kren_zamkn(mu, i, nu, b1, b3);
G2 = get_tf_kren_zamkn(mu+dmu, i, nu, b1, b3);

P = bodeoptions; % Set phase visiblity to off and frequency units to Hz in options
P.MagUnits = 'abs';
P.MagScale = 'linear';
P.FreqScale = 'linear';
P.Xlim = [0, 1.8];
P.Ylim = {[0, 1],[-90, 0]};
P.MagLowerLimMode = 'manual';
P.MagLowerLim = -100;
P.Grid = 'off';
% АФЧХ при Тмин и Тмакс
%h = bodeplot(G1,'r',G2,'b--',P)
bodeplot(G1,'r',G2,'b--',P)

w=0.8607;
[mag1,phase1] = bode(G1,w)
[mag2,phase2] = bode(G2,w)

dmag=mag2-mag1