The time-varying and incommutable character of the coefficient matrix of periodically time-varying linear systems are the bottleneck of the design for high precision direct integration methods.