% Skewed potential well. FWHM is 1, depth is 100.
% All of the nonzero potential is contained within -2.5 < x < 2.5,
% and the integral of the potential equals -2 (same as v1square).
function v1 = v1skewed(x)


steep = (x >= -1.5) & (x <= -0.5);
v1( steep ) = 50 * ( - 1 - sin( (1+x(steep)) * pi) );

gradual = (x >= -0.5) & (x <= 2.5);
v1( gradual ) = 50 * ( - 1 - sin( (1-x(gradual)) * pi/3) );

v1( (x <= -1.5) | (x >= 2.5) ) = 0;
