title: Python : generating random numbers from a power law distribution
body: <p>I want to draw a random variable between 2 to 15, from a power law distribution with negative exponent (a = -2). I found the following :</p>
  
  <pre><code>r = scipy.stats.powerlaw.rvs(a, loc = 2, scale = 13, size = 1000)
  </code></pre>
  
  <p>But it does not take negative numbers for a.</p>
  
  <p>Anyone knows a way out?</p>