求解R中的代数Riccati方程

Question

我知道在Matlab中求解方程非常简单。但是有没有办法在R中解决它？我需要它用于我的卡尔曼滤波器实现。