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