add solve-riccatti function but not convergence