%
% This defines the nonlinear relationship between heat capacity
% and temperature (Eq. 2.3.5). The format is suitable for the
% Wegstein method and is a divergent expression using direct
% substution
function f = cp(T)
f = (0.716-(15.04/sqrt(T)))/(4.257*10^(-6));