def mimf_ferrario(mi):
mf=-<NUM_LIT>*mi**<NUM_LIT:6>+<NUM_LIT>*mi**<NUM_LIT:5>-<NUM_LIT>*mi**<NUM_LIT:4>+<NUM_LIT>*mi**<NUM_LIT:3>-<NUM_LIT>*mi**<NUM_LIT:2>+<NUM_LIT>*mi+<NUM_LIT><EOL>return mf<EOL>
Curvature MiMf from Ferrario etal. 2005MNRAS.361.1131.