%
%   HT1D_SS3A.M  Function evaluation for ODE step in HT1D_SS3.M 
%     
%   given differential eqn.:  z' = Az + b  (with time dependent BCs on right)
%
      function zp = ht1d_ss3a(t,z)
      global A b mu
      bb = b;  bb(length(bb)) = bb(length(bb))*exp(-mu*t);
      zp = A*z + bb;
%
%   end of function


