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

var   x,y,z,y_int  : RealVectorType_TL;
      n,i          : integer;
      interpol_val : real;
begin
  clrscr;
  writeln('Numerische Integration 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 Koeffizienten ------ }
  spline_coeff(x,y,z,n);

  { ----- Definition der Stuetzstellen ----- }
  for i:=1 to n do y_int[i]:=(i-1)*0.09375;

  { ----- Integration ----- }
  spline_integrat_arr(x,y,y_int,n);

  writeln;
  gotoxy(3,5);writeln('Splineintegration      exakt');
  writeln;

  for i:=1 to n do writeln(y_int[i],cos(x[1])-cos(x[i]));
end.