cp = C1 + C2*T + C3*T^-2
function specificHeatCapacity input Real T; output Real cp; end specificHeatCapacity;