Analytic functions with blackbody

See #3077 
I wonder if we should not raise an exception on finding any input value smaller than 0 (0 itself is OK, just makes things go to infinity but the end result is 0 anyway).