function [v,d]=feig(w,f)
%FEIG   Eigenvalues and eigenvectors of MVFR matrix.
%       FEIG(W,F) returns the eigenvalues of an MVFR
%       matrix one row per frequency. W is the associated
%       frequency vector.
%
%       [V,D]=FEIG(W,F) returns the eigenvectors V and
%       the diagonal matrices of eigenvalues, D, in MVFR form.
%       See also FSVD

%       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
k=1:m;         % vector of rows of each matrix in F

if nargout==2   % both required eigenvectors and values required
  d=zeros(lw*m,n);   % Set up output matrix for values
  v=zeros(lw*m,n);   % Set up output matrix for vectors
  for i=1:lw     % for each frequency
     [vec,val]=eig(f(k+(i-1)*m,:));
      v(k+m*(i-1),:)=vec;
      d(k+m*(i-1),:)=val;
  end % for i=1:lw
else   % only eigenvalues required
  v=zeros(lw,n);   % Set up output matrix for values
  for i=1:lw       % for each frequency
     v(i,:)=eig(f(k+(i-1)*m,:)).';  % transpose eigenvalues
  end % for i=1:lw
end  % else