c1 = {0.560, 0.560, 0.496, 0.496, 0.670, 0.670, 0.590, 0.590, 1.100, 1.100, 0.940, 0.940, 0.760, 0.760, 0.880, 0.880, 0.668, 0.668, 0.793, 0.793}

c2 = {1.635, 0.922, 1.273, 0.718, 1.864, 0.655, 1.200, 0.245, 1.864, 0.569, 2.259, 0.0726, 2.640, 0.357, 3.337, 0.404, 1.742, 0.056, 2.980, 0.681}

c3 = {1.468, 1.468, 1.142, 1.142, 2.570, 2.570, 2.071, 2.071, 7.827, 7.827, 5.975, 5.975, 5.617, 5.617, 7.028, 7.028, 4.496, 4.496, 6.340, 6.340}

c4 = {0.765, 0.765, 0.868, 0.868, 1.053, 1.053, 0.944, 0.944, 1.122, 1.122, 0.962, 0.962, 0.864, 0.864, 0.987, 0.987, 0.671, 0.671, 0.875, 0.875}

c5 = {0.160, 0.160, 0.168, 0.168, 0.195, 0.195, 0.107, 0.107, 0.314, 0.314, 0.269, 0.269, 0.226, 0.226, 0.251, 0.251, 0.191, 0.191, 0.226, 0.226}

Mah = {0.167, 0.167, 0.147, 0.147, 0.209, 0.209, 0.184, 0.184, 0.42, 0.42, 0.35, 0.35, 0.43, 0.43, 0.478, 0.478, 0.4, 0.4, 0.5, 0.5}

g1={}

g2 = Plot[1.05,{t,0,20},PlotRange->{0,1.2},PlotStyle->Dashed];
g3 = Plot[0.95,{t,0,20},PlotRange->{0,1.2},PlotStyle->Dashed];
g4 = Plot[1,{t,0,20},PlotRange->{0,1.2}];

kol=10000

dt=20/kol

param = 1;
MTreg = {};
Msigma = {};
Mmu={};
Miv = {};
Mnu = {};
For [rezim=1, rezim<21,rezim = rezim + param,
{OmegaA=Sqrt[c1[[rezim]]*c4[[rezim]]+c2[[rezim]]]; 
TA=1/OmegaA;
Tv=1/c4[[rezim]];
Kc=c3[[rezim]]*c4[[rezim]]/(c1[[rezim]]*c4[[rezim]]+c2[[rezim]]);
ash=2*(c1[[rezim]]+c5[[rezim]]-c4[[rezim]])/c3[[rezim]];
bsh=((c1[[rezim]]+c4[[rezim]]+c5[[rezim]])^2-4*(c1[[rezim]]*c4[[rezim]]+c2[[rezim]]))/((c3[[rezim]])^2);
mu=4*Abs[-ash/2+Sqrt[(ash^2)/4-bsh]];
Kcsh=Kc/(1+mu*Kc);
TAsh=TA/Sqrt[1+mu*Kc];
If[Tv<(10*TAsh),{iv=1.2/(Kcsh*Tv); nu=0.05/(Kcsh*Tv^2)},{iv=0.5/(Kcsh*Tv); nu=5/(Kcsh*Tv^2 )}];
Tnu=iv/nu;
M1={{s+c1[[rezim]],c5[[rezim]]*s+c2[[rezim]],0,c3[[rezim]]},{-1,s+c4[[rezim]],0,0},{-1,0,s,0},{mu,0,nu*(Tnu*s+1)/s,-1}} ;
M2={{s+c1[[rezim]],c5[[rezim]]*s+c2[[rezim]],0,c3[[rezim]]},{-1,s+c4[[rezim]],0,0},{-1,0,0,0},{mu,0,nu*(Tnu*s+1)/s,-1}};
W=Det[M2]/Det[M1];
h=InverseLaplaceTransform[W/s,s,t];
hf[t_]=Abs[h];
g1 = Append[g1,Plot[h,{t,0,20},PlotRange->{0,1.2}]];
Ust=hf[10000];
hmax=hf[0];
For[q=1,q<kol+1,{If[Sign[Abs[hf[dt*(q-1)]-Ust]-0.05]≠Sign[Abs[hf[dt*(q)]-Ust]-0.05],Treg=dt*(q-0.5)],If[hmax<hf[dt*q],hmax=hf[dt*q]]};q++];
sigma=Round[100*(hmax-Ust)/Ust,0.01];
MTreg = Append[MTreg, Treg];
Msigma = Append[Msigma, sigma];
Mmu = Append[Mmu, mu];
Miv = Append[Miv, iv];
Mnu = Append[Mnu, nu];
Print["Режим №", rezim, ": tрег = ", Treg, "с, sigma = ", sigma, "%."]}]

Show[g1, g2, g3, g4]

r1 = Table[{Mah[[i]],Mmu[[i]]},{i,20}];
r2 = Table[{Mah[[i]],Miv[[i]]},{i,20}];
r3 = Table[{Mah[[i]],Mnu[[i]]},{i,20}];

gm1 =ListLinePlot[r1];
gm2 =ListLinePlot[r2];
gm3 =ListLinePlot[r3];

Mmu

Miv

Mnu