Analytic functions with blackbody

See #3077 
Okay, so I don't have to explicitly catch overflow/underflow anymore? The logic was inherited from codes that did not use Numpy.