Creating variably-faced geodesic meshes through code

I'm currently trying to write something to generate a mesh for a Goldberg Polyhedron with any number of applicable faces (a mostly hexagonal sphere, and the number of hexagons can be theoretically infinitely large), but I'm really lost on how I would appropriately automate adding the additional hexagons through code. I can hammer out the initial 12 pentagons' vertex data manually, but building off these is where I'm stumped. Any ideas on how to approach this or any resources to learn how? I'm having trouble understanding the wikipedia page's math on how I'd translate that to code.