Numerical Evaluation of Integral Involving Matrices

Hi,

I have come across an integral involving matrices while trying to develop a closed form solution for the problem I have been working on.

I am quite confident about my closed form expression. However, my closed form solution differs from solution obtained using simulation. So, I believe that I am wrong in evaluating that integral.

The integral that I am evaluating at t= 0.0300 sec (end point of Gaussian pulse :

I = Integral (I1. dt)

Where I1 = exp(-I0) ; I0 = Integral((A_mo + Alpha_mo * G(t)) dt)

Both A_mo and Alpha_mo are constant matrices (6 x 6) and G(t) is a Gaussian function.

I am also attaching my code where my attempt as well as structures of all matrices are given.

I would like to evaluate this integral for three different value of A_mo (also included in code)

I would appreciate some help regarding this.

Regards, DK