title: skew normal distribution in scipy
body: <p>Does anyone know how to plot a skew normal distribution with scipy? 
  I supose that stats.norm class can be used but I just can't figure out how.
  Furthermore, how can I estimate the parameters describing the skew normal distribution of a unidimensional dataset?</p>