clc;clear;
a=50*1e-10;%m
me=9.1*1e-31;%kg
V0=1;%eV
V_J=1.602e-19;
h_bar=1.0545*1e-34;%J*s
k0=sqrt（2*me*V0*V_J/（h_bar^2））;
z0=k0*a/2;
E=-1:0.001:0;
k=sqrt（2*me.*（E+V0）*V_J/h_bar^2）;
z=k*a/2;
y1=tan（z）;
y2=sqrt（（z0./z）.^2-1）;
figure（1）
plot（Ey1‘r--‘Ey2‘g.-‘）
legend（‘tan（z）‘‘sqrt（（z0./z）.^2-1）‘）
grid?on
xlabel（‘E（eV）‘）
axis（[-1?0?-18?18]）
y=y1-y2;
[n1m1]=find（abs（y（1:100））<1）;
E1=E（m1）;
[n2m2]=find（abs（y（100:200））<0.1）;
E2=E（100:200）;
E2=E2（m2）;
[n2m3]=find（abs（y（300:400））<0.01）;
E3=E（300:400）;
E3=E3（m3）;
EE=[E1E2E3];
figure（2）
plot（[-50?-25?-25?25?25?50][0?0?-1?-1?0?0]）;hold?on
for?i=1:3
plot（[-25?25][EE（i）?EE（i）]）;hold?on
end
axis（[-50?50?-1.2?0.2]）
text（-10-0.05‘阱寬50A‘）;%文本標(biāo)注
text（-15-1.1‘阱深1eV‘）;%文本標(biāo)注
text（-15E1‘基態(tài)-0.9870eV‘）;%文本標(biāo)

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1656??2015-04-08?14:04??Untitled.m