Kronig Band Structure Matlab

.602e-19)/hbar * (a-b)) .* ... cosh(sqrt(2*m*(V0 - E)*1.602e-19)/hbar * b) + ... ((V0 - 2*E)/ (2*sqrt(E*(V0 - E)))) .* ... sin(sqrt(2*m*E*1.602e-19)/hbar * (a-b)) .* ... sinh(sqrt(2*m*(V0 - E)*1.602e-19)/hbar * b); % Find energies E where f(E) = 0 (transcendental equ