scopesim_templates.stellar.imf.prim_power# scopesim_templates.stellar.imf.prim_power(m, power)# Take floats or arrays, but return arrays. returns: m**(power + 1) / (power + 1) and handles the case when power = -1