aluminum.m

利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤

function vEpsilon = aluminum(vLambda_in)
global C
w = C./vLambda_in*2*pi;
h_bar = 6.626068e-34/2/pi;
e = 1.60217653e-19;
Wp = 15.565*e/h_bar;
gamma = 0.608*e/h_bar;
vEpsilon = 1 - Wp^2./(w.^2-j*gamma*w);