Spring_Mass.m

%This script solves the ODE Equation for the
%Spring Mass System given in the HW

%Values of System 
stiff = [10,20,20,10];
mass= [1;2;4;];

%ODE Equations
f1=[-1*(stiff(1)+stiff(2)), stiff(2), 0];
f2=[stiff(2), -1*(stiff(2)+stiff(3)), stiff(3)];
f3=[0,stiff(3), -1*(stiff(3)+stiff(4))];

K= [f1;f2;f3];
K= -1*K;

M= diag(mass);

K_tilde= (M^(-1/2))*(K)*(M^(-1/2));

%Final lambda and omega
e=eig(K_tilde)

w= sqrt(e)