% --------------------------------------------
% Calculate H * psi for a given H=T+V and psi.
% --------------------------------------------
function  [y,a]=hpsi(psi,pote,kine)
% calculate h{psi}
a=fft(psi);
%a=fftshift(a);
a=a.*kine;
%a=fftshift(a);
b=ifft(a);
y=b+pote.*psi;
