Uses CRT;
Const
 Xc=37.0;
 Yc=-39.0;
 Xe=-43.0;
 Ye=30.66;
 R=35.0;
 Ae=35.0;
 Be=25.0;

Var
 X0c,Y0c,X0e,Y0e,X1,Y1,X2,Y2,X3,Y3,Step:Real;
 i,Nt,k:Byte;
 N:Word;
 Rect,Trap,Simp:Array[1..100] of Real;
 Np:Array[1..100] of Word;

Function Ysl(Xsl:Real):Real;
begin
 Ysl:=Y1+(Xsl-X1)/(X2-X1)*(Y2-Y1);
end;

Function Ycr1(Xcr:Real):Real;
begin
 Ycr1:=Y0c+Sqrt(Sqr(R)-Sqr(Xcr-X0c));
end;

Function Ycr2(Xcr:Real):Real;
begin
 Ycr2:=Y0c-Sqrt(Sqr(R)-Sqr(Xcr-X0c));
end;

Function Yel(Xel:Real):Real;
begin
 Yel:=Y0e-Be*Sqrt(1-Sqr((Xel-X0e)/Ae));
end;

Function F1(Xf:Real):Real;
begin
 F1:=Ycr1(Xf)-Ysl(Xf);
end;

Function F2(Xf:Real):Real;
begin
 F2:=Yel(Xf)-Ycr2(Xf);
end;

Function F3(Xf:Real):Real;
begin
 F3:=Ysl(Xf)-Ycr2(Xf);
end;

Function Rectangle1:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X1;
 b:=0;
 h:=(b-a)/N;
 Sum:=0;
 for i:=0 to N-1 do
  Sum:=Sum+h*F1(a+h*i);
 Rectangle1:=Sum;
end;

Function Rectangle2:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=0;
 b:=X2;
 h:=(b-a)/N;
 Sum:=0;
 for i:=0 to N-1 do
  Sum:=Sum+h*F2(a+h*i);
 Rectangle2:=Sum;
end;

Function Rectangle3:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X2;
 b:=X3;
 h:=(b-a)/N;
 Sum:=0;
 for i:=0 to N-1 do
  Sum:=Sum+h*F3(a+h*i);
 Rectangle3:=Sum;
end;

