Analytic functions with blackbody

See #3077 
Actually, I would immediately convert the input to `Hz`, since where you use the wavelength, the frequency is a better alternative, and you use the frequency anyway, i.e.,
```
    freq = u.Quantity(in_x, u.AA, dtype=np.float64).to(u.Hz)
```