write('d3y ');
for i:=1 to kol-3 do
d3y[i]:=d2y[i+1]-d2y[i];
for i:=1 to kol-3 do
write(d3y[i]
3,' ');
writeln;

write('d4y ');
for i:=1 to kol-4 do
d4y[i]:=d3y[i+1]-d3y[i];
for i:=1 to kol-4 do
write(d4y[i]
3,' ');

вот кусок кода

мне нужно, чтобы цикл автоматом создавал вот эти "dXy"

и соответсвующие к ним операции
точнее выполнял их

т.е. сейчас код заканчивается на d5y, но этого мало, нужно до 128, однако руками же не будешь копипастить 128 циклов))

нужно как-то автоматизировать)

d3y, d4y и все последующие - это массивы
типа real

Нет ничего проще!
Создай 2мерный статический массив (лучше, конечно, динамический, но...)

const kol = 128;
var dXy: array[kol - 1][kol] of real;

//...
for j := 2 to Kol - 1 do
begin
  write('d', i:3, 'y :');
  for i:=1 to kol-j do
  begin
    dxy[j][i]:=dxy[j-1][i+1]-dxy[j-1][i];
    write(dxy[j][i], ' ');
  end;
  writeln;
end;

Как-то так, не проверял.