function h = hermite(degree, x)


%% Function returns (0-6 degree) hermite polynomial, evaluated at x. Inputs:
%% 
%% degree: an int from 0 to 6. Degree of the polynomial.
%% x: value(s) where polynomial is to be evaluated. May be scalar or
%%    more general vector, matrix, or tensor. Entrywise calculations will
%%    be made, return value will have same dimensions.
%% 
%% Ken Gosier
%% May 2002


if 0 == degree
   h = 1;
elseif 1 == degree
   h = -x;
elseif 2 == degree
   h = x.^2 - 1;
elseif 3 == degree
   h = -x.^3 + 3*x;
elseif 4 == degree
   h = x.^4 - 6*x.^2 + 3;
elseif 5 == degree
   h = -x.^5 + 10*x.^3 - 15*x;
elseif 6 == degree
   h = x.^6 - 15*x.^4 + 45*x.^2 - 15;
else
   T = sprintf('bad degree val %d must be an int from 0 to 6', degree);
   error(T);
end