% -*- texinfo -*-
% @deftypefn {Function File} {[@var{C}]} = BIM2Areaction(@var{mesh}, @var{delta}, @var{zeta})
%
% Builds the matrix for the discretization of the LHS
% of the equation
% @iftex
% @tex
% $ \delta \zeta u = f $
% @end tex
% @end iftex
% @ifnottex
% @var{delta} * @var{zeta} * u = f
% @end ifnottex
%
% Input:
% @itemize @minus
% @item @var{mesh}: PDEtool-like mesh structure with required fields 'p', 'e', 't'.
% @item @var{delta}: element-wise constant scalar function.
% @item @var{zeta}: piecewise linear conforming scalar function.
% @end itemize
%
% @seealso{BIM2Arhs, BIM2Aadvdiff, BIM2Cmeshproperties}
% @end deftypefn