private void GetQ22(IMatrix q21, SymmetricMatrix Q)
 {
     Q22 = Q.Plus(Q.Multiply(q21).Multiply(Q11).Multiply(q21.Transposition).Multiply(Q));
 }