石墨烯(graphene)中的狄拉克锥(Dirac cone)图形。直接按照能级公式绘图，公式参考[1]中 energy of the electrons 公式。
%matlab代码 clear; gamma0=2.8; a=2.46;%晶格常数的值 xlist=linspace(-2,2,100); ylist=linspace(-2,2,100); [kx,ky]=meshgrid(xlist,ylist); %能量与波矢kx,ky的函数关系 E=sqrt((gamma0.^2).*(1+4.*(cos(ky.*a./2).^2)+4.*(cos(ky.*a./2)).*cos(sqrt(3).*kx.*a./2))); surf(xlist,ylist,E)%能量为正的部分 hold on surf(xlist,ylist,-E)%能量为负的部分 shading flat
狄拉克锥图形
俯视图
参考：
[1] https://en.m.wikipedia.org/wiki/Graphene
[2] 石墨烯晶格