Analytic functions with blackbody

See #3077 
I'm fine with either, but think the units of `in_x` and the flux should be consistent (i.e., either wavelength and flux per unit wavelength or frequency and flux per unit frequency). We could even make this automatic...