найти оптимальное количество разбиений

[Надежда Ижбулатова]

сть программа поиска решения интеграла методо трапеции симпсона и прямокгольника.программа рабочая и считает правильно.нужно только прописать чтобы программа сама находила оптимальное количество разбиений n .вот как примерно это должно выглядеть только в одном месте не знаю что писать. препод сказал через див и мод.но у меня не получается.если можете помогите

program pr;
var
h,b,a:real; s1,s2,s3:integer;

e,i,n:integer;
function F(x:real):real;
begin
F:=exp(cos(x));
end;
begin
write('vvedite a=');
readln(a);
write('vvedite b=');
readln(b);
write('vvedite e=');
readln(e);
n:=1;
while (здесь не знаю сто вставить) do n:=n+1;

h:=(b-a)/n;
s1:=0;
s2:=0;
s3:=0;

for i:=1 to n do begin
s1:=f(a)*h+s1;
s2:=(f(a)+f(a+h))/2*h+s2;
s3:=(f(a)+4*f((2*a+h)/2)+f(a+h))*(h/6)+s3;
a:=a+h;
end;
writeln('s1=',s1:10:e);
writeln('s2=',s2:10:e);
writeln('s3=',s3:10:e);
writeln('n=',n);
readln;
end.

________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.

[BDA]

Определить n до подсчета интеграла нельзя.
Делается это так:
Считаем интеграл при n1 (в s1)
Увеличиваем n в 2 раза
Считаем интеграл при n2 (в s2)
Если разница между полученными интегралами больше eps, то повторяем подсчет (заносим s2 в s1, увеличиваем n и считаем интеграл в s2)