function haxes = mpcplot(haxes,Phi,Freq)
%
%  Name: mpcplot
% 
%  Usage: haxes = mpcplot(haxes,Phi,Freq)
% 
%   Description:
%     Plots the Modal Phase Colinearity of the current complex Modal Set
%
%   inputs:
%	Phi = complex modal matrix
%
%   Following workspace variables must be defined:
%
%     

%
%	
%  Version SWD960716
%
%  Author(s): Scott W. Doebling, doebling@lanl.gov
%
%  Institution: Los Alamos National Laboratory,
%		ESA-EA
%
%  Copyright (c) 1996, The Regents of the University
%     	of California. (See top level of code for
%	full copyright statement.)


%
%  Check that the mode is complex
%

if isreal(Phi),
 	disp('ERROR: Cannot plot Modal Phase Colinearity for real modes')
else

%
% 	Compute MPC Value
%
	MPC = mpc_comp(Phi);

%
%	Draw MPC plot
%
	axes(haxes)
	polar(angle(Phi),abs(Phi),'b')
	title(['Mode at ',num2str(Freq),' Hz, MPC = ',num2str(MPC(1))])

end

return
