
 t=0:0.1:30;
 A=1; % Amplitud
 tau=3; % tiempo decaimiento
 Omega=1;  %frec. angular corregida
 t0=0;   % fase


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