Дифференциал функции

uses crt;
const n=10;//количество разбиений
var  h,h1,h2,s:real;
     a,b,x:real;
     dy1,dy2:real;
     i:integer;
     y:array[0..n] of real;
begin
writeln('Введите начальное и конечное значения a и b');
readln(a);
readln(b);
clrscr;
h:=(b-a)/n;//шаг разбиения промежутка
x:=(a+b)/2;//заданная точка
 s:=a;
 for i:=0 to n do
  begin
   y[i]:=1+1/s;
   s:=s+h;
  end;
 i:=trunc((x-a)/h+h/2);//округление до целой части с отбросом дробной
 h1:=2*h; h2:=h*h;
 dy1:=(-y[i-1]+y[i+1])/h1;//формула1-й производной
 dy2:=(y[i-1]-2*y[i]+y[i+1])/h2;//формула 2-й производной
 writeln('Значение a= ',a,' b= ',b);
 writeln('Искомая точка x= ',x);
 writeln('Шаг разбиения h= ',h);
  writeln('Значение первой производной dy1= ',dy1:5:5);
  writeln('Значение второй производной dy2= ',dy2:5:5);
end.