Computes matrix exponential integrals of the following general type:
Matrix A must be Hermitian. For further info see the paper by Charles van Loan (http://dx.doi.org/10.1109/TAC.1978.1101743).
A,B,C - the three matrices involved in the integral T - integration time
R - the resulting integral
The auxiliary matrix method is massively faster than either commutator series or diagonalisation.