function yp = murder(t,y)

% Model for Newton's Law of Cooling
% Parameters for the model 
k = log(4/3); Te = 22;
yp = -k*(y - Te);
end

% [t,y] = ode23(@murder,[-9,12],y0);
