1 2 3
background = function(mu0,day){ return (mu0 * (1 + 0.43 * sin((2*pi)/(365.24) * day + 4.36)) ) }