function [orix,oriy] = showquiv(Es,Eo)
%
[rows,cols] = size(Es);
orix = zeros(rows,cols);
oriy = zeros(rows,cols);
for r = 1:rows
   for c = 1:cols
      if Es(r,c)>0
         orix(r,c) = cos(Eo(r,c));
         oriy(r,c) = -sin(Eo(r,c));
      end     
   end
end
quiver(orix,oriy);
