FORMULE DE THOMAS SIMPSON

ARTICLE IV

DÉTERMINATION DES VALEURS NUMÉRIQUES DES SIX QUAN- TITÉS po^ ^j A', A", A', , A' (11° 39), ET DES QUATRE QUAN- TITÉS (7o, 0*2 , <y'', ^4 (n° 40), PAR LA FORMULE DE THOMAS SIMPSON.

49. Dans le cas des voûtes d'épaisseur variable, les nombres /5o) A, . . . , 0*0, 0*2, ... qui entrent dans les valeurs de Q et L ont été obtenus par des méthodes d'intégration indiquées dans les n®' 39 et 40. Les formules obtenues par ces méthodes exigent des calculs assez longs; aussi sera-t-il souvent plus simple et suffisamment exact de déterminer ces nombres par la méthode d'intégration approchée dite ybrmu/e de Thomas Simpson.

Rappelons brièvement en quoi consiste cette formule.

Elle a pour but de trouver la valeur numérique approchée d'une intégrale définie / jdx, dans laquelle j^ est la

dérivée d'une fonction de x que l'on ne sait pas déterminer sous forme finie. A cet eflfet, on calcule un certain nombre de valeurs de y correspondant à des valeurs de x croissant en progression arithmétique. On obtient ainsi les ordonnées d'un certain nombre de points de la courbe dont l'aire est égale à l'intégrale définie cherchée. En joignant ces points par des lignes droites, on obtient un polygone inscrit dont l'aire diffère d'autant moins de celle de la courbe que les points sont plus multipliés. Mais, afin d'éviter le calcul d'un trop grand nombre d'ordonnées , tout en conservant une approximation suffisante, on remplace les côtés du polygone par des arcs de paraboles du

lOO CHAPITRE III. — STABILITÉ DBS TOUTES EN ARC DE CERCLE

second degré passant par trois points consécutifs, et dont l'ordonnée , contenant trois paramètres indéterminés, est de la formey = ax- + bx -{- c.

Divisons la différence b — a des valeurs extrêmes de x en un nombre pair 2« de parties égales. Soient h Tune de ces parties, /o» J^4 7 J^a-. • • • > JK2/1 ^^s valeurs dej^ correspondant aux in-\-i valeurs de xa, a-j-A, a -h 2 A, . . ., i; Taire du polygone curviligne ainsi obtenu étant représentée par S y on aura

s=•3[ro+rî»^-4(rl^-ra-i-...-^J,„l)-^-2(rJ-f-r4-^•••-^-r«/I-2)l

dans laquelle la parenthèse est la somme des deux ordonnées extrêmes ayant pour indices zéro et 'in, plus quatre fois la somme des ordonnées d'indice impair, plus deux fois la somme des ordonnées d'indice pair.

On pourra prendre généralement la valeur de S pour celle de l'intégrale cherchée.

Lorsque la courbe présentera des parties d'une courbure prononcée, il pourra convenir de partager l'intervalle (i — a) en deux ou plusieurs parties, auxquelles on appliquera séparément la formule précédente en augmentant ou diminuant le nombre des divisions suivant que les valeurs de y varieront plus ou moins rapidement.

ARTICLE V