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