function [z dT1 dT2] = Tprof(T0,Th1,Th2,w1,w2,wh1,wh2,UA1,UA2)

z=(0.0001:0.0001:1)';

R1=(w1-wh1)./(w1*wh1);
R2=(w2-wh2)./(w2*wh2);

dT1 = (Th1-T0).*exp(R1*UA1.*z);
dT2 = (Th2-T0).*exp(R2*UA2.*z);