Analytic functions with blackbody

See #3077 
I don't think it's necessary to do the full conversion here, because it is done anyway in `blackbody_nu()` later. But I can change `isinstance(in_x, u.Quantity)` to `hasattr(in_x, 'unit')`.