Analytic functions with blackbody

See #3077 
I think B_lambda is used a bit more, since most people work with wavelengths. But I realise I am biased by optical observations... Anyway, if we stick with `AA` for the input, I think we should retrun `Blambda`. It may be slightly more elegant to then calculate B_lambda directly (don't quite understand why this would be less stable).