function y = f_exp_trig(t,p)
  y = p(1)*exp(-p(2)*t).*cos(p(3)*t + p(4));
