adding lagrange multipliers
added lagrange constraint equations
now calculate derivatives of summations of nuclei number densities
adding sum of abundance fractions equal to one as a constraint
trying a new way to allocate memory

previous subroutine allocated unnecessary ‘temp’ tables
moved index definitions

make sure to keep all initial nuclei in the distribution even if they
do not undergo a reaction