function [g] = clxbode(Z1,Z2,Z3,Z4,Z5,Z6);
%[G] = CLXBODE(SS_,IU,W) or
%[G] = CLXBODE(A,B,C,D,IU,W) produces a complex gain instead of an 
%      absolute gain.

% R. Y. Chiang & M. G. Safonov 6/85
% Copyright (c) 1988 by the MathWorks, Inc.
% All Rights Reserved.
% ----------------------------------------------------------------
%

inargs = '(a,b,c,d,iu,w)';
eval(mkargs(inargs,nargin,'ss'))

[rc,cc] = size(c);
cs = max(size(w));
%
b = b(:,iu);
d = d(:,iu);
s = sqrt(-1) * w;
g = ltifr(a,b,s);
g = c * g + diag(d) * ones(rc,cs);
%
% ----- End of CLXBODE.M ---- RYC/MGS 6/85 %









