equimap.m

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

function map=equimap(n)
x = linspace(0, 1, ceil(n^(1/3)));
[mX, mY, mZ] = ndgrid(x);
map = [mX(:), mY(:), mZ(:)];
%plot3(map(:,1), map(:,2), map(:,3), '.')