title: Séries de Minggatu
text: # Séries de Minggatu
  
  En mathématiques, les séries de Minggatu sont des développements en série de fonctions trigonométriques découverts vers 1730 par Minggatu, un mathématicien mongol au service des empereurs de la dynastie Qing. Ses résultats, parmi lesquels figure la découverte des nombres de Catalan, s'inspirent des travaux des mathématiciens européens, mais sont obtenus par des méthodes géométriques et algébriques totalement indépendantes des leurs ; ils ont été publiés (après sa mort) dans son ouvrage majeur, le Geyuan Milu Jifa (Méthode rapide pour obtenir les divisions exactes d'un cercle).
  
  ## Contexte historique
  
  En 1701, le missionnaire jésuite Pierre Jartoux (1669-1720) vint en Chine, apportant avec lui trois développements en série de fonctions trigonométriques dus à Isaac Newton et James Gregory :
  
      $\pi =3\left(1+{\frac {1}{4\cdot 3!}}+{\frac {3^{2}}{4^{2}\cdot 5!}}+{\frac {3^{2}\cdot 5^{2}}{4^{3}\cdot 7!}}+\cdots \right)$
      $\sin x=x-{\frac {x^{3}}{3!}}+{\frac {x^{5}}{5!}}-{\frac {x^{7}}{7!}}+\cdots$
      $\operatorname {vers} x=1-\cos x={\frac {x^{2}}{2!}}-{\frac {x^{4}}{4!}}+{\frac {x^{6}}{6!}}+\cdots .$
  
  Ces séries suscitèrent un grand intérêt parmi les mathématiciens chinois, le calcul de $\pi$ par ces « méthodes rapides » ne demandant que des multiplications et des divisions, et s'avérant beaucoup plus efficace que les calculs à l'aide de racines carrées de l'algorithme classique de Liu Hui. Cependant, Jartoux n'expliqua pas la méthode employée pour construire ces séries, et Minggatu le soupçonna de ne pas vouloir en partager le secret. Il travailla par intermittence sur ce problème durant trente années, et obtint un modèle géométrique de construction de séries trigonométriques, qui lui permit non seulement de retrouver les trois séries précédentes, mais également six autres, découvrant au passage et utilisant les nombres de Catalan, au moins 20 ans avant les mathématiciens européens.
  
  ## Un exemple de calcul de séries : la méthode des deux segments
  
  La construction ci-contre (figure 2) illustre la méthode de Minggatu. L'arc BCD est une portion du cercle unité (de centre A) ; on pose BC = CD = x ; construisant des triangles isocèles BGC et DCH, on a BD = 2x-GH. On pose alors EJ = EF, FK = FJ et on prolonge BE jusqu'en L tel que EL = BE ; de même, on prolonge BF jusqu'à M avec BF = MF. L, M et C sont alignés. Posant BI = BC, on a les proportionnalités
  
      $AB:BC:CI=1:x:x^{2}$
  
  On obtient de même une série de triangles semblables ABE, BEF, FJK, BLM, CMO, MOP, CGH et de plus les triangles CMO et EFJ sont égaux, aussi
  
      $AB:BE:EF:FJ:JK=1:p:p^{2}:p^{3}:p^{4}$,
      $1:BE=BE:EF;$ donc $EF=BE^{2}$
      $1:BE^{2}=x:GH$, d'où
      $GH=x\cdot BE^{2}=xp^{2}$,
      et $BD=2x-xp^{2}$.
  
  ABEC et BLIN étant semblables,
  
      $EF=LC=CM=MG=NG=IN$
      $LM+MN=CM+MN+IN=CI+OP=JK+CI$
      $\therefore AB:(BE+EC)=BL:(LM+MN)$ et $AB:BL=BL:(CI+JK)$
      Posant $BL=q$
      $AB:BL:(CI+JK)=1:q:q^{2}$
      $JK=p^{4}$, $CI=y^{2}$ et $CI+JK=q^{2}=BL^{2}=(2BE)^{2}=(2p)^{2}=4p^{2}$.
  
  Finalement, des manipulations algébriques amènent à
  
      ${\frac {(x^{2})^{2}}{16}}={\frac {(q^{2}-{\frac {q^{4}}{16}})^{2}}{16}}=\sum _{j=0}^{2}(-1)^{j}{2 \choose j}{\frac {q^{2(2+j)}}{16^{j+1}}}$
  
  et plus généralement (en itérant la construction) :
  
      ${\frac {x^{2n}}{16^{n-1}}}=\sum _{j=0}^{n}(-1)^{j}{n \choose j}{\frac {q^{2(n+j)}}{16^{n+j-1}}}$.
  
  Éliminant les termes en $q^{n}$ sauf le premier et le dernier, on aboutit à
  
      ${\begin{aligned}&x^{2}+{\frac {x^{4}}{16}}+{\frac {2x^{6}}{16^{2}}}+{\frac {5x^{8}}{16^{3}}}+{\frac {14x^{10}}{16^{4}}}+{\frac {42x^{12}}{16^{5}}}\\[10pt]&+{\frac {132x^{14}}{16^{6}}}+{\frac {429x^{16}}{16^{7}}}+{\frac {1430x^{18}}{16^{8}}}+{\frac {4862x^{20}}{16^{9}}}\\[10pt]&+{\frac {16796x^{22}}{16^{10}}}+{\frac {58786x^{24}}{16^{11}}}+{\frac {208012x^{26}}{16^{12}}}\\[10pt]&+{\frac {742900x^{28}}{16^{13}}}+{\frac {2674440x^{30}}{16^{14}}}+{\frac {9694845x^{32}}{16^{15}}}\\[10pt]&+{\frac {35357670x^{34}}{16^{16}}}+{\frac {129644790x^{36}}{16^{17}}}\\[10pt]&+{\frac {477638700x^{38}}{16^{18}}}+{\frac {1767263190x^{40}}{16^{19}}}+{\frac {6564120420x^{42}}{16^{20}}}\\[10pt]&=q^{2}+{\frac {62985}{8796093022208}}q^{24}\end{aligned}}$
  
  Les coefficients des numérateurs : 1, 1, 2, 5, 14, 42, 132 ...... (voir la figure 1) sont les nombres de Catalan ; Minggatu est le premier à les avoir $décrits^,.$
  
  Ainsi :
  
      $q^{2}=\sum _{n=1}^{\infty }C_{n}{\frac {x^{2n}}{4^{2n-2}}}$$^,,$
  
  où $C_{n}={\frac {1}{n+1}}{2n \choose n}$ est un nombre de Catalan. Minggatu introduisit ainsi l'usage de relations de récurrence dans les mathématiques chinoises, montrant en particulier que $C_{n}=\sum _{k}(-1)^{k}{nk \choose k+1}C_{nk}$.
  
  Si, dans la figure 2, on pose l'angle ${\widehat {(AB,AE)}}=\alpha$, on a alors ${\widehat {(AB,AC)}}=2\alpha$ ; les relations précédentes permettent alors à Minggatu d'en déduire la relation :
  
      $\sin(2\alpha )=2\sin \alpha -\sum _{n=1}^{\infty }C_{n}{\frac {(\sin \alpha )^{2n+1}}{4^{n-1}}}$
  
  (cette dernière formule se démontre aisément à l'aide de la série génératrice découverte par Euler 20 ans plus tard : $\sum _{n=0}^{\infty }C_{n}x^{n+1}={\frac {1-{\sqrt {1-4x}}}{2}}$).
  
  ## La limite quand le nombre de segments augmente
  
  Par des calculs analogues, Minggatu obtient des séries correspondant à des décompositions en n= 3, 5, 10, 100, 1000 et 10000 segments ; il en déduit (de manière empirique) la série correspondant au cas limite quand n tend vers l'infini. Ainsi, il réécrit
  
  $y_{100}=100a-{\frac {(100a)^{3}}{24.002400240024002400}}+{\frac {(100a)^{5}}{24.024021859697730358*80}}+$..........
  
  $y_{1000}=1000a-{\frac {(1000a)^{3}}{24.000024000024000024}}+{\frac {(1000a)^{5}}{24.000240002184019680*80}}+$..............
  
  $y_{10000}=10000a-{\frac {(10000a)^{3}}{24.000000240000002400}}+{\frac {(10000a)^{5}}{24.000002400000218400*80}}+$..................
  
  (où les y sont les longueurs de la corde et les a les longueurs des segments subdivisant l'arc) et remarque que les dénominateurs 24.000000240000002400，24.000002400000218400*80 approchent 24 et 24*80 respectivement, tandis que n*a (100a, 1000a, 10000a) tend vers la longueur de l'arc ; $ainsi^,$
  
      $corde=arc-{\frac {arc^{3}}{4*3!}}+{\frac {arc^{5}}{4^{2}*5!}}-{\frac {arc^{7}}{4^{3}*7!}}+$.....
      $=\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}*arc^{2*n-1}}{(4^{n-1}*(2*n-1)!)}}$
  
  Par inversion de séries, il exprime enfin l'arc en fonction de la corde :
  
      $arc=corde+{\frac {corde^{3}}{24}}+{\frac {3*corde^{5}}{640}}+{\frac {5*corde^{7}}{7168}}+$............
  