explde

Exponential-matrix solution of a linear differential equation based on the Magnus-exponent method.

This function numerically solves a first-order, linear, and homogeneous matrix differential equation with either constant or non-constant coefficients. The solution employs the Magnus exponent method (https://vixra.org/abs/2604.0065). Three approximation algorithm options are supported with residual errors proportional to the 3rd, 5th, and 7th power of the integration interval size. Paged operation is supported for parallel ensemble computations.