Documentation of hermite

Global Index (all files) (short | long) | Local Index (files in subdir) (short | long)

Function Synopsis

y = hermite(x, n);

Help text

  y = hermite(x, n);

Listing of function hermite

function y = hermite(x, n);

if size(x, 1) == 1;
  x = x'; rot = 1;
end;

ytem = zeros(length(x), 1);
ytem(:,1) = ones(length(x), 1);
ytem(:,2) = 2*x;

if n == 0; y = ytem(:,1); return; end;
if n == 1; y = ytem(:,2); return; end;

if n > 1;

  m = 1;
  p = 2;

  for i = 2:n;
    ynew = 2*x.*ytem(:,p)-2*(i-1)*ytem(:,m);
    tem = m;
    m = p;
    p = tem; 
    ytem(:,p) = ynew;
  end

  y = ynew;

end;