s4 - 4hs3 - (a2 - b2)2 = 0 |
Une fois cette équation résolue : a² + b² = u² + v² + 2d² = (u + v)² - 2uv + 2d² = (u + v)² - 2h(u + v) + 2d² donne la valeur de d :
d = √(a² + b² - s² + 2hs)/2 |
Le problème est soluble si 1/h = 1/u + 1/v > 1/a + 1/b, puisque bien évidemment u et v sont < a et b !
Application numérique :
Comme la résolution de l'équation du quatrième degré n'est pas quelque chose de sympathique,
on peut utiliser une méthode par approximations successives, comme la méthode de Newton :
En partant de s0 on calcule sn+1 = sn - f(sn)/f '(sn)
s0 = a + b
s = (3s4 - 8hs3 + (a2 - b2)2) / (4s3 - 12hs2) |
Avec a = 3, b = 2, h = 1, s0=5, s converge vers 4,311185197... et d=1,231186...
Un script pour effectuer ces calculs :
Choisissons alors arbitrairement deux triplets de Pythagore primitifs
x² + y² = z² et x'² + y'² = z'² et soit p le PGCD de x et x', en multipliant le premier triplet par x/p et le second par x'/p :
(xx'/p)² + (yx'/p)² = (zx'/p)² et (x'x/p)² + (y'x/p)² = (z'x/p)²
Posons xx'/p = d/m, yx'/p = u/m, zx'/p = a/m, y'x/p = v/m et z'x/p = b/m pour obtenir un système "d'équations du couloir" avec a, b, u, v, d entiers
h = uv/(u + v) = mxx'yy' / p(xy' + x'y). Soit alors q = PGCD(xx'yy'/p², (xy' + x'y)/p)
En choisissant m = (xy' + x'y)/pq, h est aussi un entier. Tout ceci donne donc des données et une solution entières pour un problème du couloir.
(r,s) entiers arbitraires, premiers entre eux et de parité opposée, de même pour (r',s')
x = r² - s², y = 2rs, z = r² + s² et de même pour x', y', z' p = PGCD(x,x'), q = PGCD(xx'yy'/p², (xy' + x'y)/p), m = (xy' + x'y)/pq a = mzx'/p, b = mz'x/p h = xx'yy'/qp² d = mxx'/p
|
Par exemple, r = 2, s = 1, r' = 3, s' = 2 donne x = 3, y = 4, z = 5, x' = 5, y' = 12, z' = 13
p = PGCD(3,5) = 1, q = PGCD(720,56) = 8, m = 7
a = 175, b = 273, h = 90 et d = 105 :
Deux échelles de 1,75m et 2,73m se croisent à 90cm du sol dans un couloir de 1,05m de large.
On peut aussi échanger x' et y' soit x = 3, y = 4, z = 5, x' = 12, y' = 5, z' = 13
p = PGCD(3,12) = 3, q = PGCD(80,21) = 1, m = 21
a = 420, b = 273, h = 80 et d = 252
La plus petite possibilité est : a = 70, b = 119, h = 30, d = 56 (u = 42, v = 105)
Obtenue à partir des triplets de Pythagore primitifs (x et y inversés) :
(4, 3, 5) r = 2, s = 1 et (8, 15, 17) r' = 4, s' = 1.
p = PGCD(8,4) = 4 donne d = 8m, b = 10m, u = 6m, a = 17m, v = 15m
donc h = uv/(u + v) = m×6×15/21 = 30m/7 entier si m = 7 (multiple de 7, mais on cherche la plus petite valeur).
Un script pour obtenir de telles valeurs rationnelles
Choisir les valeurs de r, s≠r et r', s'≠r'
Le programme en déduit les valeurs par les formules ci-dessus.