Another numerical solution of the general matrix differential equation

for X is considered without any stability condition for A
1and A
2. Like Davison\´s method, the proposed algorithm requires only some n
2words of memory and

multiplications where

and

. This new approach is well suited to solve large and possibly unstable systems. We take the opportunity to run the differential equation for various D. A very efficient technique follows to design the so-called receding horizon control problem.