function temp = coffeemodel(k,Tzero,Tenv,time)
% coffeemodel(k,T0,Tenv,time)
% Simple Newton's Law model of coffee cooling
temp = Tenv + (Tzero-Tenv).*exp(-k.*time/1000); 