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)
@trusted pure nothrow @nogc
Unqual!(CommonType!(T1, T2))
poly
(
T1
T2
)
(
T1 x
,
in T2[] A
)
if (
isFloatingPoint!T1 &&
isFloatingPoint!T2
)
Unqual!(CommonType!(T1, T2)) poly(T1 x, T2[N] A)

Parameters

x T1: the value to evaluate.
A T2[]: 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