An embedded pair of cubature formulae of degree 5 and 7 for the triangle (Q1105317): Difference between revisions

Au début du travail les AA. font l'observation que pour construire un procédé automatique pour le calcul des intégrales dans l'espace à deux dimensions de la forme \(I[f]=\iint_{R}w(x,y)f(x,y)dxdy\) rélatives à une région R, deux formules de cubature \(Q_ j[f]=\sum^{N_ j}_{i=1}w_{ij}f(x_{ij},y_{ij}),\) \(j\in \{1,2\}\) sont nécessaires. L'expression \(| Q_ 1[f]-Q_ 2[f]|\) peut être utilisée dans la qualité d'évaluation de l'erreur pour des formules de cubature plus peu exactes. Pour réduire le nombre d'évaluations de la fonction on peut essayer de déterminer \(Q_ 2\) ainsi que \(\{(x_{i1},y_{i1})| \quad i=1,...,N_ 1\}\subset \{(x_{i2},y_{i2})| \quad i=1,...,N_ 2\}\) ou par \(Q_ 1\) on désigne un procédé connu. Le problème posé dans cette manière peut être considéré comme une extension de l'idée de Kronrod. Les AA. [Numer. Math. 49, 269-274 (1986; Zbl 0574.65014)] ont construit déja des paires de formules de cubature pour des régions symétriques (rectangulaires). Dans le travail de [ACM Math. Softw. 8, 210-218 (1982; Zbl 0478.65011)], on utilise aussi une paire de formules de cubature pour le triangle. La deuxième partie du travail est consacrée à la démonstration des deux théorèmes. Le théorème 1 est utilisé dans l'application donnée dans ce traivail. Dans la dernière partie de ce travail (application au triangle) on considère l'intégrale I où \(R=T_ 2=\{(x,y),\quad x\geq 0,\quad y\geq 0,\quad x+y\leq 1\}\) et \(w(x,y)=1\). On utilise la formule de Radon aux 7 points, ayant le degré 5 d'exactitude, dans la qualité de \(Q_ 1\). \textit{D. P. Laurie} [NRIMS Tech. Rep. TWISK 78 (1979; Zbl 0445.65014)] a donnée une formule de cubature \(Q_ 2\) du degré 8 d'exactitude sur 19 noeuds. Les AA. font l'observation qu'il n'y avait pas une formule de cubature \(Q_ 2\) du degré 7 d'exactitude. Ils donnent une formule de cubature \(Q_ 2\) du degré 7 d'exactitude pour le triangle \(T_ 2\), sur 16 noeuds. Pour celle-ci on utilise les résultats du théoreme 1 de la deuxième partie du travail. La formule de cubature est donnée dans le tableau 1.