Программа вычисления интеграла методом трапеций

[shewlett]

Программа работает. Но преподаватель сказал, что надо преобразовать программу так, что бы она сама увеличивала (в 2 раза) изначально заданное число разбиений до тех пор, пока интеграл (площадь фигуры) подсчитанная при новом разбиении будет отличаться от интеграла при предыдущем разбиении станет меньше чем е=0.0001.
Может кто-нибудь помочь? (или хотя бы объяснить, что от меня хотят)

program rfr;
  var dx,x1,x2,e,i:real;
function Fx(x:real):real;
  begin
  Fx:=1.2*(900-x)*(exp(3*ln(cos*(0.05223*x))));
  end;
procedure CountViaTrap;
  var xx1,xx2,xx3:real;
  c:longint;
    begin
    writeln('Vsego interaciy:',round(abs(x2-x1)/e));
    i:=0;
    for c:=1 to round(abs(x2-x1)/e) do begin
    write('Interaciya',c,chr(13));
    xx1:=Fx(x1+c*e);
    xx2:=Fx(x1+c*e+e);
    if xx2>xx1 then xx3:=xx1 else xx3:=xx2;
    i:=i+abs(xx2-xx1)*e+abs(xx3)*e;
    end;
  writeln('Integral=', i);
  end;
      begin
      writeln('Vvedite ishodnie znacheniya:');
      write('Nachalnoe znachenie x(x1)=');readln(x1);
      write('Konechnoe znachenie x(x2)=');readln(x2);
      write('Tochnost vichisleniyya(e)=');readln(e);
      countViaTrap;
end.

[Аватар]

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

[Владимир_614]

Можете сообщить электронку на pniti2015@yandex.ru? Хотел бы ТЗ скинуть

[shewlett]

Уже сам разобрался. Всем спасибо)