maincount.m

for the parameter estimation and plot the temperature profile for the exchanger

clear
clc

solinit=bvpinit([0:0.01:1],[0.5 1.8 0.9 2.5])
sol=bvp4c(@heatexcount,@heatexcountbc,solinit);

xint=[0:0.01:1];
yint=deval(sol,xint);

plot(xint,yint)
legend('xgc','tgc','xgh','tgh');