This was done in Mathematica 11, using ParametricIPOPTMinimize
The objective function is the apothem of the hexagon, the shortest distance from the center to a face. Minimizing it is equivalent to minimizing the area of the hexagon.
Nicely fits together! Very impressive! What is optimized(minimized/maximed) exactly? Total area? or total circumference? or?