%function L = getL(K,n)
function L = getL(K,n)
A = K(1:n,1:n);
C = K((n+1):end,(n+1):end);
B = K(1:n,(n+1):end);
L = A - B*(C\B');