function  [mag,phase]=ghs(num,den,s)
%  Hadi Saadat, 1998
%  This function returns the magnitude and phase angle in radian of
%  the open-loop transfer function G(s)H(s) for the specified complex
%  variable s =