Tests if a matrix is normal.
Ais a matrix
εis a positive number
A is a matrix, then
IsNormal(A, ε) returns
A is normal, that is, iff
A commutes with its conjugate transpose
A*. Floating-point comparisons are made with epsilon
ε; if omitted, a default (non-zero) epsilon is used.