clear
basis=[
2.46?0?0
1.23?2.13?0
0?0?10?
?????]

?b=[?
0?0?0
0?1.42?0
??];

all=[];
?for?ii=-3:3
????for?jj=-3:3
????????for?kk=0
???????tb=b+repmat（[ii?jj?kk]*basissize（b1）1）;??????????
???????all=[all
??????????tb];
??????end
????end
?end

plot（all（:1）all（:2）‘.‘）
axis?equal

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

????I.A....????695323??2009-09-18?17:42??緊束縛近似方法在計算石墨烯能帶中的應用.pdf

?????文件????????311??2009-09-22?17:27??graphene.m

-----------?---------??----------?-----??----

???????????????695634????????????????????2