function [mag,phs] = mybode(sys,freqVec)
  [mag,phs1] = bode(sys,freqVec*(2*pi));
  mag = squeeze(mag)';
  phs1= squeeze(phs1)';
  abs = exp(1i*phs1*pi/180);
  phs = angle(abs)*180/pi;
end