gold.m

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

function vEpsilon = gold(vLambda_in)
load metals
cs = csaps(vLambda, vAu);
vEpsilon = ppval(cs, vLambda_in);