% -*- texinfo -*-
% @deftypefn {Function File} {[@var{b}]} = BIM2Arhs(@var{mesh}, @var{f}, @var{g})
%
% Constructs the RHS for the DAR problem
% @iftex
% @tex
% $ -div ( \alpha \gamma ( \eta \vect{\nabla} u - \vect{beta} u )) + \delta \zeta u = f g $
% @end tex
% @end iftex
% @ifnottex
% -div (@var{alpha} * @var{gamma} (@var{eta} grad u - @var{beta} u )) + @var{delta} * @var{zeta} u = @var{f}*@var{g}
% @end ifnottex
%
% Input:
% @itemize @minus
% @item @var{mesh}: PDEtool-like mesh with required field 'p', 'e', 't'.
% @item @var{g}: element-wise constant scalar function.
% @item @var{f}: piecewise linear conforming scalar function.
% @end itemize
%
% @seealso{BIM2Areaction, BIM2Aadvdiff, BIM2Cmeshproperties}
% @end deftypefn