Riccati equation