Изменение программы (многочлен Лагранжа)

uses crt;
type
  massiv=array [0..3] of real;
  mas = array [0..2] of real;
var
  x,y,a:massiv;
  a1:mas;
  i,j:byte;
  z:real;
 
procedure vvid;
  begin
    for i:=0 to 3 do
    begin
      write('X',i,' = '); readln(x[i]);
      write('Y',i,' = '); readln(y[i]);
    end;
  end;
 
procedure obch;
var ost:real;
  begin
    for i:=0 to 3 do
    begin
      z:=1;
      a1[0]:=1;
      a1[2]:=0;
      for j:=0 to 3 do
      begin
        if j<>i then
        begin
          z:=z*(x[i]-x[j]);
          a1[0]:=a1[0]*(-x[j]);
          a1[2]:=a[2]-x[j];
          ost:=-x[j];
        end;
      end;
      a1[1]:=a1[0]/ost+(a1[2]-ost)*ost;
      a[3]:=a[3]+y[i]/z;
      a[2]:=a[2]+a1[2]*y[i]/z;
      a[1]:=a[1]+a1[1]*y[i]/z;
      a[0]:=a[0]+a1[0]*y[i]/z;
    end;
  end;
 
procedure vuvid;
  begin
    writeln('y = ',a[3]:3:4,'*x^3 + (',a[2]:3:4,'*x^2) + (',a[1]:3:4,'*x) + (',a[0]:3:4,')');
  end;
 
BEGIN
clrscr;
vvid;
obch;
vuvid;
readkey;
END.