I need help implementing this bit of code into matlab. I have a known matrix lest call it Z. and I have another unknown matrix called H. The matrices are related by the following relationship. Z = HTH How do I solve for H in matlab? If it is relevant Z is a make up of other known matrices and that code in matlab is: Z = A + 2*eye(n) where A is the initial known matrix and n is just an integer. UPDATE: Yes A and Z will always be positive.