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

var x,y,z        : RealVectorType_TL;
    n,i          : integer;
    interpol_val : real;
begin
  clrscr;
  writeln('Spline-Interpolation');
  writeln;
  n:=10;

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

  { ----- Spline-Koeffizienten ----- }
  spline_coeff(x,y,z,n);

  interpol_val:=spline_interpol(pi/4.0,x,y,z,n);

  writeln;
  writeln('Der interpolierte Wert der Sinusfkt. bei x=pi/4 lautet :',interpol_val);
end.