heatexcountbc.m

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

function bc=heatexcountbc(ya,yb)

bc=[ya(1);ya(2);yb(3);yb(4)];