Вычисление площади фигуры

var
a,b,c,d,e,Sold,Snew,a1,b1 : double;
p1,p2 : boolean;
n,k,m :integer;
procedure vvod();
begin
write('Задайте погрешность вычислений: ');
readln(e);
write('Задайте длину первой полуоси эллипса: ');
readln(a);
write('Задайте длину второй полуоси эллипса: ');
readln(b);
write('Задайте первую полудиагональ ромба: ');
readln(c);
write('Задайте вторую полудиагональ ромба: ');
readln(d);
end;
 
function check(x,y :double):boolean;
begin
p1:=((sqr(x)/sqr(a)+ sqr(y)/sqr(b))<=1)and (y>=0);
p2:=(abs(x)<=c) and(abs(y)<=-d/c+d);
if p1 and p2 then check:=true else check:=false;
end;
 
procedure MonteCarlo_m(a1,b1 :double; k: longint; var m: longint);
var
i: integer;
x,y :double;
flag :boolean;
begin
m:=0;
n:= sqr(k);
for i:=1 to n do begin
x:=random*(2*a1+1)- a1;
y:=random*(2*b1+1)- b1;
flag:=check(x,y);
if flag = true then m:=m+1;
end;
end;
 
begin
randomize;
k:=1;
Snew:=10000;
vvod();
repeat
Sold:=Snew;
if (abs(a)>abs(c)) then a1:=a else a1:=c;
if (abs(b)>abs(d)) then b1:=b else b1:=d;
MonteCarlo_m(a1,b1,k,m);
Snew:=m/n*(a1*b1);
k:=k+1;
until (abs(Sold-Snew)<e);
writeln('Площадь: ', Snew);
writeln('Нажмите любую клавишу...');
read();
exit;
end.