Revision history [back]

t, s = var(' t s ')

A = matrix([[3*t,t^2],[-2*t, t]])

B = A.apply_map(lambda e: integrate(e,t,0,s))

B.exp()

I Think this is what you are looking for. Note that simple 'integrate' function won't integrate the matrix.