Python routines for numerical problems in stellar astrophysics