An efficient algorithm for the computation of the integral of the state transition matrix over the interval [

] is developed, It is assumed that there exists some

such that the state transition matrix and the integral of the state transition matrix can be computed as accurately as required over the interval [

].