function mf = fit_mf(varargin)
%function mf = fit_mf(varargin)
%similar to fit_zpk, but in Matt's format

sys = fit_zpk(varargin{:});
mf = filtZPKs(sys.z{1},sys.p{1},sys.k);