function linearTE_g

global ks ind kn Ny Nx E dx dy ind Hr

handle=findobj('Tag','vid');
nr=str2num(get(handle,'String'));

for n=1:(Ny+1)
  Hm(n,:) = Hr(1+(n-1)*(Nx+1):n*(Nx+1),ind(nr))';
end

for n=2:Nx
  for m=2:Ny
    Ey(m,n) = 0.5*(Hm(m,n+1)-Hm(m,n-1))/dx(n);
    Ex(m,n) = 0.5*(Hm(m+1,n)-Hm(m-1,n))/dy(m);
  end
end
Emax = max( max( sqrt( abs(Ex*Ey'))));

figure;  surf( Hm/max(max(abs(Hm))));
title(['Hz, kn = ',num2str( ks( nr))]);
figure;  quiver( Ex/Emax, Ey/Emax);
title(['Et, kn = ',num2str( ks( nr))]);
figure;  quiver( Ey/Emax, Ex/Emax);
title(['Ht, kn = ',num2str( ks( nr))]);