Analytic functions with blackbody

See #3077 
OK, @pllim, no automatic units! And yes, let's go with SI here (maybe no default unit for the input in that case) and two wrappers with predefined flux unit and defaults for input of Angstrom and Hz.