{$I DECL.TLD}
{$I MAT2.TLD}
{$I SPLINES.TLD}

var x,y,xs,y_int : RealVectorType_TL;
    n_int,n,i    : integer;
    interpol_val : real;
begin
  clrscr;
  writeln('Spline-Interpolation eines Felds');
  writeln;
  n:=10;

  { ----- Erzeugung der Daten ----- }
  for i:=1 to n do
  begin
    x[i]:=i;
    y[i]:=x[i]*x[i];
  end;

  n_int:=15;

  { ----- Definition der Stuetzstellen ----- }
  for i:=1 to n_int do xs[i]:=2.25+(i-1)*0.5;

  { ----- Interpolation ----- }
  spline_interpol_arr(x,y,n,xs,y_int,n_int);

  writeln('        x                  interpoliert                    exakt');
  writeln;
  for i:=1 to n_int do
     writeln(xs[i],'      ',y_int[i],'         ',xs[i]*xs[i]);
end.