function [m,p,w,FocusInfo,sdm,sdp] = magphaseresp_(sys,grade,wspec)
% Compute mag and phase data of frequency response
% (magnitude in abs value, phase in radians)

% Copyright 1986-2014 The MathWorks, Inc.
[m,p,w,FocusInfo] = freqresp(sys.Data_,grade,wspec,false);
m = permute(m,[2 3 1]);
p = permute(p,[2 3 1]);
sdm = []; sdp = []; % default for systems with no covariance info