In this talk I will discuss PRyMordial: a program dedicated to the computation of observables in the early universe with a focus on the cosmological era of Big Bang Nucleosynthesis (BBN). The code is the first of its kind written in python and offers fast and precise evaluation of both the BBN light-element abundances and the effective number of relativistic degrees of freedom. PRyMordial was created for both state-of-the-art analyses in the Standard Model as well as for general investigation of New Physics present in the early universe. In this talk, I will review the physics implemented in PRyMordial and provide a short guide on how to use the code for applications in the Standard Model and beyond.