function result = deltaTm(Th_in, Th_out, Tc_in, Tc_out)
%Ideal counter-current heat exchanger
result = ((Th_in-Tc_out)-(Th_out-Tc_in))/log((Th_in-Tc_out)/(Th_out-Tc_in));

end