function SharifPWE_3DplotH(N,RES,r,epa,epb)
%
%The ouput is the 3D H-polarization band structure (3 first band surfaces)
%of a square-lattice photonic crystal with unit cell consisting of a
%substrate of "epb" relative permittivity with circular rods(holes) of
%"epa" relative permittivity and radius "r".
%
%Usage is
% SharifPWE_3DplotH(N,RES,r,epa,epb)
%
% N: N in PWE algorithm (number at which the summations are truncated)
% RES: number of the divisions on every side of BZ#1
% r: normalized radius of the circle of holes/rods (0