function h = steamEnthalpy(p,T,c)
% steamEnthalpy  Use an interpolant to compute the enthalpy of steam
%                The interpolant is of the form
%
%                    h = c_1 + c_2 * p + c_3 * T + c_4 * p * T
%
%                where h is the enthalpy, p is the pressure, and T is the
%                temperature
%
% Synopsis:  h = steamEnthalpy(p,T,c)
%
% Input:  p,T = scalar or vector values of pressure and temperature
%         c   = vector of coefficients for the interpolant
%
% Output:  h = scalar or vector values of enthalpy at the given p and T

h = c(1) + c(2)*p(:) + c(3)*T(:) + c(4)*p(:).*T(:);