Analytic functions with blackbody

See #3077 
I did a quick test, passing in `Column` gives wrong answer because it is not being converted to Hz properly. I think `Column` support for unit conversion is out of scope of this PR.