for
𝐵
=
1
to whatever you wish
for
𝜇
=
0
to int (
𝐵
/
√
3
)
evaluate
𝐻
=
𝐵
2
+
𝜇
2
for
𝑘
=
2
𝜇
to int(
√
𝐻
)
if
𝑛
=
𝐻
/
𝑘
is an integer and gcd
(
𝐵
,
𝑘
,
𝑛
)
=
1
then the quintet of main curvatures is
(
−
𝐵
,
𝐵
+
𝑘
,
𝐵
+
𝑛
,
𝐵
+
𝑘
+
𝑛
−
2
𝜇
,
𝐵
+
𝑘
+
𝑛
+
2
𝜇
)
next
𝑘
next
𝜇
next
𝐵
Algorithm 1:
An algorithm that produces all integral disk packings.