% 
% Written by:
% -- 
% John L. Weatherwax                2005-08-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

omega = 1.6; 
2/(1+omega^2)

T=1; 
T*( 1 - exp(-2*T) )/( 1 - 2*exp(-T)*cos( omega*T ) + exp( -2*T ) )

T=0.1; 
T*( 1 - exp(-2*T) )/( 1 - 2*exp(-T)*cos( omega*T ) + exp( -2*T ) )