main.m

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

clear
clc

[time,sol]=ode45(@heatexco,[0:0.01:1],[0 0 0 0]);

sol
plot(time,sol)
legend('xgc','tgc','xgh','tgh');