function [T] = Tinlet(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 > h), T = Thot; else, T = Tcool; end