expmint2.m

Computes the nested matrix exponential double integral:

  Integrate[expm(-i*A*(T-t))*B*
  Integrate[expm(-i*C*(t-x))*D*expm(-i*E*x),{x,0,t}],{t,0,T}]

This corresponds to the (1,3) block of the exponential of the auxiliary matrix (http://dx.doi.org/10.1109/TAC.1978.1101743).

Syntax

I=expmint2(spin_system,A,B,C,D,E,T)

Arguments

A,B,C,D,E - square matrices

   T         - upper limit of the outer integral

Outputs

I         - the integral as above

See also

expmint.m, dirdiff.m, trapdiff.m, Numerical infrastructure

Version 2.11, authors: Aditya Dev, Ilya Kuprov