
 t=0:0.1:30;
 A=1; 					% Amplitud
 tau=10; 				% tiempo decaimiento
 T=7; 					% periodo
 w0=2*pi/T;				% frecuencia natural

 Omega=sqrt(w0*w0-1/(4*tau*tau));  	%frec. angular corregida
 t0=0;   				% fase



 phi=A*exp(-t/(2*tau)).*cos(Omega*(t-t0));
 plot(t,phi,'b')

 