This code was used to create Pourbaix diagrams from elemental lead and native lead oxides in aqueous solutions. It serves as a base for other thermodynamic phase diagrams throughout the work. Code was created and written by Alex T. Tai, with the help of Lauren N. Walters.