% Plot several of the potentials for handy comparison.
function v1potentialsTest

x = linspace(-4, 4, 1000);

plot(x, v1square(x), 'k',  x, v1rounded(x), 'b',...
     x, v1skewed(x), 'g')
axis([-4,4,-120,20])
title('Potential Wells','fontsize',20)
xlabel('{x} (L=1)','fontsize',18)
ylabel('V_1(x)  ( h^2/8\pi^2mL^2 = 1 )','fontsize',18)
legend('v1square','v1rounded','v1skewed')
legend('location','southeast')