R,W,VM,B,FI,C,NG,NM    :Array[1..jmax] of single;

   Q,H,SV,NS,OK,QK,GK,PN,A,DVAL,DVT,V01,D0,V0,R1,B1,U1,HT,

   VU2,U21,R21,VM1,BET1,BETL1,VM21,A0,BET2,Z1,HTB,VM2,U2,

   R2,PSI2,PSI1,W1,W2,DVM,DW,DR,BETA,DFI,z,jn:single;

   J,N    :Integer;

Begin

   Assign(F00,'KRREZ');

   Assign(F01,'CON');

   Assign(f0,'nastcen.dat'); Rewrite(f0);

    Assign(f0isl,'isl.gr2'); Rewrite(f0isl);

   Rewrite(F00); Rewrite(F01);

   (*Ввод  исходных данных*)

   QN:=180;         (*Подача насоса, м3/час*)

   HN:=340;          (*Напор насоса, м*)

   NO:=1475;        (*Частота вращения, об/мин*)

   RO:=1000;        (*Плотность жидкости, кг/м3*)

   KP:=1;           (*Число потоков*)

   ISt:=8;           (*Число ступеней*)

   S:=1;            (*Толщина лопастей,мм*)

   TAU:=1500;       (*Допустимые напряжения кручения, н/см2*)

   EK:=0.96;        (*Механический кпд*)

   DELTA:=0.166;  (*Угол атаки, рад*)

   PSI10:=0.73;     (*Коэф. стеснения на входе*)

   PSI20:=0.7973;   (*Коэф. стеснения на выходе*)

   K1:=0.06;     K2:=0.8;      K3:=1.0;    K4:=0.55;

   jn:=jmax*1.0;

   Q:=QN/KP/3600;

   H:=HN/ISt;

   SV:=3.14*NO/30;

40:   NS:=3.65*NO*Sqrt(Q)/Exp(0.75*Ln(H));

     If NS > 140 then

       Begin

         Writeln(F01,' Коэффициент быстроходности больше 140');

         Goto 1;

     End;

   OK:=1/(1+0.68*Exp(-0.666*Ln(NS)));

   QK:=Q/OK;

   GK:=1-0.42/Sqr(0.43429*Ln(4500*Exp(0.333*Ln(Q/NO)))-0.172);

   PN:=RO*9.81*QN*HN/1000/OK/GK/EK/3600;

   A:=Exp(0.333*Ln(16*975000*1.1*PN/NO/3.14/TAU))*10;

   DVAL:=(Int(A/5)+1)*5;

   DVT:=Int(DVAL*1.5);

   V01:=K1*Exp(0.333*Ln(QK*Sqr(NO)));