Writing a new distribution

I’m exactly trying to write power-law and Schechter distributions!! :slight_smile:
They are so common in astronomy and should be available in pymc3…