A dynamical computation of a multiple integral involved in the neutron star theory
This study consists in computing dynamically using the CADNA library an integral involved in the neutron star theory. Massive stars collapse with an explosion ejecting all the matter, except a central residue: a neutron star or a black hole, depending on the details of the scenario. Neutron stars are mainly done of neutrons with a typical proportion of protons smaller than 10%. With time, equilibrium between protons and neutrons takes place. Time to reach this equilibrium is obtained by calculating high order integrals involving several variables for any particle appearing in the reactions.
where and is a normalization function.
with and .
Each one-dimensional integral is evaluated by two iterative processes:
Using this strategy, the exact singificant digits of the approximation of
obtained are those of the exact value of the
integral, up to 6.
To get the FORTRAN source code.