modulus.m

plane wave expansion for 2D and 3D photonic crystal band gap calculations

function z=modulus(x)
%z=sqrt(x(:,1).^2+x(:,2).^2+x(:,3).^2);
[m n]=size(x);
if n==3,
	z=sqrt(x(:,1).^2+x(:,2).^2+x(:,3).^2);
else
   z=sqrt(x(1,:).^2+x(2,:).^2+x(3,:).^2);
end