Лабораторные работы / Лабораторные работы - 2003 / programed

      Line 1     Col 1   Insert Indent          Unindent
program stat;
uses crt;
const
     np=9; nk=6; k=10;
var
    vm,v,ak:array[1..nk] of real; p:array[1..np] of real;
    hg,h:array[1..2] of real;
    a,b,c,e,ro,pn,g,x:real;
    i,kl:integer;
    bu:boolean;
FUNCTION SGN(ab:real):integer;
begin
             if ab>0 then sgn:=1;
             if ab<0 then sgn:=-1;
             if ab=0 then sgn:=0;
end;
FUNCTION FUNC (x:real):real;
label 300,400;
var  vm:array[1..6] of real;fx:real;
begin
        h[1]:=x;
        p[9]:=pn * hg[1] / (hg[1] - h[1]);
        p[6]:=p[9] + ro * g * h[1] *0.000001;
        v[1]:=ak[1] * SGN(p[1] - p[6]) * SQRt(ABS(p[1] - p[6]));
        v[3]:=ak[3] * SGN(p[6] - p[3]) * SQRT(ABS(p[6] - p[3]));
        v[6]:=v[1] - v[3];
        p[7]:=p[6] - SGN(v[6]) * sqr(v[6] / ak[6]);
        v[2]:=ak[2] * SGN(p[2] - p[7]) * SQRT(ABS(p[2] - p[7]));
        v[5]:=ak[5] * SGN(p[7] - p[5]) * SQRt(ABS(p[7] - p[5]));
        v[4]:=ak[4] * SGN(p[7] - p[4]) * SQRT(ABS(p[7] - p[4]));
        fx:=(v[2] + v[6] - v[4]-v[5]) * ro;
        FOR i:=1 TO 6 do
        vm[i]:=v[i] * ro;
        IF kl = 0 THEN FUNC:=fx;
        IF kl = 1 THEN writeln('x = ', x, ' fx = ', fx);
        end;
PROCEDURE MPD(a,b,eps:real; var bu:boolean; var xcon:real);
          var fa,fb,x,fx:real;
          label met;
          begin
          fa:=func(a);
          fb:=func(b);
          if fa*fb>0 then
          begin
          writeln(' РЕШЕНИЯ НЕТ ');  bu:=false; goto met
          end;
          repeat
            x:=(a+b)/2;  fx:=func(x);
            if fx*fa<0 then b:=x else a:=x;
          until abs(a-b)<eps;
          xcon:=abs(a+b)/2;  bu:=true;
          met:
          end;
   begin
  writeln('OTHOCuTELbHA9I POGRESHNOSTb e = (%)');
  read(e);
  writeln('BbICOTA EMOKOCTu 1(М) = ');
  read( hg[1]);
  writeln('BbICOTA EMKOCTu 2(M) = ');
  read( hg[2]);
  writeln('PLOTNOSTb ro = (кг/м3)');
  read(ro);
  writeln('NACHALNOE DAVLENIE pn = (МPа) ');
  read(pn);
  for i := 1 to 5 do
  begin
  writeln ('DAVLENIE p[',i,'] = ');
  read(p[i]);
  end;
  for i:=1 to 6 do
 begin
  writeln (' KOEFICIENT PROPUCKNOI CPOCOBNOCTI ak[',i,'] = ');
  read(ak[i]);
  end;
  writeln('VIVOD PROMEGHUTOCHNIH RESULTATOV : NET - 0  CHAST - 1  POLN
  read(kl);
  IF kl = 2 THEN writeln('    h         p(5-7)         vm');
  g:= 9.815;
  e:= e / 100;
  a:= 0;
  b:= hg[1] * (1 - e);
  MPD(a, b, e,bu, x);
  if bu then
 begin
  a:= ro * g * 0.000001;
  b:= p[7] + ro * g * hg[2] *0.000001;
  c:= (p[7] - pn) * hg[2];
  h[2]:= (b - SQRt(b *b - 4 * a * c)) / 2 / a;
  p[8]:= pn * hg[2] / (hg[2] - h[2]);
  FOR i := 1 TO 6 do
  vm[i]:= v[i] * ro;
  writeln('VIVOD REZULTATA');
  for i:=1 to 2 do
  writeln('VICOTA h[',i,'] = ',h[i]:k:k);
  writeln('DAVLENIE p[4] =',p[4]:k:K);
  for i:=7 to 9 do
  writeln('DAVLENIE p[',i,'] =',p[i]:k:k);
  for i:=1 to 6 do
  writeln('PACXOD vm[',i,'] =',vm[i]:k:k);
  end
 else
 begin
  kl:=2;
  writeln('a=',a:k:k);
  writeln('    h         p(5-7)          vm');
  x:=func(a);
  writeln('fa= ',x);
  writeln('b=',b);
  writeln('    h         p(5-7)          vm');
  x:=func(b);
  writeln('fb= ',x);
  end;
  repeat until  keypressed


END.