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

x = linspace(-3.5, 4.5, 1000);

plot(x, v1doubledip(x), 'g', x,   v1localization(x), 'b',...
   [-3,3], [-10.25,-10.25], 'r-o')
% 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('v1doubledip','v1localization','V = -10.25')
legend('location','east')