function phi = efld_mag2( x )
%
%  return the magnitude squared of the electric field at point x

e = efld( x );
phi = e(1)^2 + e(2)^2;  %   phi = |E|^2 = E_x^2 + E_y^2
