Вывод функции interp1 в Matlab

Uses crt;
type massiv=array[1..20] of real;
const x:array[0..10] of real=(0,1,2,3,4,5,6,7,8,9,10);
      y:array[0..10] of real=(9.9,8.5,6.2,4.4,6.5,6.9,7.4,8.3,8.6,9.1,9.8);
var
a,b,c,d,h,k,l,r,s:massiv;
i,n:integer;
begin
Clrscr;
write('n='); read(n);
for i:=2 to n-1 do begin
k[1]:=0;l[1]:=0;
h[i-1]:=x[i-1]-x[i-2]; h[i]:=x[i]-x[i-1];
s[i]:=2*(h[i]+h[i-1]);
r[i]:=3*((y[i]-y[i-1])/h[i]-(y[i-1]-y[i-2])/h[i-1]);
k[i]:=(r[i]-h[i-1]*k[i-1])/(s[i]-h[i-1]*l[i-1]);
l[i]:=h[i]/(s[i]-h[i-1]*l[i-1]);
writeln('k[',i,']=',k[i]:2,' l[',i,']=',l[i]:2);
end;
for i:=n-2 downto 2 do begin
c[n-1]:=k[n-1];
c[i]:=k[i]-l[i]*c[i+1];
end;
for i:=1 to n-1 do begin
h[i]:=x[i]-x[i-1];
a[i]:=y[i-1];
b[i]:=(y[i]-y[i-1])/h[i]-h[i]*(2*c[i]+c[i+1])/3;
d[i]:=(c[i+1]-c[i])/3*h[i];
writeln('a[',i,']=',a[i]:2,' b[',i,']=',b[i]:2,' c[',i,']=',c[i]:2,' d[',i,']=',d[i]:2);
end;
readkey;
end.

x=[0.1 0.2 0.4 0.5 0.6 0.8 1.2];
y=[-3.5 -4.8 -2.1 0.2 0.9 2.3 3.7];
n=5;
k(1)=0;
l(1)=0;
for i=3:n-1
    h(i-1)=x(i-1)-x(i-2);
    h(i)=x(i)-x(i-1);
    s(i)=2*(h(i)+h(i-1));
    r(i)=3*((y(i)-y(i-1))/h(i)-(y(i-1)-y(i-2))/h(i-1));   
k(i)=(r(i)-h(i-1)*k(i-2))/(s(i)-h(i-1)*l(i-2));
l(i)=h(i)/(s(i)-h(i-1)*l(i-2));
end;
for i=n-2:3
    c(n-1)=k(n-1);
c(i)=k(i)-l(i)*c(i+1);
end;
for i=2:n-1
h(i)=x(i)-x(i-1);
a(i)=y(i-1);
b(i)=(y(i)-y(i-1))/h(i)-h(i)*(2*c(i-1)+c(i))/3;
d(i)=(c(i)-c(i-1))/3*h(i);
end;