added flux_envelope_3d()
- calculate the flux envelope for a triplet of reactions