equimap.m

来自「利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤」· M 代码 · 共 8 行

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), '.')