aboutsummaryrefslogtreecommitdiff
path: root/sys/vops/apol.gx
blob: 04d162c5843bd21e3fbcc0df7036987fb9f5fed1 (plain) (blame)
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
# 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.

PIXEL procedure apol$t (x, coeff, ncoeff)

PIXEL	x		# point at which the polynomial is to be evaluated
PIXEL	coeff[ncoeff]	# coefficients of the polynomial, lower orders first
int	ncoeff

int	i
PIXEL	pow, sum

begin
	sum = coeff[1]
	pow = x

	do i = 2, ncoeff {
	    sum = sum + pow * coeff[i]
	    pow = pow * x
	}

	return (sum)
end