% bundna.m  
% övn. uppg. 4.6 i 'Kvantvärldens fenomen'

close all
clear all

hbarC=197.3271; % [e.V.nm]
b=1/2; % [nm] OBS! halva brunnens bredd.
mC2=510999; % [e.V.]

konst=sqrt(2*mC2*b*b/(hbarC*hbarC))

V_0=5; % [e.V.]         

E=linspace(-V_0+5*eps,-5*eps,100);

f_jamn= tan(konst*sqrt(E+V_0) ) ; % villkor for jamn paritet, sidan 74 i 'Kvantvarldens fenomen'.
g_jamn= sqrt(-E./(E+V_0) ); 

subplot(2,1,1)
plot(E,f_jamn)
hold on
plot(E,g_jamn,'r')
xlabel('E [e.V.]','fontsize',15)
title('Jämn paritet','fontsize',15)
axis([-5.5 0.5 -10 10])

f_udda= tan(konst*sqrt(E+V_0) ) ; % villkor for udda paritet, sidan 75 i 'Kvantvarldens fenomen'.
g_udda= -sqrt( (E+V_0)./(-E) ) ;

subplot(2,1,2)
plot(E,f_udda)
hold on
plot(E,g_udda,'r')
xlabel('E [e.V.]','fontsize',15)
title('Udda paritet','fontsize',15)
axis([-5.5 0.5 -10 10])




