%
% This defines the nonlinear relationship between heat capacity 
% and temperature (Eq. 2.3.2). The format is suitable for the  
% bisection and tangent methods:
%
function f = heat(T)
f = 0.716-4.257*T*10^(-6)-15.04/sqrt(T);
