c1=0.82;
c2=11.0;
c3=8.8;
c4=1.2;
c5=0.3;
c6=7.3;

OmegaA=Sqrt[c1*c4+c2]

TA=1/OmegaA

Tv=1/c4

Kc=c3*c4/(c1*c4+c2)

ash=2*(c1+c5-c4)/c3

bsh=((c1+c4+c5)^2-4*(c1*c4+c2))/(c3^2)

mu=-ash/2+Sqrt[(ash^2)/4-bsh]

Kcsh=Kc/(1+mu*Kc)

TAsh=TA/Sqrt[1+mu*Kc]

iv=1/(Kcsh*Tv)

nu=0.075/(Kcsh*Tv^2)

Tnu=iv/nu

MatrixForm[M1={{s+c1,c5*s+c2,0,c3},{-1,s+c4,0,0},{-1,0,s,0},{mu,0,nu*(Tnu*s+1)/s,-1}}]

MatrixForm[M2={{s+c1,c5*s+c2,0,c3},{-1,s+c4,0,0},{-1,0,0,0},{mu,0,nu*(Tnu*s+1)/s,-1}}]

W=Simplify[Det[M2]/Det[M1]]

h=Simplify[InverseLaplaceTransform[W/s,s,t]]

hf[t_]=Abs[h]

Plot[h,{t,0,50},PlotRange->{0,1.8}]

kol=10000;
dt=50/kol;

Ust=hf[10000];
hmax=hf[0];
For[i=1,i<kol+1,{If[Sign[Abs[hf[dt*(i-1)]-Ust]-0.05]≠Sign[Abs[hf[dt*(i)]-Ust]-0.05],Treg=dt*(i-0.5)],If[hmax<hf[dt*i],hmax=hf[dt*i]]};i++]
sigma=Round[100*(hmax-Ust)/Ust,0.01];
Print["Treg=",Treg,"  с."]
Print["Sigma=",sigma,"%."]

Soch=8 ;
koefstep=(1.1-0.9)/(Soch-2)

ocen=0;
VTsr={};
VSsr={};
For[i=1,i<Soch,
{For[j=1,j<Soch,
{For[k=1,k<Soch,
{For[l=1,l<Soch,
{For[m=1,m<Soch,
{c1c=c1*(0.9+koefstep*(i-1));
c2c=c2*(0.9+koefstep*(j-1));
c3c=c3*(0.9+koefstep*(k-1));
c4c=c4*(0.9+koefstep*(l-1));
c5c=c5*(0.9+koefstep*(m-1));
M1={{s+c1c,c5c*s+c2c,0,c3c},{-1,s+c4c,0,0},{-1,0,s,0},{mu,0,nu*(Tnu*s+1)/s,-1}};
M2={{s+c1c,c5c*s+c2c,0,c3c},{-1,s+c4c,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];
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];
If[Treg>5|| sigma>5,{Print["Режим с коэффициентами с1=",c1c,", c2=",c2c,", c3=",c3c,", c4=",c4c,", c5=",c5c," не удовлетворяет требованиям к переходному процессу, т.к. Treg=",Treg," c, Sigma=",sigma,"%."],ocen++,VTsr=Append[VTsr,Treg],VSsr=Append[VSsr,sigma]}]};
m++]};
l++]};
k++]};
j++]};
i++]
Ver=(ocen/((Soch-1)^5))*100.;
Tregsr=Mean[VTsr];
sigmasr=Mean[VSsr];
Print["Вероятность того, что система самолет – рассчитанный автопилот на данном режиме полета не будет удовлетворять заданным характеристикам переходного процесса равна ",Ver,"%."]
If[Ver>0,Print["При этом наиболее вероятное время регулирования – ",Tregsr,"с, а наиболее вероятное перерегулирование - ",sigmasr,"%. Для решения этой проблемы следует уменьшить число nu(передаточное число АП по сигналу интеграла угла тангажа), либо увеличить число iv(передаточное число АП по сигналу угла тангажа), либо уменьшать nu и увеличивать iv одновременно, но на меньшие величины, чем в первых двух случаях, т.е. таким образом сдвнуть Tnu ещё дальше от Tv - увеличить Tnu."]]