Uses CRT;
Const
 R=25.0;
 Ae=16.0;
 Be=37.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 Ycr(Xcr:Real):Real;
begin
 Ycr:=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:=-Ycr(Xf);
end;

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

Function F3(Xf:Real):Real;
begin
 F3:=Yel(Xf)-Ycr(Xf);
end;

Function F4(Xf:Real):Real;
begin
 F4:=Ysl(Xf)-Ycr(Xf);
end;

Function Rectangle1:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X0c-Sqrt(sqr(R)-sqr(Y0c));
 b:=X1;
 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:=X1;
 b:=0;
 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:=0;
 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 Rectangle4:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X3;
 b:=X2;
 h:=(b-a)/N;
 Sum:=0;
 for i:=0 to N-1 do
  Sum:=Sum+h*F4(a+h*i);
 Rectangle4:=Sum;
end;

Function Trapezium1:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X0c-Sqrt(sqr(R)-sqr(Y0c));
 b:=X1;
 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:=X1;
 b:=0;
 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:=0;
 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 Trapezium4:Real;
var
 i:Integer;
 Sum,h,a,b:Real;
begin
 a:=X3;
 b:=X2;
 h:=(b-a)/N;
 Sum:=(F4(a)+F4(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F4(a+h*i);
 Trapezium4:=Sum;
end;

Function Simpson1:Real;
var
 i,Ns:Integer;
 Sum,h,a,b:Real;
begin
 a:=X0c-Sqrt(sqr(R)-sqr(Y0c));
 b:=X1;
 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:=X1;
 b:=0;
 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:=0;
 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;

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

BEGIN
 ClrScr;
 X0c:=-4.0; Y0c:=-15.0;
 X0e:=-7.0; Y0e:=8.0;
 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:=0;
 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:=X0c+Sqrt(sqr(R)-sqr(Y0c));
 Writeln('X2= ',X2:9:5,' - X координата второй точки пересечения');
 Writeln('Y2= ',Y2:9:5,' - Y координата второй точки пересечения');
 Step:=1;
 X3:=0;
 For i:=1 to 5 do
  begin
   Step:=Step/10;
   Repeat
    X3:=X3+Step;
    Y3:=Ysl(X3);
   Until Sqr((X3-X0e)/Ae)+Sqr((Y3-Y0e)/Be)>1;
   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+Rectangle4;
    Trap[k]:=Trapezium1+Trapezium2+Trapezium3+Trapezium4;
    Simp[k]:=Simpson1+Simpson2+Simpson3+Simpson4;
    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.