% square well potential. Width is 1, depth is 100.
function v1 = v1square(x)

inside = ( abs(x) < 1 );
v1(inside) = -100;
v1(!inside) = 0;
