add function for midpoint integration (output matches c code to rtol ~ 1.2 * 10**-3)
add function for midpoint integration (output matches c code to rtol ~ 1.2 * 10**-3)

jesford/smd_nfw,jesford/smd-nfw