function sm = StabMode(F,Vt,LLD,LON,SBA,DIM)
%	Stability and Modes of Motion
%	=============================
	sm	=	1;
		[Eigenvectors,EigenvalueMatrix]	=	eig(F)
		kv	=	length(F);
		'Eigenvector Amplitudes'
		for k = 1:kv		% Eigenvector Amplitudes
			a	=	Eigenvectors(:,k);
			Magnitude	=	real(sqrt(ctranspose(a) * a));
			k
			EvectorAmplitude	=	sqrt(a .* conj(a)) / Magnitude
		end
		'Velocity-weighted Eigenvector Amplitudes'	
		for k = 1:kv		% Velocity-Weighted Eigenvector Amplitudes
			a	=	Eigenvectors(:,k);
			if LLD >= 1
				if LON >= 1			% Reduced-Order Model
					if SBA >= 1		% Longitudinal Model
						if DIM >= 6	% Stability/Hybrid Axes
							a	=	[a(1)/Vt;a(2:4);a(5:6)/Vt];	
														% 6th-Order										
						else
							a	=	[a(1)/Vt;a(2:4)];	% 4th-Order						
						end
					else
						if DIM >= 6	% Body Axes
							a	=	[a(1:2)/Vt;a(3:4);a(5:6)/Vt];	
														% 6th-Order										
						else
							a	=	[a(1:2)/Vt;a(3:4)];	% 4th-Order						
						end			
					end	
				else
					if SBA >= 1		% Lateral-Directional Model
						if DIM >= 6		% Stability/Hybrid Axes
							a	=	[a(1:5),a(6)/Vt];	% 6th-Order										
						else
							a	=	a;					% 4th-Order						
						end
					else
						if DIM >= 6		% Body Axes
							a	=	[a(1)/Vt;a(2:5);a(6)/Vt];	% 6th-Order										
						else
							a	=	[a(1)/Vt;a(2:4)];	% 4th-Order						
						end			
					end			
				end
			else
				a	=	[a(1:6)/Vt;a(7:12)];		% Original Model					
			end
			Magnitude	=	real(sqrt(ctranspose(a) * a));
			k
			EvectorAmplitude	=	sqrt(a .* conj(a)) / Magnitude
		end