title: unitary diagonlization in Mathematica
body: <p>is there a built-in function or made already code in mathematica that computes a unitary diagonlization of a given matrix?</p>
  
  <p>I.e, if the matrix is: A then I want to find P and D such that:</p>
  
  <p>A=P* D P , where PP* = I, and * is the complex conjugate transpose.</p>