Вычисление интеграла методом трапеции.

[Айдар]

Добрый день помогите с данной задачи. надо вычислить интеграл c точностью 0,001 (любой самый простейший можно) вот набрасал код:

var a,b,s,x,e,n:real;
    i:integer;
begin
readln(a,b,n);
h:=(b-a)/n;
s:=0;
e:=0,001;
for i:=1 to n do
 begin
{    
вот эту часть не могу записать, т.к. с инегралами только начал знакомиться, тут еще как я понял надо будит увеличивать цикл (N)
в 2 раза, и сравнивать разность s и sk 
}
 end;
writeln(s);
end.

Вот формула нахождения интеграла
s=[1/2*(F(a)+F(b)+F(a+h)+F(a+2h)+...+F(a +(N-1)*h))]*h

[L_M]

А обязательно ли пользоваться этой формулой? Интеграл - это площадь. То есть если разбить на h равных кусков длинной x каждый, то площади будет просто x*(f(x1)+f(x2)+...). А точность увеличивается с увеличением количества кусков, на которые разбивается. И вообще, как предполагается считать точность? Значит должно быть изначально дано точное решение?

[Айдар]

нет находим 1 решение интеграла, далее делим трапецию на 2 раза меньше (n/2),считаем, если их разность равна 0.001 или меньше, то выводим ответ и кол-во отрезков (n)

Вот написал, все правильно считает, но надо сделать так чтобы он мне выводил ответ при N повторов и 2*N повторов, у меня почемуто все затеряется и выводит ответ только при 2*n повторов, что делать?

var
n,i:integer;
a,b,e,shag,sum,itog,itog2:real;
function F(x:real):real;
begin
f:=x*x+3;
end;
begin
write ('a=');
readln (a);
write ('b=');
readln (b);
write ('N=');
readln (n);
e:=0.001;
writeln ('e=0.001');
{shag:=(b-a)/n;}
sum:=0;
repeat
shag:=(b-a)/n;
sum:=0;
for i:=1 to n-1 do
begin
sum:=sum+f(a+i*shag);
end;
sum:=sum+(f(a)+f(b))/2;
itog:=(b-a)/n*sum;
itog2:=itog;
n:=2*n;
writeln (itog:10:5);
writeln (itog2:10:5)
until itog2-itog<e;
{writeln ('Otvet=',itog:0:5);
writeln ('Otvet2=',itog2:0:5);}
writeln (n);
readln;
end.