function s=fcsod(w,f)
%FCSOD  Column sum of abs(off-diagonals) of MVFR matrix.
%       FCSOD(W,F) returns the column sum of the absolute
%       values of off diagonals of the component matrices of
%       the MVFR matrix, F, one row per component matrix.
%       W is the associated frequency vector.
%       See also FRSOD,FCDOM,FRDOM

%       Dr M.P. Ford 4th August 1987
% Copyright (c) 1987 by GEC Engineering Research Centre & Cambridge Control Ltd


[m,n]=fsize(w,f);
lw=length(w);
if m~=n
  error('Not a square system');
end
s=ones(lw,n)*inf;   % Set up output matrix
k=1:m;              % vector of rows of each matrix in F
for i=1:lw          % for each frequency
  fm=abs(f(k+(i-1)*m,:));
  s(i,:)=sum(fm)-diag(fm)';
end % for i=1:lw