Помогите пожалуйста исправить недоработку в коде программы.

program Hermite;
var i,j,s,n,k,i1,m,j1,i2,j2: integer;
b: array[1..100,1..100] of real;
x: array[1..100] of real;
var a: array[1..100] of integer;
y: array[1..100] of real;
z: array[1..100] of real;
f,ml,kl,l,h: real;
var
zl: array[1..16,1..2] of real;
begin
writeln('vvedite chislo dannyh tochek');
readln(m);
for i:=1 to m do
begin
write('vvedite kratnost',i,'tochki i=');
readln(a[i]);
end;
n:=0;
for i:=1 to m do begin
n:=n+a[i];
end;
k:=1;
{задаём значения x узловых точек, вводим функцию, вычисляем значения функции в искомых точках}
for i:=1 to m do begin
write('x[',i,']=');
readln(x[i]);
y[i]:=1/25*cos(5*x[i]);
for j:=k to k+a[i]-1 do begin
b[j,1]:=x[i];
b[j,2]:=y[i];
end;
k:=k+a[i];
end;
{задаём значения производных в узловых точках, заносим эти данные в таблицу}
k:=1;
for i:=1 to m do begin
for j:=3 to a[i]+1 do begin
write(j-2,'proizv',i,'tochki=');
readln(z[j-2]);
for j1:=k to k+a[i]-j+1 do begin
b[j1,j]:=z[j-2];
end;
end;
k:=k+a[i];
end;
{Вычисляем остальные разделённые разности, не заданные изначально и заносим их в таблицу}
for j:=3 to n+1 do begin
for i:=1 to n-j+2 do begin
if b[i,1]<>b[i+j-2,1] then
b[i,j]:=(b[i,j-1]-b[i+1,j-1])/(b[i,1]-b[i+j-2,1]);
end;
end;
[COLOR="RoyalBlue"]{Задаём точки, в которых будет вычисляться значения интерполяционного многочлена Эрмита}
l:=b[1,1];
for i2:=1 to 13 do begin
zl[i2,1]:=1;
h:=0.25;
l:=1+h;[/COLOR]
end;
{используя формулы интерполяции, вычисляем значения интерполяционного многочлена Эрмита в заданных точках функции}
for i2:=1 to 13 do begin
f:=b[1,2];
for j:=3 to n+1 do
begin
ml:=1;
i:=1;
while i<=j-2 do begin
ml:=ml*(zl[i2,1]-b[i,1]);
i:=i+1;
end;
f:=f+b[1,j]*ml;
end;
zl[i2,2]:=f;
end;
{вывод таблицы разделенных разностей}
writeln('vyvod tabl radelennih raznostei');
for i1:=1 to n do begin
for j1:=1 to n+1 do begin write(b[i1,j1]:3:2,'');
end;
writeln;
end;
writeln('-----------------');
{вывод матрицы, первый столбец которой содержит значения точек, а второй значения интерполяционного многочлена Эрмитчках}
writeln('znacheniya pribligaushei funkcii');
writeln('x Hn(x)');
for i2:=1 to 13 do begin
for j2:=1 to 2 do begin
write(zl[i2,j2]:3:4,'');
end;
writeln;
end;
readln;
end.