Метод Монте-Карло(геметрически)

[4ainegg]

Помогите написать программу.Условие:подсчитать площадь фигуры,образованной осями координат и функцией y=a*e^(-b*x)+c,a,b,c- клавиатуры.После подсчитать погрешность в зависимости от кол-ва точек.За эталон использовать площадь,найденную интегралом.Проделать операцию(нахождение пл. по М-К 5 раз и вывести ответ в таблице для дальнейшей оценки точности метода)

вот то,что я смог сделать...

program poguch;
uses crt;
const a=-5;
b=-0.2;
c=10;
var h,i,n,nf,np,x,y:longint;
Sf,Sp,Si,p,u:real;
begin
clrscr;
randomize;
writeln('Metod Monte-Karlo');
writeln('y=a*e^(-bx) + c');
writeln('Primem a=-5,b=-0.2,c=10');
write('Vvedite kolichestvo tochek:');
readln(n);

nf:=0;
np:=0;
h:=0;
u:=0;
write(' np');
write(' nf');
write(' Pl pryamoygolnika');
write(' Pl figuri');
write(' Pl integr');
writeln(' Pogreshnost');
writeln;
for h:=1 to 5 do
begin
for i:=1 to n do
begin
x:=random(50);
y:=random(50);
if (x<=-ln(-c/a)/b) and (y<=(a+c)) then
inc(np);
if (x<=-ln(-c/a)/b) and (y<=a*exp(-b*x)+c) then
inc(nf);
end;

Sp:=(a+c)*(-ln(-c/a)/b);
Sf:=(nf/np)*Sp;
Si:=-25*exp(ln(2))+10*5*ln(2)+25;
p:=abs((Sf-Si)/Si)*100;
u:=u+p;
write(' ',np:3);
write(' ',nf:3);
write(' ',Sp:2:4);
write(' ',Sf:3:2);
write(' ',Si:3:3);
writeln(' ',p:3:3,'%');
readln;
end;
u:=u/5;
writeln('sredna pogreshnost=',u:7:5,'%');
readln
end.

np - кол-во точек,попавших в прямоугольник,описанный вокруг искомой фигуры.
nf - кол-во точек,попавших в искомую фигуру.