title: Numpy transpose multiplication problem
body: <p>I tried to find the eigenvalues of a matrix multiplied by its transpose but I couldn't do it using numpy.</p>
  
  <pre><code>testmatrix = numpy.array([[1,2],[3,4],[5,6],[7,8]])
  prod = testmatrix * testmatrix.T
  print eig(prod)
  </code></pre>
  
  <p>I expected to get the following result for the product:</p>
  
  <pre><code>5    11    17    23
  11    25    39    53
  17    39    61    83
  23    53    83   113
  </code></pre>
  
  <p>and eigenvalues:</p>
  
  <pre><code>0.0000
  0.0000
  0.3929
  203.6071
  </code></pre>
  
  <p>Instead I got <code>ValueError: shape mismatch: objects cannot be broadcast to a single shape</code> when multiplying <code>testmatrix</code> with its transpose.</p>
  
  <p>This works (the multiplication, not the code) in MatLab but I need to use it in a python application.</p>
  
  <p>Can someone tell me what I'm doing wrong?</p>