uses MtxExpr,
Math387,
MtxVec, MtxVecEdit,
MtxVecTee,
Polynoms;
procedure TForm1.Button1Click(Sender: TObject);
var X,Y,Y2: Vector;
PP :
TPiecePoly;
i : Integer;
YVal : double;
begin
PP := TPiecePoly.Create;
try
// generate function - note that X values are monotonical
X := Ramp(100,0,1);
Y := RandUniform(100,0,50) + Ramp(100,100,0.25);
// construct cubic splines, but do not evaluate them
Spline1D(X,Y,PP);
X := Ramp(800,0,0.125);
//get interpolation points
PP.Evaluate(X,Y2);
// evaluate
DrawIt(Y,'Original');
DrawIt(Y2,'Interpolated');
finally
PP.Free;
end;
end;