Ну что ты мучаешься, пройдись отладчиком и посмотри где, что не так. У тебя проблемы в выходе за пределы диапазона в массивах.
Код:
for I := 1 to (n-j) do
y2[i] := (((y[i + 2] - y[i + 1]) / (x[i + j] - x[i])) - (y[i + 1] - y[i]) / (x[i + j] - x[i])) / (x[i + 2] - x[i]);
Вот здесь у тебя выход x[i + 2], здесь y2[i]