function dT_dt = tempDE(t, T)

k = 0.02;
Te = 20; % 20 degree environment

dT_dt = -k * (T - Te);
