Extended Arrhenius equation in Python

Project dedicated to estimate prototype stability based on the extended Arrhenius equation - so called ASAP (Accelerated Stability Assessment Program).
Experience in: coding, statistical analysis, kinetics, stability as well mathematical models evaluation.

What I need:
Code written in Python

Main goal: estimation of time required to reach specification limit for extended Arrhenius equation (called ASAP - Accelerated Stability Assessment Program - isoconversion)

addition 1: Monte-Carlo simulation required to support probability of reaching specification limit (by CI at 95%)

addition 2: solve equation based on multiple conditions results (chemical test for impurities)

addition 3: statistical analysis of equation and "best" model fitting

addition 4: additional option to estimate influence of temperature excursion on equation result (eg. how 48h at 50C influence on the long term stability of sample at conditions 25C/60% RH)

addition 5: option to add the "B" parameter value for already estimated packaging material
Supporting materials: https://freethinktech.com/freethink-presentations/