Bé, en Matemàtica Aplicada ja se sap que ens ve de gust utilitzar l'ordinador per a cercar aproximacions a les solucions dels problemes.

Podríem fer un primer programa en MatLab utilitzant "força bruta" per analitzar tots els conjunts de 6 punts en un quadrat on tant l'abcisa com l'ordenada varien d'1 a N. Però aquest programa seria notòriament ineficient a l'haver d'analitzar N¹² conjunts de punts.

Un enfocament més eficient seria fixar 2 punts amb distància mínima m=1, i tenint en compte que en un hexàgon regular la distància màxima és M=2·m=2, escollir successivament els altres 4 punts dins de la intersecció de corones circulars amb radis 1 i 2 centrades en els punts anteriors. D'aquesta manera, l'espai per a escollir els punts restants es va reduint progressivament.

Podem utilitzar un segon programa en MatLab per a fer això, utilitzant com a paràmetre el valor n dels increments d'abcises i ordinades que prenem dins de l'espai de variació per a cadascuna.

Amb punts inicials [2 , 2] i [3 , 2] i n=0'1, el mínim valor trobat de M/m és 1'943066 amb els punts

2.0000    2.0000     3.0000    2.0000     2.3000    1.0265     1.2000    2.6282     2.3000    2.9579     1.2000    1.3876

que podeu veure gràficament i que s'aproximen a un pentàgon regular més el seu centre. Si prenguérem exactament un pentàgon regular més el seu centre, tindríem M/m = 4·cos(0'2·π)·cos(0'3·π) ≈ 1'902 .

Naturalment, prenent valors menors dels increments n obtindríem una millor aproximació a la solució òptima amb valor mínim de M/m, que en tot cas no pot diferir de l'obtingut més de 4n√2 .