Francesco Bonazzi

I'm looking at the stats code... it definitely needs some overhaul.

Apart from bugs raising the errors... `JointPSpace` does currently not implement three methods necessary for this computation: 1. `where` 2. `conditional_space` 3. `probability`

SymPy does requires the code to have these conditions: 1. no line can end with whitespaces (tests are failing because of this), 2. exactly one newline at the end of...

> can I take up this issue? if you can find a way to solve this, feel free to open a PR.

@oscarbenjamin feel free to review :)

OK, I'll try to address all these problems in the next days.

One big question arises: what is the shape of nested array/matrix-like objects? For example, a matrix whose elements are matrix symbols should be shape `(n, n)` or `(n, n, n,...

> Whatever the answer is it should be consistent with the above or the above should be changed to be consistent with it. The problem is that such a change...

There is already a `Density` class in `sympy.stats.rv`. `Density.doit()` returns the PDF. One possibility the following: 1. in case of finite and discrete random symbols, `Density.doit()` should not evaluate or...

Please do not edit every distribution... that would be a huge mess.