title: NumPy - Create Upper Triangular Matrix with Diagonals of Increasing Power
body: <p>I'd like to create a square upper triangular matrix that is defined as follows for some float <code>c</code> and some dimension <code>N</code>:</p>
  <pre><code>[[1 , c , c^2, ... c^N],
   [0,  1 ,   c, ... c^{N-1}],
   [0,  0 ,   1, ... c^{N-2}],
   .
   .
   .
   [0,  0 ,   0, ....    1]]
  </code></pre>
  <p>For concreteness, if <code>N=2</code>, then the matrix should be</p>
  <pre><code>[[1, c],
   [0, 1]]
  </code></pre>
  <p>And if <code>N=3</code>, then the matrix should be:</p>
  <pre><code>
  [[1, c, c^2],
   [0, 1,   c],
   [0, 0,   1]]
  </code></pre>
  <p>How can I do this?</p>