Setup

import math
import numpy as np
import pandas as pd

hex_sides = [0.5]
for hex in range(2, 10):
    
    A = sum(hex_sides)
    
    a = 13
    b = 24*A
    c = 12*A**2 - 4
    d = b**2 - (4*a*c)

    #print(-1*b, d, 2*a)
    
    x = (-1*b + math.sqrt(d)) / (2*a)
    
    hex_sides.append(x)
    
hex_sides

[0.5,
 0.07692307692307693,
 0.0004271790970461572,
 1.3169502753275818e-08,
 0.0,
 0.0,
 0.0,
 0.0,
 0.0]

(-90 + math.sqrt(7*19*61)) / 13**2

0.00042717909704610463

Rohan Lewis

2022.01.23