Непонятен смысл не которых строк в коде - объяснить (дифференцирование функций)

type
diff=array[-1..10,0..2] of real;
vector=array [0..10] of real;
raznosty=array[0..5,0..5] of real;
 
function factorial(n: word): word; 
begin
{???}
if n=0 then Factorial := 1 else Factorial:= n*Factorial(n-1); 
end;
 
procedure Q(x:real; var razn:raznosty;  var y:real; var n:byte; var x0:vector);
var i,j:byte; t,ti:real;
begin
for j:=1 to n-1 do
for i:=0 to n-j-1 do
{???}
razn[i,j]:=razn[i+1, j-1]-razn[i,j-1];
t:=(x-x0[0])/(x0[1]+x0[0]);
ti:=t;
y:=razn[0,0];
for i:=1 to n-1 do
  begin
  {???}
  y:=y+(ti*razn[0,i])/factorial(i);
  ti:=ti*(t-i);
  end;
end;
 
var n,i,j:byte; h,y0,yN:real; deltaY:raznosty; y:diff; x:vector;
begin
writeln('Введите кол-во точек в которых известно значение функции');
readln(n);
for i:=0 to n-1 do
  begin
  writeln('Введите x',i,' и f(x',i,')');
  readln(x[i], y[i,0]);
  {???}
  deltaY[i,0]:=y[i,0];
  end;
{Шаг}
h:=x[1]-x[0];
{???}
Q(x[0]-h, DeltaY, y[-1,0],n,x);
Q(x[n-1]+h, DeltaY, y[n,0],n,x);
{Первые производная для крайних разностей}
y[0,1]:=(y[1,0]+y[0,0])/h;
y[n-1,1]:=(y[n-1,0]-y[n-2,0])/h;
{Первая производная для центральной разности}
for i:=1 to n-2 do
y[i,1]:=(y[i+1,0]-y[i-1,0])/(2*h);
{Вторая производная}
for i:=0 to n-1 do
y[i,2]:=(y[i-1,0]-2*y[i,0]+y[i-1,0])/(h*h);
for j:=1 to 2 do
  begin
  writeln('Значения производных ',j,'-го порядка');
  for i:=0 to n-1 do
  writeln('Для x=',x[i],'  y=',y[i,j]:4:2); 
  end;
end.

function factorial(n: word): Real; 
var
  i : word; r : Real;
begin
  r := 1;
  for i := 2 to n do r := r * i;
  factorial := r;
end;

ti:=ti*(t-i);

ti:=ti*(t-i)/(i+1);