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

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

  n:=10;

  for i:=1 to n do
  begin
    x[i]:=i*0.1875;
    y[i]:=sin(x[i]);
  end;

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

  { ------ Berechnung der Ableitung ------ }
  interpol_val:=spline_deriv(pi/8.0,x,y,z,n);

  writeln;
  writeln('Berechnung des interpolierten Werts der Ableitung aus den Stuetzstellen');
  writeln;
  writeln('Der Wert der Ableitung der Sinusfkt. bei pi/8 lautet :',interpol_val);
  writeln('Der exakte Wert lautet : ',cos(pi/8.0));
end.