Procedure FPOLY(X:real; var P:array of real; NP:integer);
var
    J:integer;
begin
    P[1]:=1;
    For J:=2 To NP do
      P[J]:=P[J - 1] * X;
end;