...Wert eines Polynoms n-ten Grades an einem gegebenen Punkt bestimmen?

// Simultaneous evaluation of a given polynomial and its first derivative at a given point
// Simultane Berechnung des Wertes eines Polynoms n-ten Grades und seiner Ableitung

type
  TPolynomArray = array of Double;

procedure Horner(Polynom: TPolynomArray; X: Extended; var FX, derivation: Extended);
var
  i: Integer;
  H: Integer;
begin
  H := High(Polynom);
  FX := Polynom[H];
  derivation := 0.0;
  for i := H - 1 downto 0 do
  begin
    derivation := derivation * X + FX;
    FX := FX * X + Polynom[i];
  end;
end;

{Beispiel / Sample code }

procedure TForm1.Button1Click(Sender: TObject);
var
  X, FX, derivation: Extended;
begin
 (* Horner's scheme give an algorithm for the simultaneous evaluation
   of a given polynomial and its first derivative at a given point *)

 (* Hornerschema zur Berechnung eines Polynoms n-ten Grades an einem
    bestimmten Punkt *)

  (* f(x) = 3 x^5 + 4 x^4 + 13 x^3 - 59 x^2 + 19 x - 97 *)

  X := 2.5;
  Horner(VarArrayOf([-97, 19, - 59, 13, 4, 3]), X, FX, derivation);
  ShowMessage(Format('x = %n'#13#10'f(x) = %n'#13#10'f''(x) = %n', [X, FX, derivation]));
end;

printed from
www.swissdelphicenter.ch
developers knowledge base