niladridas

# Take a covariance matrix nxn
P = [2 1 1; 1 3 1;1 1 5]
C = eye(3,3)
R = [1 0 0; 0 2 0;0 0 300000]
K = P*C'*inv(C*P*C'+R)