one reason is:
a matrix calculation reduces the computation of composition of linear transformations, as well as the computation of image elements under a linear transformation, to arithmatic operations in the underlying field. that is:
conceptual--->numerical.
sometimes, this is preferrable for getting "actual answers" in a physical application, where some preferred basis (coordinate system) might already be supplied.
for example, the differentiation operator is a linear transformation from Pn(F) to Pn(F).
actually "computing a derivative" IS just computing the matrix product [D]B[p]B = [p']B:
for n = 2, and F = R, we have for the basis B = {1,x,x2}, that [D]B=
[0 1 0]
[0 0 2]
[0 0 0],
or that if p(x) = a + bx + cx2,
p'(x) = b + 2cx.
of course, this would be just as easy using D(p) = p' using the calculus definition,
but it's not so clear what happens if you want to use THIS basis: {1+x,1-x,1-x2}, using the calculus definition, whereas the matrix form makes it transparent.