function [K] = Kinlet(y,params) Uhot = params(1); Ucool = params(2); Thot = params(3); Tcool = params(4); kg = params(5); kw = params(6); L = params(7); tw = params(8); h = params(9); if (y > 0), K = kg; elseif (y < -tw), K = kg; else, K = kw; end