Matrix Riccati code



I have written some code to solve the matrix riccati differential equation,
which is useful for control systems.
It is proof of concept, and has several bugs.
It is attached.  I would appreciate any comments.
The biggest issue is setting the solution constants,
which for control systems, the solution is zero at time T.
Dan Stanger

Attached file: odematriccati.mac
Attached file: array_symbols.mac