T = 298:2000;

kappa = 8.8527 + 7.082E-3 * T + 2.5329E-6 * T.^2 + 2.9918E3 ./ T;

plot (T, kappa);

xlabel ('T [K]');
ylabel ('\kappa [W/m/K]');
title ('Zirconium Thermal Conductivity (best fit, std. dev. 5 to 9.5%)');

print -deps zirconium_kappa.eps

