From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- sys/vops/lz/apold.x | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100644 sys/vops/lz/apold.x (limited to 'sys/vops/lz/apold.x') diff --git a/sys/vops/lz/apold.x b/sys/vops/lz/apold.x new file mode 100644 index 00000000..885ed4fe --- /dev/null +++ b/sys/vops/lz/apold.x @@ -0,0 +1,25 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# APOL -- Evaluate a polynomial at X, given the coefficients of the polynomial +# in COEFF and returning the computed value as the function value. + +double procedure apold (x, coeff, ncoeff) + +double x # point at which the polynomial is to be evaluated +double coeff[ncoeff] # coefficients of the polynomial, lower orders first +int ncoeff + +int i +double pow, sum + +begin + sum = coeff[1] + pow = x + + do i = 2, ncoeff { + sum = sum + pow * coeff[i] + pow = pow * x + } + + return (sum) +end -- cgit