tofindzeros.m

衬底对量子点性质影响计算。。。。。。。。。。。。。。。

function
y=tofindzeros(E,n,m,V,d,a)
me=9.109e-31;
h=1.054e-34;
k=(2*E*me/h^2-(m*pi/d)^2).^0.5;
q=(2*V*me/h^2-k.^2).^0.5;
y=-besselj(n,k.*a)*(besselk(n-1,q.*a)-besselk(n+1,q.*a))/2.*q+...
    besselk(n,q.*a)*(besselj(n-1,k.*a)-besselj(n+1,k.*a))/2.*k;