%EA3-Fall 1999
% 
%%---------------------------------------------------------
%%DETERMINATION OF NORMAL MODES
%%---------------------------------------------------------


clear all;
close all;


%Define Y (reduced set of state variables for second-order form)
%Y(1) = rsp1
%Y(2) = rsp2

%Enter parameters of system
K1 = 8e6;
m2 = 50;
K3 = 8e6;
m5 = 800;

%Enter B matrix (Ydotdot = BY -- second order form)

B = [-(K1+K3)/m2		K3/m2;
	  K3/m5		-K3/m5];
	  
disp('The B-matrix is:'), disp(B)
%calculate eigenvalues and eigenvectors
[Evector,Evalue] = eig(B);

disp('The Eigenvalues are:'), disp(diag(Evalue)')
freq=sqrt(-diag(Evalue));
disp('The normal mode angular frequencies are:'), disp(freq')

disp('The Eigenvectors are:'), disp(Evector)