Matlab code for the solution to Riccati matrix difference equations associated with the Kalman filter