%
% 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));

