задача по методам прямоугольника и трапеции

[Flur]

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

program integrirov;
uses crt,GRAPH;
const
a=1;
b=4;
var
dr,md:integer;
hs,hsp,hspp:real;
n:integer;
integral:real;
pramoyg:real;
trap:real;
sch,snch:real;
x:real;
i:integer;
x0,y0,mx,my,px,py,px2,py2,pysr,pxsr:integer;
y,y2,ysr,dx,m,ap,bp:real;
s:string;
          function Fx(g:real):real;
          begin
          Fx:=((ln(x)*ln(x))/x);
          end;
begin
clrscr;
dr:=detect;
writeln('Metod pramoygolnikov');
writeln('Vvedite chislo razbienii');
readln(n);
Hsp:=(b-a)/n;
hspp:=a+Hsp/2;
x:=a;
for i:=1 to n-1 do begin
x:=hspp+i*hsp;
pramoyg:=pramoyg+Fx(x)*hsp;
end;
writeln('Chislennoe znachenie po metody pramoygolnikov=',pramoyg);
writeln('Chislennoe znachenie po metody trapessi=');
writeln('Chislo rozbienii=',n);
Hs:=(b-a)/n;
sch:=0;
for i:=1 to n-1 do
sch:=sch+fx(a+hs*i);
trap:=hs*((fx(a)+fx(b))/2+s);
writeln('Chislennoe znachenie po metody Trapesii=',trap);
writeln('Nacmite enter');
readln;
initGraph(dr,md,'');
mx:=100; {masshtab po X}
my:=100; {masshtab po Y}
x0:=10;
y0:=getmaxY-2;
line(x0,0,x0,y0+2);
line(0,y0,getmaxX-1,y0);
for i:=1 to 10 do begin
line(x0+i*mx,y0-3,x0+i*mx,y0+3); {zasechki po osi OX}
str(i,s);
outtextXY(x0+i*mx,y0+10,s);
line(x0-3,y0-(i*my),x0+3,y0-(i*my)); {zasechki po osi OY}
str(i,s);
outtextXY(x0+10,y0-i*my,s);
end;
{grafic}
x:=a;
dx:=0.001;
while x<=b do begin
x:=x+dx;
px:=x0+round(x*mx);
y:=Fx(x);
py:=y0-round(y*my);
putpixel(px,py,10);
end;
{razbienie na pramoygolniki}
x:=a;
while x<=b do begin
px:=x0+round(x*mx);
y:=Fx(x+hsp/2);
py:=y0-round(y*my);
if x<b then begin
px2:=x0+round((x+hsp)*mx);
pxsr:=x0+round((x+hsp/2)*mx);
ysr:=Fx(x+hsp/2);
y2:=Fx(x+hsp);
pysr:=y0-round(ysr*my);
py2:=y0-round(y2*my);
end;
if (x=a) or (y=b) then
line(px,y0,px,py) else line(px,y0,px,pysr);
line(px,pysr,px2,pysr);
x:=x+hsp;
end;

{closeGraph;}
readln;
end.

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

[rubius2008]

Ответы по идеи должны совпадать с некоторой погрешностью. Возьмите линейную функцию a*x+b, ответы по идеи должны быть близки. Еще можно проверить, уменьшая шаг (например, в 10 раз), если ответы начинают совпадать значит - хорошо.

[Flur]

СПАСИБО всем кто думал над задачей!!! сегодня сдал зачет.

[Stilet]

Я не думал. У меня с математикой отношения натянутые.