function res = newtonWrapper(t)
% a wrapper for newtonCooling to make it a function of t
res = newtonCooling(93.5, 24.5, 0.0426, t);