Analytic functions with blackbody

See #3077 
This should also be renamed to something more specific - maybe ``blackbody_sample_wl``?