interpolate.m

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

function varray=interpolate(vbegin,vend,n,inclend)
varray=zeros(n+1,3);
step=(vend-vbegin)/n;
for i=1:3,
   if step(i) ==0,
      varray(:,i)=vbegin(i);
   else
      varray(:,i)=(vbegin(i):step(i):vend(i))';
   end
end
if inclend==0
   varray(n+1,:)=[];
end