The general solution in the region (n−1)a < x < nais
ψn(x) =Ansin(k[x−na]) +Bncos(k[x−na])k=sqrt2 mE
̄h^2Now lets look at the boundary conditions atx=na. Continuity of the wave function gives
ψn(na) = ψn+1(na)
Ansin(0) +Bncos(0) = An+1sin(−ka) +Bn+1cos(−ka)
Bn = −An+1sin(ka) +Bn+1cos(ka)Bn+1 =Bn+An+1sin(ka)
cos(ka).
The discontinuity in the first derivative is
dψn+1
dx∣
∣
∣
∣
na−
dψn
dx∣
∣
∣
∣
na=
2 maV 0
̄h^2ψn(na)k[An+1cos(−ka)−Bn+1sin(−ka)]−k[Ancos(0)−Bnsin(0)] =
2 maV 0
̄h^2Bnk[An+1cos(ka) +Bn+1sin(ka)−An] =2 maV 0
̄h^2BnSubstitutingBn+1from the first equation
k[An+1cos(ka) + [Bn+An+1sin(ka)] tan(ka)−An] =2 maV 0
̄h^2BnAn+1(cos(ka) + sin(ka) tan(ka)) +Bntan(ka)−An=
2 maV 0
̄h^2 kBncos^2 (ka) + sin^2 (ka)
cos(ka)An+1=2 maV 0
̄h^2 kBn−Bntan(ka) +AnAn+1=2 maV 0
̄h^2 kBncos(ka)−Bnsin(ka) +Ancos(ka)Plugging this equation forAn+1back into the equation above forBn+1we get
Bn+1 =Bn+An+1sin(ka)
cos(ka)Bn+1 =Bn+( 2 maV 0
̄h^2 k Bncos(ka)−Bnsin(ka) +Ancos(ka))
sin(ka)
cos(ka)Bn+1 =Bn
cos(ka)+
(
2 maV 0
̄h^2 kBnsin(ka)−Bnsin^2 (ka)
cos(ka)+Ansin(ka))
Bn+1 =Bn
cos(ka)+
(
2 maV 0
̄h^2 kBnsin(ka)−Bn(
1
cos(ka)−cos(ka))
+Ansin(ka))
Bn+1 =2 maV 0
̄h^2 kBnsin(ka) +Bncos(ka) +Ansin(ka).