clear;clc

e2p=0;
t=-0.5;
s=0.13;
a=2.46;

n=-2:0.01:2;
[kxky]=meshgrid（n）;
w=sqrt（1+4*cos（sqrt（3）*kx*a/2）.*cos（a*ky/2）+4*cos（a*ky/2