Function Trapezium1:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X1;
 b:=0;
 h:=(b-a)/N;
 Sum:=(F1(a)+F1(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F1(a+h*i);
 Trapezium1:=Sum;
end;

Function Trapezium2:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=0;
 b:=X2;
 h:=(b-a)/N;
 Sum:=(F2(a)+F2(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F2(a+h*i);
 Trapezium2:=Sum;
end;

Function Trapezium3:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X2;
 b:=X3;
 h:=(b-a)/N;
 Sum:=(F3(a)+F3(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F3(a+h*i);
 Trapezium3:=Sum;
end;

Function Simpson1:Real;
var
 i,Ns:Integer;
 Sum,h,a,b:Real;
begin
 a:=X1;
 b:=0;
 Ns:=N div 2;
 h:=(b-a)/(2*Ns);
 Sum:=(F1(a)+F1(b))/3*h;
 for i:=1 to Ns do
  Sum:=Sum+4.0*h/3.0*F1(a+h*(2*i-1));
 for i:=2 to Ns do
  Sum:=Sum+2.0*h/3.0*F1(a+h*(2*i-2));
 Simpson1:=Sum;
end;

Function Simpson2:Real;
var
 i,Ns:Integer;
 Sum,h,a,b:Real;
begin
 a:=0;
 b:=X2;
 Ns:=N div 2;
 h:=(b-a)/(2*Ns);
 Sum:=(F2(a)+F2(b))/3*h;
 for i:=1 to Ns do
  Sum:=Sum+4.0*h/3.0*F2(a+h*(2*i-1));
 for i:=2 to Ns do
  Sum:=Sum+2.0*h/3.0*F2(a+h*(2*i-2));
 Simpson2:=Sum;
end;

Function Simpson3:Real;
var
 i,Ns:Integer;
 Sum,h,a,b:Real;
begin
 a:=X2;
 b:=X3;
 Ns:=N div 2;
 h:=(b-a)/(2*Ns);
 Sum:=(F3(a)+F3(b))/3*h;
 for i:=1 to Ns do
  Sum:=Sum+4.0*h/3.0*F3(a+h*(2*i-1));
 for i:=2 to Ns do
  Sum:=Sum+2.0*h/3.0*F3(a+h*(2*i-2));
 Simpson3:=Sum;
end;

BEGIN
 ClrScr;
 X0c:=2.0;
 Y0c:=Yc+Sqrt(sqr(R)-sqr(X0c));
 Y0e:=-5.66;
 X0e:=-8;
 Writeln('X0c=',X0c:9:5,' - X Є®®а¤Ё**в* жҐ*ва* ®Єаг¦*®бвЁ');
 Writeln('Y0c=',Y0c:9:5,' - Y Є®®а¤Ё**в* жҐ*ва* ®Єаг¦*®бвЁ');
 Writeln('X0e=',X0e:9:5,' - X Є®®а¤Ё**в* жҐ*ва* н««ЁЇб*');
 Writeln('Y0e=',Y0e:9:5,' - Y Є®®а¤Ё**в* жҐ*ва* н««ЁЇб*');
 Writeln;
 Step:=1;
 X1:=X0e-Ae;
 For i:=1 to 5 do
  begin
   Step:=Step/10;
   Repeat
    X1:=X1+Step;
    Y1:=Y0e+Be*Sqrt(1-sqr((X1-X0e)/Ae));
   Until Sqr(X1-X0c)+Sqr(Y1-Y0c)>Sqr(R);
   X1:=X1-Step;
  end;
 Y1:=Y0e+Be*Sqrt(1-sqr((X1-X0e)/Ae));
 Writeln('X1= ',X1:9:5,' - X Є®®а¤Ё**в* ЇҐаў®© в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Writeln('Y1= ',Y1:9:5,' - Y Є®®а¤Ё**в* ЇҐаў®© в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Y2:=0;
 X2:=X0e+Ae*Sqrt(1-sqr(Y0e/Be));
 Writeln('X2= ',X2:9:5,' - X Є®®а¤Ё**в* ўв®а®© в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Writeln('Y2= ',Y2:9:5,' - Y Є®®а¤Ё**в* ўв®а®© в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Step:=1;
 X3:=X2;
 For i:=1 to 5 do
  begin
   Step:=Step/10;
   Repeat
    X3:=X3+Step;
    Y3:=Ysl(X3);
   Until Sqr(X3-X0c)+Sqr(Y3-Y0c)>Sqr(R);
   X3:=X3-Step;
  end;
 Y3:=Ysl(X3);
 Writeln('X3= ',X3:9:5,' - X Є®®а¤Ё**в* ваҐв쥩 в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Writeln('Y3= ',Y3:9:5,' - Y Є®®а¤Ё**в* ваҐв쥩 в®зЄЁ ЇҐаҐбҐзҐ*Ёпя');
 Writeln;
 Writeln('НЌ*¦¬ЁвҐ "Enter" ¤«п Їа®¤®«¦Ґ*Ёп...');
 Readln;
 ClrScr;
 k:=0;
 Repeat
  Write('‚ўҐ¤ЁвҐ Є®«ЁзҐбвў® в®зҐЄ (0 ¤«п ўл室* Ё§ Їа®Ја*¬¬л) ');
  Readln(N);
  If N>0 then
   Begin
    ClrScr;
    Inc(k);
    Np[k]:=N;
    For i:=1 to 58 do Write('-');
    Writeln;
    Writeln('| —Ёб«® а*§ЎЁҐ*Ё© |              ђҐ§г«мв*в               |');
    Writeln('|                 |--------------------------------------|');
    Writeln('|       n         |      A     |      B     |     C      |');
    For i:=1 to 58 do Write('-');
    Writeln;
    Rect[k]:=Rectangle1+Rectangle2+Rectangle3;
    Trap[k]:=Trapezium1+Trapezium2+Trapezium3;
    Simp[k]:=Simpson1+Simpson2+Simpson3;
    For i:=1 to k do
     Writeln('|     ',Np[i]:4,'        |',Rect[i]:12:6,'|',Trap[i]:12:6,'|',Simp[i]:12:6,'|');
    For i:=1 to 58 do Write('-');
    Writeln;
   End;
 Until N=0;
END.