function [ abs, phase ] = mnbode_list( sys, freq, varargin )
%mnbode_list make gain, phase list of system with frequency list
%   [abs, phase] = mnbode_list(sys,freq)
%           options: 'c','complex' : output complex list on abs
%                    'd','delay'   : include time delay

[temp_abs,temp_phase] = bode(sys,freq*2*pi);
abs = squeeze(temp_abs);
phase = squeeze(angle(exp(1i*temp_phase*pi/180)))*180/pi;
for i=1:nargin-2;
    if strcmp(varargin{i},'d')||strcmp(varargin{i},'delay')
        phase = squeeze(angle(exp(1i*(squeeze(temp_phase)*pi/180-varargin{i+1}*freq*2*pi))))*180/pi;
    end
end
for i=1:nargin-2;
    if strcmp(varargin{i},'complex')||strcmp(varargin{i},'c')
        abs = abs.*exp(1i*phase*pi/180);
    end
end
end