The numerical solution of the general matrix differential equation

for

is considered where A
1and A
2are stable matrices. The algorithm proposed requires only

words of memory (for large

) and converges in approximately

s where μ is the multiplication time of the digital computer, and

where

. The algorithm is particularly suitable for systems where

is large (e.g,

).