Analytic functions with blackbody

See #3077 
Actually, since `HC` is only used once, I would just calculate it in the BB function; alternatively, convert to `Hz` a little earlier, and just calculate `const.h * nu` instead of `const.h * const.c / w`.