Analytic functions with blackbody

See #3077 
Sounds OK to me, though let me also be a bit more explicit about the automation:
```
if flux_unit is None:
    flux_unit = FLAMBDA if getattr(in_x, 'unit', u.AA).is_equivalent(u.AA) else FNU
```