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

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

hold off
plot(x, v1dkhead(x), 'b', x, v1dktail(x), 'k')
axis([-6,6,-200,50])
title('DK potentials','fontsize',20)
xlabel('{x} (L=1)','fontsize',18)
ylabel('V_1(x)  ( h^2/8\pi^2mL^2 = 1 )','fontsize',18)
legend('v1dkhead','v1dktail')
legend('location','southeast')