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

u1=[1,0,0]

u2=[0,1,0]

u3=[0,0,1]

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