Invert the matrix lam*t(lam)+diag(u) via Woodbury identity
woodbury(lam, u)
lam
p by k matrix
u
p dimensional vector
The matrix inverse