title: Symbolic integration with Python
body: <p>I'd like to compute a distribution in a symbolic way with Python.
  The distribution is</p>
  
  <pre><code>f(z) = \frac{1}{(b-a)^2}\int^b_a z\mathbb{1}_{[a,b]}(zw)dw
  </code></pre>
  
  <p>How to compute  this integral knowing that </p>
  
  <pre><code>[a,b]\bigcap [a/z,b/z] \neq \emptyset
  </code></pre>