%Exercise for Section 4.4, overerdamped harmonic motion
% x'' + 6*x' + 4*x=0
% with solution y(t)=c1*exp((-3-2*sqrt(5))*t) + c2*exp((-3=2*sqrt(5))*t) 
% The damping constant is c=6/2=3, 
% the natural frequency omega0=sqrt(4)=2.
% The constants depending on the initial conditions
c1=1;c2=2;

t = linspace(0,5,50);
y = c1*exp((-3-sqrt(5))*t) + c2*exp((-3+sqrt(5))*t) ;
Amplitude=norm([c1 c2]);

figure(2)
hold all
plot(t,y,'-r')
legend('y(t)')
title('Overdamped harmonic motion c=3>2=\omega_0')
