Errors working with a largematrix

[ptc-5412730]

Z is a 51 row X 3 column matrix. I need to obtain the 51 X 51 matrix H = Z(Z'Z)^(-1)Z', where ' denotes transpose and ^(-1) means inverse. H is idempotent, so as a check, either H'H or HH' should yield an identity matrix. Straightforward evaluation of H evidently yields an incorrect result, as H'H is not an identity matrix.. Can MathCad handle matrices this large? Or are there special procedures for dealing with such matrices?

[Werner_E]

Matrices of that size should be no problem for Mathcad or for Prime and there are no special treatments.

Why don't you tell us which version of Mathcad or Prime you are using and why not attach your worksheet showing what you think is wrong?

[Werner_E]

Sorry, I just read your post more carefully.

H is idempotent, so as a check, either H'H or HH' should yield an identity matrix.

Why? Doesn't idempotent simply mean that H^2 = H?

This is true for H: