def imf(m):
m1 = <NUM_LIT>; m2 = <NUM_LIT><EOL>a1 = <NUM_LIT>; a2 = <NUM_LIT>; a3 = <NUM_LIT><EOL>const2 = m1**-a1 -m1**-a2 <EOL>const3 = m2**-a2 -m2**-a3 <EOL>if m < <NUM_LIT>:<EOL><INDENT>alpha = <NUM_LIT><EOL>const = -const2 -const3<EOL><DEDENT>elif m < <NUM_LIT>:<EOL><INDENT>alpha = <NUM_LIT><EOL>const = -const3<EOL><DEDENT>else:<EOL><INDENT>alpha = <NUM_LIT><EOL>const = <NUM_LIT:0.0><EOL><DEDENT>return m**-alpha + const<EOL>
Returns
-------
N(M)dM
    for given mass according to Kroupa IMF, vectorization
    available via vimf()