Pythagoräische Quadrupel
a² + b² + c² = d² mit ganzzahligen Werten.
a = 2u
b = 2v
c = (u² + v² - t²) / t
d = (u² + v² + t²) / t
wobei t ein echter Teiler von u²+v² mit t² < u²+v².
Für t wird ein zufälliger Teiler ausgewählt. Andere taugliche Werte selber ermitteln.