clear all;
set(0,'DefaultFigureWindowStyle','docked');

mue = 1;
V = @(x) cos(x); % potential
VH = @(x) - x.^2 / (2 * mue); % hermite correction
N = 250;
Nx = 500;

[x w] = hermpts(Nx);
B = hermiteS(0:N-1,x);

serr = @(A, B) fprintf('absolute error: %g\n', ...
                       max(abs(A(:)-B(:))));