function r = beta_SEIR(t)

   global alpha beta0 beta1 gamma mu

   r = beta0*(1+beta1*cos(2*pi*t));
