Analytic functions with blackbody

See #3077 
Perhaps the function should be named `planck_blackbody` or `planck_bb`?  The user will typically see it in the `analytic_functions` namespace, so `planck` might be slightly unclear? (also not a terribly strong opinion if someone else really wants it to be short)