function y = sine_fn(t)

y = sin(0.1*t).*exp(-0.001*t) + 3*cos(0.3*t).*sin(0.05*t).*exp(-0.003*t) + 5*sin(0.5*t).*cos(0.01*t).*exp(-0.0005*t);