Johnson SU Distribution

There are two shape parameters $a\in \mathbb{R}$ and $b>0$, and the support is $x\in \mathbb{R}$.

Implementation: scipy.stats.johnsonsu