pop.m

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

function [out, vLIFO] = pop(vLIFO, N)
out = vLIFO(1:min(end, N));
if length(vLIFO) == N
    vLIFO = [];
else
    vLIFO = vLIFO(N+1:end);
end