function [x] = k_ice(T) % temperature in °C and k in W/m/K temp = [0 -10 -20 -30 -40]; k = [2.25 2.25 2.25 2.25 2.25]; x = interp1(temp,k,T,'pchip'); end