Analytic functions with blackbody

See #3077 
:+1: for `planck_blackbody`--still pretty short and reasonably unambiguous.