1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
{ ******************************************************************
Polynomial evaluations for special functions.
Translated from C code in Cephes library (http://www.moshier.net)
****************************************************************** }
unit upolev;
interface
uses
utypes;
type
TabCoef = array[0..9] of Float;
function PolEvl(var X : Float; Coef : TabCoef; N : Integer) : Float;
function P1Evl(var X : Float; Coef : TabCoef; N : Integer) : Float;
implementation
function PolEvl(var X : Float; Coef : TabCoef; N : Integer) : Float;
{ ------------------------------------------------------------------
Evaluates polynomial of degree N:
2 N
y = C + C x + C x +...+ C x
0 1 2 N
Coefficients are stored in reverse order:
Coef[0] = C , ..., Coef[N] = C
N 0
The function P1Evl() assumes that Coef[N] = 1.0 and is
omitted from the array. Its calling arguments are
otherwise the same as PolEvl().
------------------------------------------------------------------ }
var
Ans : Float;
I : Integer;
begin
Ans := Coef[0];
for I := 1 to N do
Ans := Ans * X + Coef[I];
PolEvl := Ans;
end;
function P1Evl(var X : Float; Coef : TabCoef; N : Integer) : Float;
{ ------------------------------------------------------------------
Evaluate polynomial when coefficient of X is 1.0.
Otherwise same as PolEvl.
------------------------------------------------------------------ }
var
Ans : Float;
I : Integer;
begin
Ans := X + Coef[0];
for I := 1 to N - 1 do
Ans := Ans * X + Coef[I];
P1Evl := Ans;
end;
end.
|