function l = lambda_t(r, f, t);
% LAMBDA_T  Periodic function with avarage rate r and frequency f 
     
l = r*(1+cos(2*pi*f*t));