C++ zero temperature Chiral Mean Field implementation with spinodal solutions

The CMF model, which is used to describe neutron stars, neutron star mergers, and heavy-ion collisions, is a very successful model for the equation of state at large baryon densities. This proprietary was initially written in Fortran77 and relies on legacy routines. In this work, we report the updates of a C++ implementation of the CMF model for zero temperature within the MUSES collaboration. GNU GSL has been adopted for integration and multidimensional root-finding algorithms. In our code, metastable solutions are reachable; therefore, spinodal solutions of the field can be studied. Finally, this new code allows the computation of higher-order derivatives required for the nuclear matter hydrodynamical evolution codes.