function y = psolve( p, x )
% PSOLVE Polynomial solver
% Solves a polynomial (given a coeff. vector) for a vector input
% J. Watlington, 11/18/95

[m,num_samples] = size( x );
[m,poly_order] = size( p );

y = zeros( [ 1, num_samples ]);     %  preallocate
x_power = ones( [ 1, poly_order ]);
for loc = 1:num_samples
	for power = 2:poly_order
		x_power( power ) = x_power( power - 1 ) * x( loc );
	end
	y( loc ) = p * x_power';
end