Analytic functions with blackbody

See #3077 
Consensus was *not* to silence Numpy (see above).