Информатика / !704LAB2

uses crt;
const e=0.0001;
var x,y,s1,dx,x2,x1,t1,t2:real;
    i,n:integer;
    xi,yi:array[1..7] of real;
function f1(x:real):real;
begin
 f1:=0.67*x+19.325;
end;
function f2(x:real):real;
begin
 f2:=f1(x)-35*sqrt(1-sqr((x-13.01)/25))-1;
end;
function f4(x:real):real;
begin
 if x>10 then x:=10;
 f4:=-sqrt(400-sqr(x+10))-6.68+35*sqrt(1-sqr((x-13.01)/25))+1;
end;
function f5(x:real):real;
begin
 f5:=35*sqrt(1-sqr((x-13.01)/25))-1;
end;
begin
 clrscr;
 x:=-20;
 y:=0.67*x+19.325;
 repeat
  x:=x+e;
  y:=0.67*x+19.325;
 until sqr((x-13.01)/25)+sqr((y+1)/35)<=1;
 xi[1]:=x; yi[1]:=y;
 x:=-40;
 repeat
  x:=x+e;
  y:=0.67*x+19.325;
 until y>=0;
 xi[2]:=x; yi[2]:=0;
 x:=-20;
 repeat
  x:=x+e;
 until sqr((x-13.01)/25)+sqr(1/35)<1;
 xi[3]:=x; yi[3]:=0;
 x:=0;
 repeat
  x:=x+e;
 until sqr(x+10)+sqr(6.68)>=400;
 xi[4]:=x; yi[4]:=0;
 x:=20;
 repeat
  x:=x+e;
 until sqr((x-13.01)/25)+sqr(1/35)>=1;
 xi[5]:=x; yi[5]:=0;
 y:=-5;
 repeat
  y:=y-e;
 until sqr(10)+sqr(y+6.68)>=400;
 xi[6]:=0; yi[6]:=y;
 y:=-10;
 repeat
  y:=y-e;
 until sqr((-13.01)/25)+sqr((y+1)/35)>=1;
 xi[7]:=0; yi[7]:=y;
{/////////////////////////////////////////////////////////////////////////}
 writeln('                  ‚лзЁб«Ґ­Ёп ¤«п ®Ў« бвЁ Ђ                   ');
 writeln(' ----------------------------------------------------------- ');
 writeln('| —Ёб«® а §ЎЁҐ­Ё© |                 ђҐ§г«мв в               |');
 writeln('|                 |-----------------------------------------|');
 writeln('|        n        |      A      |      B      |      C      |');
 writeln('|-----------------|-----------------------------------------|');
 writeln('|       100       |             |             |             |');
 writeln('|      1000       |             |             |             |');
 writeln('|     10000       |             |             |             |');
 writeln(' ----------------------------------------------------------- ');
{1. ЊҐв®¤ Їаאַ㣮«м­ЁЄ®ў}
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   dx:=(xi[3]-xi[2])/n;
   x:=xi[2];
   repeat
    s1:=s1+dx*f1(x);
    x:=x+dx;
   until x>=xi[3];
   dx:=(xi[1]-xi[3])/n;
   x:=xi[3];
   repeat
    s1:=s1+dx*(f2(x));
    x:=x+dx;
   until x>=xi[1];
   gotoxy(23,6+i);
   write(s1:0:4,'  ');
   s1:=0;
  end;
{2. ЊҐв®¤ ва ЇҐжЁ©}
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   x2:=xi[3]; x1:=xi[2];
   dx:=(x2-x1)/n;
   x:=x1+dx;
   s1:=dx*(f1(x1)+f1(x2));
   repeat
    s1:=s1+dx*f1(x);
    x:=x+dx;
   until x>=x2-dx;
   x2:=xi[1]; x1:=xi[3];
   dx:=(x2-x1)/n;
   x:=x1+dx;
   s1:=s1+dx*(f2(x2)+f2(x1));
   repeat
    s1:=s1+dx*(f2(x));
    x:=x+dx;
   until x>=x2-dx;
   gotoxy(37,6+i);
   write(s1:0:4);
   s1:=0;
  end;
{3. ЊҐв®¤ ‘Ё¬Їб®­ }
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   x2:=xi[3]; x1:=xi[2];
   dx:=(x2-x1)/n;
   x:=x1+dx; t1:=0; t2:=0;
   repeat
    t1:=t1+f1(x);
    t2:=t2+f1(x-dx/2);
    x:=x+dx;
   until x>=x2-dx;
   t2:=t2+f1(x2-dx/2);
   s1:=dx/6*(f1(x1)+f1(x2)+2*t1+4*t2);
   x1:=xi[1]; x2:=xi[3];
   dx:=(x1-x2)/n;
   x:=x2+dx; t1:=0; t2:=0;
   repeat
    t1:=t1+f2(x);
    t2:=t2+f2(x-dx/2);
    x:=x+dx;
   until x>=x1-dx;
   t2:=t2+f2(x1-dx/2);
   s1:=s1+dx/6*(f2(x1)+f2(x2)+2*t1+4*t2);
   gotoxy(50,6+i);
   write(s1:0:4);
   s1:=0;
  end;
 gotoxy(1,12);
 writeln('                  ‚лзЁб«Ґ­Ёп ¤«п ®Ў« бвЁ B                   ');
 writeln(' ----------------------------------------------------------- ');
 writeln('| —Ёб«® а §ЎЁҐ­Ё© |                 ђҐ§г«мв в               |');
 writeln('|                 |-----------------------------------------|');
 writeln('|        n        |      A      |      B      |      C      |');
 writeln('|-----------------|-----------------------------------------|');
 writeln('|       100       |             |             |             |');
 writeln('|      1000       |             |             |             |');
 writeln('|     10000       |             |             |             |');
 writeln(' ----------------------------------------------------------- ');
{1. ЊҐв®¤ Їаאַ㣮«м­ЁЄ®ў}
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   dx:=10/n;
   x:=0;
   repeat
    s1:=s1+dx*f4(x);
    x:=x+dx;
   until x>=10;
   dx:=(xi[5]-10)/n;
   x:=10+dx;
   repeat
    s1:=s1+dx*f5(x);
    x:=x+dx;
   until x>=xi[5];
   gotoxy(23,17+i);
   s1:=s1-100;
   write(s1:0:4,'  ');
   s1:=0;
  end;
{2. ЊҐв®¤ ва ЇҐжЁ©}
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   dx:=10/n;
   x:=dx;
   s1:=dx*(f4(0)+f4(10));
   repeat
    s1:=s1+dx*f4(x);
    x:=x+dx;
   until x>=10-dx;
   dx:=(xi[5]-10)/n;
   x:=10+dx;
   s1:=dx*(f4(10)+f4(xi[5]));
   repeat
    s1:=s1+dx*f4(x);
    x:=x+dx;
   until x>=xi[5]-dx;
   gotoxy(37,17+i);
   write(s1:0:4);
   s1:=0;
  end;
{3. ЊҐв®¤ ‘Ё¬Їб®­ }
 s1:=0; n:=10;
 for i:=1 to 3 do
  begin
   n:=n*10;
   dx:=10/n;
   x:=dx; t1:=0; t2:=0;
   repeat
    t1:=t1+f4(x);
    t2:=t2+f4(x-dx/2);
    x:=x+dx;
   until x>=10-dx;
   t2:=t2+f4(10-dx/2);
   s1:=dx/6*(f4(0)+f4(10)+2*t1+4*t2);
   dx:=(xi[5]-10)/n;
   x:=10+dx; t1:=0; t2:=0;
   repeat
    t1:=t1+f4(x);
    t2:=t2+f4(x-dx/2);
    x:=x+dx;
   until x>=xi[5]-dx;
   t2:=t2+f4(xi[5]-dx/2);
   s1:=dx/6*(f4(10)+f4(xi[5])+2*t1+4*t2);
   gotoxy(50,17+i);
   write(s1:0:4);
   s1:=0;
  end;
 readln;
end.