poly

Evaluate polynomial A(x) = a₀ + a₁x + a₂x² + a₃x³; ...

Uses Horner's rule A(x) = a₀ + x(a₁ + x(a₂ + x(a₃ + ...)))

Unqual!(CommonType!(T1, T2)) poly(T1 x, T2[] A)
Unqual!(CommonType!(T1, T2)) poly(T1 x, T2[N] A)
@safe pure nothrow @nogc
Unqual!(CommonType!(T1, T2))
poly
(
T1
T2
int N
)
(
T1 x
,
ref const T2[N] A
)
if (
isFloatingPoint!T1 &&
isFloatingPoint!T2
&&
N > 0
&&
N <= 10
)

Parameters

x T1: the value to evaluate.
A T2[N]: array of coefficients a₀, a₁, etc.

Examples

real x = 3.1L;
static real[] pp = [56.1L, 32.7L, 6];

assert(poly(x, pp) == (56.1L + (32.7L + 6.0L * x) * x));

Meta

Source

See Implementation

std math algebraic functions