function misa=fmisalg(w,f)
%FMISALG        Calculates misalignment angles
%       FMISALG(W,F) returns a matrix whose columns are
%       the misalignment angles of F, an MVFR matrix, at
%       each frequency. There is one row for each frequency.
%       W is the associated frequency vector.

%       Dr M.P. Ford 14th September 1987
% Copyright (c) 1987 by GEC Engineering Research Centre & Cambridge Control Ltd
%       MRN0019
[m,n]=fsize(w,f);
if m~=n
   error('System is not square')
end
lw=length(w);
% set up output matrix
misa=zeros(lw,n);

id=eye(n);
infinity=ones(1,n)*2/eps;   % Note +and- inf is always returned
			    % as the minimun from the MIN function
k=1:m;

for j=1:lw
  [v,d]=eig(f((j-1)*m+k,:));
  vte=abs(v'*id);
  nv=sqrt(abs(diag(v'*v)));
  for i=1:n
     nvte(i,:)=vte(i,:)./nv(i);
  end
  for i=1:n
     [dmax,di]=max(nvte);
     [dmax,dj]=max(dmax);
     misa(j,dj)=nvte(di(dj),dj);
     nvte(di(dj),:)=-infinity;
     nvte(:,dj)=-infinity';
  end     % for i=i:n

end  % for j=1:lw
misa=acos(misa).*(180/pi);