Численное решение ДУ методом Эйлера + правило рунге - Turbo Pascal

Формулировка задачи:

попытался написать, но работает некорректно может кто подсказать в чем ошибка?

program MethosEyl; //Метод Эйлера для ДУ
uses
  Classes;
const a=0; b=1;E=0.0001;
          function F(x,y:real):real;
       begin
         F:=2*x-sqr(y);
       end;
var x,y,h,s1,s2:real; n,i:integer;
begin
  y:=0.1;
  n:=10;
  s1:=0;
  s2:=0.1;
  repeat
    s1:=s2;
 h:=(b-a)/n;
 x:=a;
 s2:=0.1;
 for i:=0 to n do begin
   x:=a+i*h;
   write('x[',i,']=',x:0:4);
write('   ');
writeln('y[',i,']=',y:0:4);
y:=y+h*F(x,y);
s2:=s2+y;
 end;
n:=2*n;
  until ((1/15)*abs(s2-s1))<=E;
readln;
end.

Решение задачи: «Численное решение ДУ методом Эйлера + правило рунге»

Листинг программы

until ((1/15)*abs(s2-s1))>=E;

Объяснение кода листинга программы

Переменная s1 и s2 - это значения, которые сравниваются в данном выражении.
(1/15)*abs(s2-s1) - это вычисление разности между s2 и s1, деленное на 15, а затем абсолютированное. Результат этого выражения будет числом от 0 до 1.
(1/15)*abs(s2-s1)>=E - это условие, которое проверяет, что результат сравнения больше или равен E.
while - это ключевое слово, которое обозначает цикл или условие продолжения.
do - это ключевое слово, которое обозначает блок кода, который будет выполняться до тех пор, пока условие истинно.
until - это ключевое слово, которое обозначает условие окончания цикла.
((1/15)*abs(s2-s1))>=E; - это полное выражение, которое будет выполнено в цикле до тех пор, пока условие истинно.