[bs] The Surface Plasmon

% The surface plasmon. Review Paper Ramakrishna : Metamaterials. 

% stacks.iop.org/RoPP/68/449

clear all
b0 = -pi;
b =pi;
a0 = -2*pi;
a =2*pi;
n = 0;
x = a;
i = sqrt(-1);
dx = 0.01*(a-a0);
dz = 0.01*(b-b0);
while x>=a0;
   n = n+1;
   m = 0;
   z = b;
   T = 0;
   while z>=0;
   m = m+1;
   dT =-real(exp(i.*x-z)).*dz;
   T = T-dT;
   Xaxis(m,n) = x;
   Yaxis(m,n) = z;
   Zaxis(m,n) = T;
   z=z-dz;
   end
   while z>=b0;
   m = m+1;
   dT =real(exp(i.*x+z)).*dz;
   T = T-dT;
   Xaxis(m,n) = x;
   Yaxis(m,n) = z;
   Zaxis(m,n) = T;
   z=z-dz;
   end
   x=x-dx;
end

%pcolor(Xaxis,Yaxis,Zaxis)
meshc(Xaxis,Yaxis,Zaxis)
%mesh(Xaxis,Yaxis,Zaxis)
%grid on
xlabel('X (tangential)')
ylabel('Z (normal)')
zlabel('Field density')
title('The surface plasmon')