function out = k_steel(T) % k_steel defines a temperature dependent thermal conductivity for a steel % material. out = 71.12-0.115*T+1.16E-4*T.^2-4.25E-8*T.^3;