Analytic functions with blackbody

See #3077 
The following can be greatly simplified:
```
with u.add_enabled_equivalencies(u.spectral() + u.temperature() ):
    w = u.Quantity(in_x, u.AA, dtype=np.float64)
    temp = u.Quantity(temperature, u.K, dtype=np.float64)
```