function c = dinamicac_unod(c, D, ~, alphac, U, gammac, cstar, qc)
% Aggiornamento della sostanza chimica (solo da U)
% ~: placeholder per compatibilitÃ  (era I, ora ignorato)

c_right = [c(2:end); 0];
c_left = [0; c(1:end-1)];

% Produzione chimica solo da cellule U
c = c + D * (c_left + c_right - 2 * c) + (alphac * 0 + gammac * U) .* (cstar - c) - qc * c;
c(1) = c(2);
c(end) = c(end - 1);
