The symsum(M) function computes the element-wise symbolic sum of M.  The result, R, is formed as R[i,j] = symsum(M[i,j], v, a, b).

F = symsum(f) is the indefinite symbolic sum of the scalar f. If f is a constant the variable of summation is n.

symsum(f,v) is the indefinite symbolic sum of the scalar f with respect to v.

symsum(f,a,b) is the definite symbolic sum of the scalar f in the range a..b.

symsum(f,v,a,b) is definite symbolic sum of the scalar f in the range a..b with respect to v.

$\mathrm{with}\left(\mathrm{MTM}\right)\:$

$M≔\mathrm{Matrix}\left(2\,3\,'\mathrm{fill}'=\exp \left(x\right)+3x^2+5\right)\:$

$\mathrm{symsum}\left(M\right)$

$\left[\begin{array}{ccc}\frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}& \frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}& \frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}\\ \frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}& \frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}& \frac{x\left(2x^2-3x+11\right)}{2}+\frac{{\ⅇ}^x}{\ⅇ-1}\end{array}\right]$