function out = k_steel_dT(T) % Derivative with respect of T of the function k_steel. out = -0.115+2.32e-4*T-12.75e-8*T.^2;