Вычисление интеграла;

[user100]

Помогите добавить в код вывод относительной и абсолютной погрешности.
И можно ли эту програмку реализовать без процедур и функций?

USES crt;
function f(x:real):real;
begin
f:=x*x+2*x+1;
end;
procedure rect (a,b:real; n:integer; var s:real);
var
i : integer;
h,x: real;
begin
h:=(b-a)/n;
x:=a+h/2;
s:=0.0;
for i:=1 to n do
begin
s:=s+f(x);
x:=x+h;
end;
s:=s*h;
end;
var a,b,e,s,pa,po:real;
    n:integer;
BEGIN
write('Начало интервала a=');
readln(a);
repeat
write('Конец интервала >',a:0:1,' b=');
readln(b);
until b>a;
repeat
write('Точность определения 0<e<1 e=');
readln(e);
until (e>0)and(e<1);
n:=round((b-a)/e);
rect (a,b,n, s);
write('Интеграл = ',s:4:4);
readln;
END.

[rubius2008]

Помогите добавить в код вывод относительной и абсолютной погрешности.

а что с чем сравнивается? аналитическое решение с численным?

И можно ли эту програмку реализовать без процедур и функций?

можно

[KobolD]

Что то ты с точностью напутал.
Обычно в таких задачах делается так:
интервал от а до b разбивается на n частей, считается интеграл (сумма), потом тот же самый интервал разбивается на 2*n частей и снова считают интеграл, а потом смотрят на отношение ( или разность) этих двух интегралов и уже его сравнивают с погрешностью, удовлетвоняет такая точность или нет. Если нет то увеличивают n и все по новой.