function DVS
% -----------ХИМИЯ-----------------
O = 0;
C = 0.855;
H = 0.145;
% ---------------------------------
D = 14; % Диаметр поршня
Vc = 30; % Объём камеры сгорания
S = 15;
Pz = 3;
K = 1.67;
Tau = 2; %Число тактов рабочего цикла
i = 1; % Число цилиндров.
Vc = 5; % (мм) Объём камеры сгорания
Ro0 = 1; % Плотность воздуха
Hu = 44*10^6; % низкая тепло отводная способность топлива
Fi0 = 0.94; % Коэффициент продувки
T0 = 288 ; % (К)Температура окружающей среды.
Ta = 370 ; % (К)температура заряда в цилиндре в конце наполнения
Ty = 1050 ; % (К) температура остаточных газов
P0 = 0.101; % (МПа) - давление окр. среды.
Vh = (( pi*D^2)/4)*S;
Va = Vh+Vc;
Eps = 1+(Vh/Vc);
Pa = (0.7 / 0.95)* P0;
Pc = Pa * Eps^K;
L0 = (1/0.21)*((C/12)+(H/4)-(O/32));
Alpha = 21-(0.855*4*0.145-0.5*2*0.145);
Py = (1.05/1.2)*P0;
NUv = (T0/(P0*(Eps-1)))*((Pa*Eps/Ta)-Py/Ty);
Lam=Pz/Pc;
for k=1:25
   n(k) = 3750+250*k;
   Pm(k) = 0.035+0.00005*n(k);
   n2(k) = 1.22 - (130 / n(k));
   n1(k) = 1.41 - (100 / n(k));
   Ppi(k) = ((Pa*(Eps^n1(k)))/(Eps-1))*((Lam/(n2(k)-1))*(1-(1/(Eps^n2(k)-1)))-(1/(n1(k)-1))*(1-(1/Eps^(n1(k)-1))));
   Pi(k) = Ppi(k)*Fi0;
   Nui(k) = Pi(k)*(L0/Hu)*(Alpha/NUv)*(1/Ro0);
   Pe(k) = Pi(k)-Pm(k);
   Num(k) = Pe(k)/Pi(k);
   Gi(k) = 3600/(Nui(k)*Hu);
   Ne(k) = (i*Vh*Pe(k)*n(k))/(30*Tau);
   Ge(k) = Gi(k)*(1/Num(k));
   Gt(k) = Ge(k)*Ne(k);
   Me(k)= 159.1*(Ne(k)/n(k));
end
figure(1);
plot(Ne,n,'.-');
figure(2);
plot(Me,n,'.-');
figure(3);
plot(Gt,n,'.-');
disp([num2str(n)]);