这个程序可以用于三角晶格空气孔光子晶体的能带计算。给出了解析的折射率分布

u1=[1,0,0]

u2=[1/2,sqrt(3)/2,0]

u3=[0,0,1]

-This code can be used in band structure of triangular lattice rod photonic crystal. It give out the dielectric distribution of circular air holes.