function  y = fung691 (t,q)
%------------------------------------------------------------------------
% Heat Exchanger 
%------------------------------------------------------------------------
   if t > q(3)
      a = -(t - q(3))/max([q(2),0.1]);
      y = q(1)*(1 - exp(a));
   else
      y = 0;
   end