интеграл методом трапеций

[_SatanA_]

Доброго времени суток, уважаемые форумчане.
Столкнулся с проблемой.

Надо посчитать методом трапеций определенный интеграл с границами от 2 до 10

Вот код:

function f(x: Double): Double;
begin
  Result := Power(Ln(x), 3) + Sin(Exp(-x)) + x + Power(Tan(x*x*x),2)/Ln(x);
end;

procedure T(g: TChart; a, b: double; n: integer; var Sum: double);
var
 i: integer;
 h, x, x1, y, y_n: double;
begin
 Sum:=0;
 h:=(b-a)/n;
 for i:=0 to n-1 do
 begin
   x:=a+i*h;
   x1:=a+(i+1)*h;
   g.Series[0].AddXY(x, f(x));
   g.Series[0].AddXY(x1, f(x));
   Sum:=Sum+f(x);
 end;
 y:=f(a);
 y_n:=f(b);
 Sum:=(Sum+(y+y_n)/2)*h;
end;

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

Вот собсна интегральная функция

[denisbrain]

Сообщение от _SatanA_

Доброго времени суток, уважаемые форумчане.
Столкнулся с проблемой.

Надо посчитать методом трапеций определенный интеграл с границами от 2 до 10

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

Вот собсна интегральная функция

это к задачи о длине
береговой линии к фрактальной геометрий вопросы. чем меньше точка тем точнее длина. не хватает чисел для точного расчета Sin/Tan как вариант

[_SatanA_]

denisbrain, я в геометрии не силен, поэтому не совсем понял, что вы имели в виду )
проверил результат в маткаде, там при разном количестве разбиений получается примерно одинаковый результат