Метод Симпсона на чётном отрезке разбиения (Delphi)

[Atherakhia]

Здравствуйте, написал программу для нахождения определённого интервала методом Симпсона, но она работает неправильно. Прошу помочь в нахождении ошибки.

program project1;

var
  x0, yn, y0, a, b, h, s1, s2, int: real;
  i, n: byte;

function F(x: real): real;
begin
  F:=2*x;
end;

begin
write('Vvedite nachalnoe znacenie funkcii=');
readln(a);
write('Vvedite konecnoe znacenie funkcii=');
readln(b);
write('Vvedite tchislo otrezkov razbienia=');
readln(n);
if (n mod 2)<>0 then begin
  n:=n+1;
end;
h:=(b-a)/n;
x0:=a;
y0:=F(x0);
for i:=1 to (n-1) do begin
  if (i mod 2)<>0 then
  begin
    y0:=F(x0);
    x0:=x0+h;
    s1:=y0;
  end;
end;
for i:=2 to (n-2) do begin
  if (i mod 2)=0 then
  begin
    y0:=F(x0);
    x0:=x0+h;
    s2:=y0;
  end;
end;
for i:=1 to n do begin
  y0:=F(x0);
  x0:=x0+h;
  yn:=y0;
end;
int:=(h/3)*(y0+(4*s1)+(2*s2)+yn);
writeln('int=', int:6:3);
readln
end.

[newerow1989]

program project1;

var
  x0, yn, y0, a, b, h, s1, s2, int: real;
  i, n: byte;

function F(x: real): real;
begin
  F:=2*x;
end;

begin
write('Vvedite nachalnoe znacenie funkcii=');
readln(a);
write('Vvedite konecnoe znacenie funkcii=');
readln(b);
write('Vvedite tchislo otrezkov razbienia=');
readln(n);
if (n mod 2)<>0 then begin
  n:=n+1;
end;
h:=(b-a)/n;
x0:=a;
y0:=F(x0);
for i:=1 to (n-1) do begin
  if (i mod 2)<>0 then
  begin
    y0:=F(x0);
    x0:=x0+h;
    s1:=s1+y0;
  end;
end;
x0:=a;
for i:=2 to (n-2) do begin
  if (i mod 2)=0 then
  begin
    y0:=F(x0);
    x0:=x0+h;
    s2:=s2+y0;
  end;
end;
for i:=1 to n do begin
  y0:=F(x0);
  x0:=x0+h;
  yn:=yn+y0;
end;
int:=(h/3)*(y0+(4*s1)+(2*s2)+yn);
writeln('int=', int:6:3);
readln
end.

Проверяй, я точно не знаю в чем состоит метод Симпсона!