Some practical quadrature formulas

Original title (in Romanian)

Câteva formule practice de cuadratura

Abstract

Authors

D.V. Ionescu
Institutul de Calcul

Dumitru V. Ionescu

Keywords

Quadrature formulas;

Paper coordinates

D.V. Ionescu, Some practical quadrature formulas. (Romanian) Com. Acad. R. P. Romîne 13 1963 689–695.

PDF

About this paper

Journal

Comunicările Academiei Republicii Populare Romîne

Publisher Name

?

DOI

–

Print ISSN

–

Online ISSN

–

??

Paper (preprint) in HTML form

1963-IonescuDV-Some-practical-formulas-of-quadrature-in-Communications-Acad.RPRvol.13nr.8-august1963pp html,body { width: 100%; height: 100%; } *, *::before,*::after { box-sizing: border-box; } body { margin: 0; color: #1E2029; font-size: 14px; line-height: normal; } hr { box-sizing: content-box; height: 0; overflow: visible; } h1, h2, h3, h4, h5, h6 { margin-top: 0; margin-bottom: 0.5em; color: rgba(0, 0, 0, 0.85); font-weight: 500; } p { margin-top: 0; margin-bottom: 1em; } ol, ul, dl { margin-top: 0; margin-bottom: 1em; } ol ol, ul ul, ol ul, ul ol { margin-bottom: 0; } dt { font-weight: 500; } dd { margin-bottom: 0.5em; margin-left: 0; } blockquote { margin: 0 0 1em; } dfn { font-style: italic; } b, strong { font-weight: bolder; } small { font-size: 80%; } sub, sup { position: relative; font-size: 75%; line-height: 0; vertical-align: baseline; } sub { bottom: -0.25em; } sup { top: -0.5em; } a { color: #0B93ff; text-decoration: none; background-color: transparent; outline: none; cursor: pointer; transition: color 0.3s; } a:hover { color: #33aaff; } a:active { color: #0070d9; } a:active, a:hover { text-decoration: none; outline: 0; } a[disabled] { color: rgba(0, 0, 0, 0.25); cursor: not-allowed; pointer-events: none; } pre, code, kbd, samp { font-size: 1em; } pre { margin-top: 0; margin-bottom: 1em; overflow: auto; } figure { margin: 0 0 1em; } img { vertical-align: middle; border-style: none; } svg:not(:root) { overflow: hidden; } table { border-collapse: collapse; } caption { padding-top: 0.75em; padding-bottom: 0.3em; color: rgba(0, 0, 0, 0.45); text-align: left; caption-side: bottom; } th { text-align: inherit; } mjx-container[jax="SVG"] { direction: ltr; } mjx-container[jax="SVG"] > svg { overflow: visible; min-height: 1px; min-width: 1px; } mjx-container[jax="SVG"] > svg a { fill: blue; stroke: blue; } mjx-assistive-mml { position: absolute !important; top: 0px; left: 0px; clip: rect(1px, 1px, 1px, 1px); padding: 1px 0px 0px 0px !important; border: 0px !important; display: block !important; width: auto !important; overflow: hidden !important; -webkit-touch-callout: none; -webkit-user-select: none; -khtml-user-select: none; -moz-user-select: none; -ms-user-select: none; user-select: none; } mjx-assistive-mml[display="block"] { width: 100% !important; } mjx-container[jax="SVG"][display="true"] { display: block; text-align: center; margin: 1em 0; } mjx-container[jax="SVG"][display="true"][width="full"] { display: flex; } mjx-container[jax="SVG"][justify="left"] { text-align: left; } mjx-container[jax="SVG"][justify="right"] { text-align: right; } g[data-mml-node="merror"] > g { fill: red; stroke: red; } g[data-mml-node="merror"] > rect[data-background] { fill: yellow; stroke: none; } g[data-mml-node="mtable"] > line[data-line], svg[data-table] > g > line[data-line] { stroke-width: 70px; fill: none; } g[data-mml-node="mtable"] > rect[data-frame], svg[data-table] > g > rect[data-frame] { stroke-width: 70px; fill: none; } g[data-mml-node="mtable"] > .mjx-dashed, svg[data-table] > g > .mjx-dashed { stroke-dasharray: 140; } g[data-mml-node="mtable"] > .mjx-dotted, svg[data-table] > g > .mjx-dotted { stroke-linecap: round; stroke-dasharray: 0,140; } g[data-mml-node="mtable"] > g > svg { overflow: visible; } [jax="SVG"] mjx-tool { display: inline-block; position: relative; width: 0; height: 0; } [jax="SVG"] mjx-tool > mjx-tip { position: absolute; top: 0; left: 0; } mjx-tool > mjx-tip { display: inline-block; padding: .2em; border: 1px solid #888; font-size: 70%; background-color: #F8F8F8; color: black; box-shadow: 2px 2px 5px #AAAAAA; } g[data-mml-node="maction"][data-toggle] { cursor: pointer; } mjx-status { display: block; position: fixed; left: 1em; bottom: 1em; min-width: 25%; padding: .2em .4em; border: 1px solid #888; font-size: 90%; background-color: #F8F8F8; color: black; } foreignObject[data-mjx-xml] { font-family: initial; line-height: normal; overflow: visible; } mjx-container[jax="SVG"] path[data-c], mjx-container[jax="SVG"] use[data-c] { stroke-width: 3; } #setText > div { justify-content: inherit; margin-top: 0; margin-bottom: 1em; } #setText div:last-child { margin-bottom: 0 !important; } #setText > br, #preview-content br { line-height: 1.2; } #preview-content > div { margin-top: 0; margin-bottom: 1em; } .proof > div, #preview-content .proof > div, #setText .proof > div, .theorem > div, #preview-content .theorem > div, #setText .theorem > div { margin-top: 1rem; } mjx-container { text-indent: 0; overflow-y: hidden; overflow-x: auto; padding-top: 1px; padding-bottom: 1px; } .math-inline mjx-container { display: inline-block !important; page-break-inside: avoid; max-width: 100%; padding: 0; line-height: 0; } .math-inline[data-overflow="visible"] mjx-container { overflow: visible; } .math-inline mjx-container mjx-assistive-mml { max-width: 100%; } .math-block { align-items: center; page-break-after: auto; page-break-inside: avoid; margin: 0; display: block; } .math-inline { display: inline-flex; max-width: 100%; } .math-block[data-width="full"] { overflow-x: auto; display: flex; } .math-block[data-width="full"] mjx-container[jax=SVG][display=true] { display: flex; flex: 1 1; justify-content: center; } svg .math-inline { display: initial; max-width: initial; } svg .math-inline mjx-container { max-width: initial; } svg mjx-container { overflow: visible; padding: 0; } svg .math-block[data-width="full"] { overflow: visible; } .math-block, .math-inline { --mmd-highlight-color: rgba(0, 147, 255, 0.25); --mmd-highlight-text-color: #1e2029; } .math-block[data-highlight-color] mjx-container[jax="SVG"] > svg { background-color: var(--mmd-highlight-color); } .math-block[data-highlight-text-color] mjx-container[jax="SVG"] > svg { color: var(--mmd-highlight-text-color); } .math-inline[data-highlight-color] mjx-container[jax="SVG"] { background-color: var(--mmd-highlight-color); } .math-inline[data-highlight-text-color] mjx-container[jax="SVG"] { color: var(--mmd-highlight-text-color); } .math-block p { flex-shrink: 1; } .math-block mjx-container { margin: 0 !important; } .math-error { background-color: yellow; color: red; } #preview-content img, #setText img { max-width: 100%; } #preview-content blockquote, #setText blockquote { page-break-inside: avoid; color: #666; margin: 0 0 1em 0; padding-left: 3em; border-left: .5em solid #eee; } #preview-content pre, #setText pre { border: none; padding: 0; overflow: auto; font-size: 85%; line-height: 1.45; border-radius: 6px; box-sizing: border-box; background: #f8f8fa; } #preview-content pre code, #setText pre code { padding: 1rem; display: block; overflow-x: auto; line-height: 24px; } #preview-content table, #setText table { display: table; overflow: auto; max-width: 100%; border-collapse: collapse; page-break-inside: avoid; margin-bottom: 1em; } #preview-content table th, #setText table th { text-align: center; font-weight: bold; } #preview-content table td, #setText table td, #preview-content table th, #setText table th { border: 1px solid #dfe2e5; padding: 6px 13px; } #preview-content table tr, #setText table tr { background-color: #fff; border-top: 1px solid #c6cbd1; } #preview-content table tr:nth-child(2n), #setText table tr:nth-child(2n) { background-color: #f6f8fa; } .main-title, #preview-content .main-title, #setText .main-title { text-align: center; line-height: 1.2; margin: 0 auto 1em auto; } .author, #preview-content .author, #setText .author { text-align: center; margin: 0 auto; display: flex; justify-content: center; flex-wrap: wrap; } .author p, #preview-content .author p, #setText .author p { min-width: 30%; max-width: 50%; padding: 0 7px; } .author > p > span, #preview-content .author > p > span, #setText .author > p > span { display: block; text-align: center; } .section-title, #preview-content .section-title, #setText .section-title { margin-top: 1.5em; } .abstract, #preview-content .abstract, #setText .abstract { text-align: justify; margin-bottom: 1em; } .abstract p, #preview-content .abstract p, #setText .abstract p { margin-bottom: 0; } #preview-content sup, #setText sup { top: -.5em; position: relative; font-size: 75%; line-height: 0; vertical-align: baseline; } .text-url, #preview-content .text-url, #setText .text-url { color: #0B93ff; cursor: text; pointer-events: none; } .text-url a:hover, #preview-content .text-url a:hover, #setText .text-url a:hover { color: #0B93ff; } mark, #preview-content mark, #setText mark { background-color: #feffe6; } span[data-underline-type] mark, #preview-content span[data-underline-type] mark, #setText span[data-underline-type] mark { background: inherit; background-color: #feffe6; padding-top: 0; padding-bottom: 0; } *[data-has-dotfill] { position: relative; overflow: hidden; } *[data-has-dotfill] .dotfill::after { position: absolute; padding-left: .25ch; content: " . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . " ". . . . . . . . . . . . . . . . . . . . . . . "; text-align: right; } .smiles, #preview-content .smiles, #setText .smiles { text-align: center; } div.svg-container, #preview-content > div.svg-container, #setText > div.svg-container { display: flex; justify-content: center; } @media print { #preview { font-size: 10pt !important; } .math-block svg, .math-inline svg { shape-rendering: crispEdges; margin-top: 1px; } #preview-content img, #setText img { display: block; } #preview-content .figure_img img, #setText .figure_img img { display: inline; } #preview-content h1, #preview-content h2, #preview-content h3, #preview-content h4, #preview-content h5, #preview-content h6, #setText h1, #setText h2, #setText h3, #setText h4, #setText h5, #setText h6 { page-break-inside: avoid; position: relative; border: 2px solid transparent; } #preview-content h1::after, #preview-content h2::after, #preview-content h3::after, #preview-content h4::after, #preview-content h5::after, #preview-content h6::after, #setText h1::after, #setText h2::after, #setText h3::after, #setText h4::after, #setText h5::after, #setText h6::after { content: ""; display: block; height: 100px; margin-bottom: -100px; position: relative; } } #preview-content code, #setText code { font-family: Inconsolata; font-size: inherit; display: initial; background: #f8f8fa; } #preview-content .mmd-highlight code, #setText .mmd-highlight code, #preview-content pre.mmd-highlight code, #setText pre.mmd-highlight code { background-color: transparent; } #preview-content pre code, #setText pre code { font-family: 'DM Mono', Inconsolata, monospace; color: #333; font-size: 15px; } .hljs-comment, .hljs-quote, #preview-content .hljs-comment, #setText .hljs-comment, #preview-content .hljs-quote, #setText .hljs-quote { color: #998; font-style: italic; } .hljs-command, #preview-content .hljs-command, #setText .hljs-command { color: #005cc5; } .hljs-keyword, .hljs-selector-tag, .hljs-subst, #preview-content .hljs-keyword, #setText .hljs-keyword, #preview-content .hljs-selector-tag, #setText .hljs-selector-tag, #preview-content .hljs-subst, #setText .hljs-subst { color: #d73a49; font-weight: bold; } .hljs-number, .hljs-literal, .hljs-variable, .hljs-template-variable, .hljs-tag .hljs-attr, #preview-content .hljs-number, #setText .hljs-number, #preview-content .hljs-literal, #setText .hljs-literal, #preview-content .hljs-variable, #setText .hljs-variable, #preview-content .hljs-template-variable, #setText .hljs-template-variable, #preview-content .hljs-tag .hljs-attr, #setText .hljs-tag .hljs-attr { color: #005cc5; } .hljs-string, .hljs-doctag, #preview-content .hljs-string, #setText .hljs-string, #preview-content .hljs-doctag, #setText .hljs-doctag { color: #24292e; } .hljs-title, .hljs-section, .hljs-selector-id, #preview-content .hljs-title, #setText .hljs-title, #preview-content .hljs-section, #setText .hljs-section, #preview-content .hljs-selector-id, #setText .hljs-selector-id { color: #6f42c1; font-weight: bold; } .hljs-subst, #preview-content .hljs-subst, #setText .hljs-subst { font-weight: normal; } .hljs-type, .hljs-class .hljs-title, #preview-content .hljs-type, #setText .hljs-type, #preview-content .hljs-class .hljs-title, #setText .hljs-class .hljs-title { color: #458; font-weight: bold; } .hljs-tag, .hljs-name, .hljs-attribute, #preview-content .hljs-tag, #setText .hljs-tag, #preview-content .hljs-name, #setText .hljs-name, #preview-content .hljs-attribute, #setText .hljs-attribute { color: #000080; font-weight: normal; } .hljs-regexp, .hljs-link, #preview-content .hljs-regexp, #setText .hljs-regexp, #preview-content .hljs-link, #setText .hljs-link { color: #009926; } .hljs-symbol, .hljs-bullet, #preview-content .hljs-symbol, #setText .hljs-symbol, #preview-content .hljs-bullet, #setText .hljs-bullet { color: #990073; } .hljs-built_in, .hljs-builtin-name, #preview-content .hljs-built_in, #setText .hljs-built_in, #preview-content .hljs-builtin-name, #setText .hljs-builtin-name { color: #24292e; } .hljs-meta, #preview-content .hljs-meta, #setText .hljs-meta { color: #999; font-weight: bold; } .hljs-meta-keyword, #preview-content .hljs-meta-keyword, #setText .hljs-meta-keyword { color: #d73a49; } .hljs-meta-string, #preview-content .hljs-meta-string, #setText .hljs-meta-string { color: #032f62; } .hljs-deletion, #preview-content .hljs-deletion, #setText .hljs-deletion { background: #fdd; } .hljs-addition, #preview-content .hljs-addition, #setText .hljs-addition { background: #dfd; } .hljs-emphasis, #preview-content .hljs-emphasis, #setText .hljs-emphasis { font-style: italic; } .hljs-strong, #preview-content .hljs-strong, #setText .hljs-strong { font-weight: bold; } .table_tabular, #preview-content .table_tabular, #setText .table_tabular { overflow-x: auto; padding: 0 2px 0.5em 2px; } .tabular, #preview-content .tabular, #setText .tabular { display: inline-table !important; width: auto; table-layout: auto; border-collapse: collapse; border-spacing: 0; margin: 0 0 1em; font-size: inherit; height: fit-content; } .tabular th, #preview-content .tabular th, #setText .tabular th { border: none !important; padding: 6px 13px; background-color: transparent; font-weight: bold; } .tabular tr, #preview-content .tabular tr, #setText .tabular tr { border-top: none !important; border-bottom: none !important; background-color: transparent; } .tabular td, #preview-content .tabular td, #setText .tabular td { border-style: none !important; background-color: #fff; border-color: #000 !important; word-break: keep-all; padding: 0.1em 0.5em !important; } .tabular td > p, #preview-content .tabular td > p, #setText .tabular td > p { margin-bottom: 0; margin-top: 0; } .tabular td._empty, #preview-content .tabular td._empty, #setText .tabular td._empty { height: 1.3em; } .tabular td .f, #preview-content .tabular td .f, #setText .tabular td .f { opacity: 0; } #preview-content .sub-table table, #setText .sub-table table { margin-bottom: 0; } html[data-theme="dark"] .tabular tr, html[data-theme="dark"] #preview-content .tabular tr, html[data-theme="dark"] #setText .tabular tr, html[data-theme="dark"] .tabular td, html[data-theme="dark"] #preview-content .tabular td, html[data-theme="dark"] #setText .tabular td { background-color: #202226; border-color: #fff !important; } .figure_img, #preview-content .figure_img, #setText .figure_img { margin-bottom: 0.5em; overflow-x: auto; } div.figure_img img, #preview-content div.figure_img img, #setText div.figure_img img { display: inline; margin: 0; } ol.enumerate, ul.itemize, #preview-content ol.enumerate, #setText ol.enumerate, #preview-content ul.itemize, #setText ul.itemize { padding-inline-start: 40px; margin: 0 0 1em 0; } li > ol.enumerate, li > ul.itemize, #preview-content li > ol.enumerate, #setText li > ol.enumerate, #preview-content li > ul.itemize, #setText li > ul.itemize { margin: 0; } ul.itemize > li, #preview-content ul.itemize > li, #setText ul.itemize > li { position: relative; min-height: 1.4em; margin-bottom: 0; } .enumerate > .li_enumerate, #preview-content .enumerate > .li_enumerate, #setText .enumerate > .li_enumerate { margin-bottom: 0; } .itemize > li > span.li_level, .li_enumerate.not_number > span.li_level, #preview-content .itemize > li > span.li_level, #setText .itemize > li > span.li_level, #preview-content .li_enumerate.not_number > span.li_level, #setText .li_enumerate.not_number > span.li_level { position: absolute; right: 100%; white-space: nowrap; width: max-content; display: flex; justify-content: flex-end; padding-right: 10px; box-sizing: border-box; } .li_enumerate.not_number, #preview-content .li_enumerate.not_number, #setText .li_enumerate.not_number { position: relative; display: inline-block; list-style-type: none; min-height: 1.4em; } .itemize > li > span.li_level .math-inline, .enumerate > li > span.li_level .math-inline, #preview-content .itemize > li > span.li_level .math-inline, #setText .itemize > li > span.li_level .math-inline, #preview-content .enumerate > li > span.li_level .math-inline, #setText .enumerate > li > span.li_level .math-inline { display: inline-block; } #preview { font-family: 'CMU Serif', 'Georgia', Helvetica, Arial, sans-serif; font-size: 17px; visibility: visible; word-break: break-word; padding: 2.5em; max-width: 800px; margin: auto; box-sizing: content-box; } #preview h1, #preview h2, #preview h3, #preview h4, #preview h5, #preview strong { font-family: 'CMU Serif Bold', 'Georgia', Helvetica, Arial, sans-serif; } #preview i, #preview em { font-family: 'CMU Serif Italic', 'Georgia', Helvetica, Arial, sans-serif; } .mmd-menu { max-width: 320px; position: absolute; background-color: white; color: black; width: auto; padding: 5px 0; border: 1px solid #E5E6EB; margin: 0; cursor: default; font: menu; text-align: left; text-indent: 0; text-transform: none; line-height: normal; letter-spacing: normal; word-spacing: normal; word-wrap: normal; white-space: nowrap; float: none; z-index: 201; border-radius: 5px; -webkit-border-radius: 5px; -moz-border-radius: 5px; -khtml-border-radius: 5px; box-shadow: 0 10px 20px #808080; -webkit-box-shadow: 0 10px 20px #808080; -moz-box-shadow: 0 10px 20px #808080; -khtml-box-shadow: 0 10px 20px #808080; } .mmd-menu:focus { outline: none; } .mmd-menu.mmd-menu-sm { max-width: 100vw; padding-bottom: 34px; border-radius: 0; -webkit-border-radius: 0; -moz-border-radius: 0; -khtml-border-radius: 0; } .mmd-menu-item-icon { color: #1E2029; margin-left: auto; align-items: center; display: none; flex-shrink: 0; } .mmd-menu-item { padding: 8px 1.25rem; display: flex; background: transparent; height: 52px; max-height: 52px; } .mmd-menu-item:focus { outline: none; } .mmd-menu-item.active { background-color: #e1e0e5; } .mmd-menu-item.active .mmd-menu-item-icon { display: flex; } .mmd-menu-item-container { overflow: hidden; } .mmd-menu-item-title { color: #1E2029; white-space: nowrap; text-overflow: ellipsis; overflow: hidden; font-size: 14px; line-height: 20px; } .mmd-menu-item-value { color: #7d829c; white-space: nowrap; text-overflow: ellipsis; overflow: hidden; font-size: 12px; line-height: 16px; } html[data-theme="dark"] .mmd-menu-item-title { color: #ebefe7; } html[data-theme="dark"] .mmd-menu-item.active .mmd-menu-item-title { color: #1E2029; } html[data-theme="dark"] .mmd-menu { background-color: #33363a; } .mmd-context-menu-overlay { background: rgba(0, 0, 0, 0.56); } .ClipboardButton { padding: 0; margin: 0.5rem; display: inline-block; cursor: pointer; color: rgb(36, 41, 47); background: rgb(246, 248, 250); border-radius: 6px; border: 1px solid rgba(31, 35, 40, 0.15); box-shadow: rgba(31, 35, 40, 0.04) 0 1px 0 0, rgba(255, 255, 255, 0.25) 0 1px 0 0 inset; position: relative; } .ClipboardButton:hover { background-color: rgb(243, 244, 246); border-color: rgba(31, 35, 40, 0.15); transition-duration: .1s; } td .mmd-clipboard-copy-container { line-height: 0; } td .ClipboardButton { margin: 1px; line-height: 0; } td .ClipboardButton .mmd-clipboard-icon { margin: 1px !important; transform: scale(0.6); } .mmd-clipboard-icon { fill: currentColor; vertical-align: text-bottom; } .mmd-clipboard-copy-icon { color: rgb(101, 109, 118); } .mmd-clipboard-check-icon { color: rgb(26, 127, 55); } .mmd-tooltipped-no-delay:hover::before, .mmd-tooltipped-no-delay:hover::after { animation-delay: 0s; } .mmd-tooltipped:hover::before, .mmd-tooltipped:hover::after { display: inline-block; text-decoration: none; animation-name: tooltip-appear; animation-duration: .1s; animation-fill-mode: forwards; animation-timing-function: ease-in; animation-delay: .4s; } .mmd-tooltipped-w::before { top: 50%; bottom: 50%; left: -7px; margin-top: -6px; border-left-color: rgb(36, 41, 47); } .mmd-tooltipped::before { position: absolute; z-index: 1000001; display: none; width: 0; height: 0; color: rgb(36, 41, 47); pointer-events: none; content: ""; border: 6px solid transparent; opacity: 0; } .mmd-tooltipped-w::after { right: 100%; bottom: 50%; margin-right: 6px; transform: translateY(50%); } .mmd-tooltipped::after { position: absolute; z-index: 1000000; display: none; padding: 0.5em 0.75em; font: normal normal 11px/1.5 'CMU Serif', 'Georgia', Helvetica, Arial, sans-serif; -webkit-font-smoothing: subpixel-antialiased; color: white; text-align: center; text-decoration: none; text-shadow: none; text-transform: none; letter-spacing: normal; word-wrap: break-word; white-space: pre; pointer-events: none; content: attr(aria-label); background: rgb(36, 41, 47); border-radius: 6px; opacity: 0; }

COMMUNICATIONS
OF THE ACADEMY OF THE ROMANIAN PEOPLE'S REPUBLIC

Volume XIII, No. 8
August 1963
SOME PRACTICAL QUADRATURE FORMULAS
BY
D. V. IONESCU
Communication presented by academician T. Popoviciu in the journal of May 18, 1963
G. Coulmy [1] considered the following quadrature formula
∫ x n x n f ( x ) d x = 13 h 36 [ f ( x 0 ) + f ( x n ) ] + 7 h 6 [ f ( x 1 ) + f ( x n − 1 ) ] + h [ f ( x n ) + f ( x n − 2 ) ] + (C) + 35 h 36 [ f ( x 3 ) + f ( x n − 3 ) ] + h [ f ( x 4 ) + f ( x 5 ) + ⋯ + f ( x n − 4 ) ] ∫ x n x n   f ( x ) d x = 13 h 36 f x 0 + f x n + 7 h 6 f x 1 + f x n − 1 + h f x n + f x n − 2 + (C) + 35 h 36 f x 3 + f x n − 3 + h f x 4 + f x 5 + ⋯ + f x n − 4 {:[int_(x_(n))^(x_(n))f(x)dx=(13 h)/(36)[f(x_(0))+f(x_(n))]+(7h)/(6)[f(x_(1))+f(x_(n-1))]+h[f(x_(n))+f(x_(n-2))]+],[(C)quad+(35 h)/(36)[f(x_(3))+f(x_(n-3))]+h[f(x_(4))+f(x_(5))+cdots+f(x_(n-4))]]:}\begin{align*} & \int_{x_{n}}^{x_{n}} f(x) \mathrm{d} x=\frac{13 h}{36}\left[f\left(x_{0}\right)+f\left(x_{n}\right)\right]+\frac{7 h}{6}\left[f\left(x_{1}\right)+f\left(x_{n-1}\right)\right]+h\left[f\left(x_{n}\right)+f\left(x_{n-2}\right)\right]+ \\ & \quad+\frac{35 h}{36}\left[f\left(x_{3}\right)+f\left(x_{n-3}\right)\right]+h\left[f\left(x_{4}\right)+f\left(x_{5}\right)+\cdots+f\left(x_{n-4}\right)\right] \tag{C} \end{align*}∫xnxnf(x)dx=13h36[f(x0)+f(xn)]+7h6[f(x1)+f(xn−1)]+h[f(xn)+f(xn−2)]+(C)+35h36[f(x3)+f(xn−3)]+h[f(x4)+f(x5)+⋯+f(xn−4)]
where the nodes x 0 , x 1 , … , x n x 0 , x 1 , … , x n x_(0),x_(1),dots,x_(n)x_{0}, x_{1}, \ldots, x_{n}x0,x1,…,xnare in arithmetic progression with the ratio h h hhh, and applied this formula to the problem of tides.
It is good to compare formula (C) with the following well-known escapement formulas:
(T) ∫ x 0 x n f ( x ) d x = h 2 [ f ( x a ) + f ( x n ) ] + h [ f ( x 1 ) + f ( x 2 ) + … + f ( x n − 1 ) ] ∫ x 1 x n f ( x ) d x = 5 h 12 f ( x 0 ) + f ( x n ) ] + 13 h 12 [ f ( x 1 ) + f ( x n − 1 ) + h [ f ( x 2 ) + f ( x 3 ) + … (D) … + f ( x n − 2 ) ] ∫ x 0 x n f ( x ) d x = 3 h 8 [ f ( x 0 ) + f ( x n ) ] + 7 h 6 [ f ( x 1 ) + f ( x n − 1 ) ] + (L) + 23 h 24 [ f ( x 2 ) + f ( x n − 2 ) ] + h [ f ( x 3 ) + f ( x 4 ) + … + f ( x n − 3 ) ] (T) ∫ x 0 x n   f ( x ) d x = h 2 f x a + f x n + h f x 1 + f x 2 + … + f x n − 1 ∫ x 1 x n   f ( x ) d x = 5 h 12 f x 0 + f x n + 13 h 12 f x 1 + f x n − 1 + h f x 2 + f x 3 + … (D) … + f x n − 2 ∫ x 0 x n   f ( x ) d x = 3 h 8 f x 0 + f x n + 7 h 6 f x 1 + f x n − 1 + (L) + 23 h 24 f x 2 + f x n − 2 + h f x 3 + f x 4 + … + f x n − 3 {:[(T)int_(x_(0))^(x_(n))f(x)dx=(h)/(2)[f(x_(a))+f(x_(n))]+h[f(x_(1))+f(x_(2))+dots+f(x_(n-1))]],[{:int_(x_(1))^(x_(n))f(x)dx=(5h)/(12)f(x_(0))+f(x_(n))]+(13 h)/(12)[f(x_(1))+f(x_(n-1))+h[f(x_(2))+f(x_(3))+dots:}],[(D){: dots+f(x_(n-2))]],[int_(x_(0))^(x_(n))f(x)dx=(3h)/(8)[f(x_(0))+f(x_(n))]+(7h)/(6)[f(x_(1))+f(x_(n-1))]+],[(L)quad+(23 h)/(24)[f(x_(2))+f(x_(n-2))]+h[f(x_(3))+f(x_(4))+dots+f(x_(n-3))]]:}\begin{gather*} \int_{x_{0}}^{x_{n}} f(x) \mathrm{d} x=\frac{h}{2}\left[f\left(x_{a}\right)+f\left(x_{n}\right)\right]+h\left[f\left(x_{1}\right)+f\left(x_{2}\right)+\ldots+f\left(x_{n-1}\right)\right] \tag{T}\\ \left.\int_{x_{1}}^{x_{n}} f(x) \mathrm{d} x=\frac{5 h}{12} f\left(x_{0}\right)+f\left(x_{n}\right)\right]+\frac{13 h}{12}\left[f\left(x_{1}\right)+f\left(x_{n-1}\right)+h\left[f\left(x_{2}\right)+f\left(x_{3}\right)+\ldots\right.\right. \\ \left.\ldots+f\left(x_{n-2}\right)\right] \tag{D}\\ \int_{x_{0}}^{x_{n}} f(x) \mathrm{d} x=\frac{3 h}{8}\left[f\left(x_{0}\right)+f\left(x_{n}\right)\right]+\frac{7 h}{6}\left[f\left(x_{1}\right)+f\left(x_{n-1}\right)\right]+ \\ \quad+\frac{23 h}{24}\left[f\left(x_{2}\right)+f\left(x_{n-2}\right)\right]+h\left[f\left(x_{3}\right)+f\left(x_{4}\right)+\ldots+f\left(x_{n-3}\right)\right] \tag{L} \end{gather*}(T)∫x0xnf(x)dx=h2[f(xa)+f(xn)]+h[f(x1)+f(x2)+…+f(xn−1)]∫x1xnf(x)dx=5h12f(x0)+f(xn)]+13h12[f(x1)+f(xn−1)+h[f(x2)+f(x3)+…(D)…+f(xn−2)]∫x0xnf(x)dx=3h8[f(x0)+f(xn)]+7h6[f(x1)+f(xn−1)]+(L)+23h24[f(x2)+f(xn−2)]+h[f(x3)+f(x4)+…+f(xn−3)]
Formula (T) is derived from the classical trapezoid formula. Formulas (D) and (L) are cited in the book by H. Mineur [2] and are called the Durand and Lacroix formulas.
In this paper we will give the expressions of the ratio in formulas (C), (D) and (L) from which we will draw conclusions regarding the application of one or the other in practice.
I. THE REMAINDER IN FORMULA (C)
To find its remainder, it is assumed that the function f ( x ) f ( x ) f(x)f(x)f(x)it is classy C 9 [ x n , x n ] C 9 x n , x n C^(9)[x_(n),x_(n)]C^{9}\left[x_{n}, x_{n}\right]C9[xn,xn]and it is a method used consistently in the work [3]. It is found
(1) R c = ∫ e c e ∗ φ ( x ) f ′ ′ ( x ) d x , (1) R c = ∫ e c e ∗   φ ( x ) f ′ ′ ( x ) d x , {:(1)R_(c)=int_(e_(c))^(e^(**))varphi(x)f^('')(x)dx",":}\begin{equation*} R_{c}=\int_{e_{c}}^{e^{*}} \varphi(x) f^{\prime \prime}(x) \mathrm{d} x, \tag{1} \end{equation*}(1)Rc=∫andcand∗f(x)f′′(x)dx,
where the function φ ( x ) φ ( x ) varphi(x)\varphi(x)f(x)coincides on the intervals [ x 6 , x 1 ] , [ x 1 , x 2 ] , … x 6 , x 1 , x 1 , x 2 , … [x_(6),x_(1)],[x_(1),x_(2)],dots\left[x_{6}, x_{1}\right],\left[x_{1}, x_{2}\right], \ldots[x6,x1],[x1,x2],…with some natural polynomials of the second degree, which we write only for the first half of the interval [ x 0 , x 0 ] x 0 , x 0 [x_(0),x_(0)]\left[x_{0}, x_{0}\right][x0,x0] :
φ 1 ( x ) = ( x − x 0 ) 2 2 − 13 h 36 ( x − x 0 ) , φ 1 ( x ) = ( x − x 1 ) Ω 2 − 19 h 36 ( x − x 1 ) + 5 h 2 36 φ 1 ( x ) = x − x 0 2 2 − 13 h 36 x − x 0 , φ 1 ( x ) = x − x 1 Ω 2 − 19 h 36 x − x 1 + 5 h 2 36 varphi_(1)(x)=((x-x_(0))^(2))/(2)-(13 h)/(36)(x-x_(0)),quadvarphi_(1)(x)=((x-x_(1))^(Omega))/(2)-(19 h)/(36)(x-x_(1))+(5h^(2))/(36)\varphi_{1}(x)=\frac{\left(x-x_{0}\right)^{2}}{2}-\frac{13 h}{36}\left(x-x_{0}\right), \quad \varphi_{1}(x)=\frac{\left(x-x_{1}\right)^{\Omega}}{2}-\frac{19 h}{36}\left(x-x_{1}\right)+\frac{5 h^{2}}{36}f1(x)=(x−x0)22−13h36(x−x0),f1(x)=(x−x1)Oh2−19h36(x−x1)+5h236,
φ 2 ( x ) = ( x − x 2 ) ′ 2 − 19 h 36 ( x − x 2 ) + 4 h ′ 36 , φ k ( x ) = ( x − x k − 1 ) ′ 2 − h 2 ( x − x k − 1 ) + + h 2 12 ( k = 4 , 5 , … ) φ 2 ( x ) = x − x 2 ′ 2 − 19 h 36 x − x 2 + 4 h ′ 36 , φ k ( x ) = x − x k − 1 ′ 2 − h 2 x − x k − 1 + + h 2 12 ( k = 4 , 5 , … ) varphi_(2)(x)=((x-x_(2))^('))/(2)-(19 h)/(36)(x-x_(2))+(4h^('))/(36),quadvarphi_(k)(x)=((x-x_(k-1))^('))/(2)-(h)/(2)(x-x_(k-1))++(h^(2))/(12)(k=4,5,dots)\varphi_{2}(x)=\frac{\left(x-x_{2}\right)^{\prime}}{2}-\frac{19 h}{36}\left(x-x_{2}\right)+\frac{4 h^{\prime}}{36}, \quad \varphi_{k}(x)=\frac{\left(x-x_{k-1}\right)^{\prime}}{2}-\frac{h}{2}\left(x-x_{k-1}\right)+ +\frac{h^{2}}{12}(k=4,5, \ldots)f2(x)=(x−x2)′2−19h36(x−x2)+4h′36,fk(x)=(x−xk−1)′2−h2(x−xk−1)++h212(k=4,5,…).
If you find cå
∫ x 1 α 1 φ 1 ( x ) d x = − h 3 72 , ∫ x 1 x 1 φ 3 ( x ) d x = 3 h 3 72 , ∫ x 2 y 2 φ 3 ( x ) d x = h 8 72 , ∫ x 1 α 1   φ 1 ( x ) d x = − h 3 72 , ∫ x 1 x 1   φ 3 ( x ) d x = 3 h 3 72 , ∫ x 2 y 2   φ 3 ( x ) d x = h 8 72 , int_(x_(1))^(alpha_(1))varphi_(1)(x)dx=-(h^(3))/(72),quadint_(x_(1))^(x_(1))varphi_(3)(x)dx=(3h^(3))/(72),quadint_(x_(2))^(y_(2))varphi_(3)(x)dx=(h^(8))/(72),\int_{x_{1}}^{\alpha_{1}} \varphi_{1}(x) \mathrm{d} x=-\frac{h^{3}}{72}, \quad \int_{x_{1}}^{x_{1}} \varphi_{3}(x) \mathrm{d} x=\frac{3 h^{3}}{72}, \quad \int_{x_{2}}^{y_{2}} \varphi_{3}(x) \mathrm{d} x=\frac{h^{8}}{72},∫x1a1f1(x)dx=−h372,∫x1x1f3(x)dx=3h372,∫x2and2f3(x)dx=h872,
∫ k − 1 x k φ k ( x ) d x = 0 ( k = 4 , 5 , … ) . ∫ k − 1 x k   φ k ( x ) d x = 0 ( k = 4 , 5 , … ) . int_(k-1)^(x_(k))varphi_(k)(x)dx=0quad(k=4,5,dots).\int_{k-1}^{x_{k}} \varphi_{k}(x) \mathrm{d} x=0 \quad(k=4,5, \ldots) .∫k−1xkfk(x)dx=0(k=4,5,…).
The degree of accuracy of formula (C) is 1, because in formula (1) we have
∫ − 4 x n φ ( x ) d x = h 3 12 . ∫ − 4 x n   φ ( x ) d x = h 3 12 . int_(-4)^(x_(n))varphi(x)dx=(h^(3))/(12).\int_{-4}^{x_{n}} \varphi(x) \mathrm{d} x=\frac{h^{3}}{12} .∫−4xnf(x)dx=h312.
From formulas (1) an evaluation of s is deduced R c R c R_(c)R_{c}Rcand name
(2) | R c | ⩽ K c M 2 h 3 , (2) R c ⩽ K c M 2 h 3 , {:(2)|R_(c)| <= K_(c)M_(2)h^(3)",":}\begin{equation*} \left|R_{c}\right| \leqslant K_{c} M_{2} h^{3}, \tag{2} \end{equation*}(2)|Rc|⩽KcM2h3,
where 8-noted
(3) M 2 = sup [ x n = x n ] | f ′ ′ ( x ) | (3) M 2 = sup x n = x n   f ′ ′ ( x ) {:(3)M_(2)=s u p_([x_(n)=x_(n)])|f^('')(x)|:}\begin{equation*} M_{2}=\sup _{\left[x_{n}=x_{n}\right]}\left|f^{\prime \prime}(x)\right| \tag{3} \end{equation*}(3)M2=sup[xn=xn]|f′′(x)|
and
(4) R c h s = ∫ α 0 x n | φ ( x ) | d ω (4) R c h s = ∫ α 0 x n   | φ ( x ) | d ω {:(4)R_(c)h^(s)=int_(alpha_(0))^(x_(n))|varphi(x)|d omega:}\begin{equation*} R_{c} h^{s}=\int_{\alpha_{0}}^{x_{n}}|\varphi(x)| d \omega \tag{4} \end{equation*}(4)Rchs=∫a0xn|f(x)|doh
some practical formulas for quadrature
Wind calculations can be found
∫ A 0 x 1 | φ 1 ( x ) | d x = 1 34 711 992 h 8 , ∫ x 4 x 2 | φ e ( x ) | d x = 1459 34992 h 2 ∫ x 4 x 2 | φ 1 ( x ) | d x = 486 + 73 73 34992 h 3 , ∫ x 3 z 4 | φ 4 ( x ) | d x = 3 54 h 2 . ∫ A 0 x 1   φ 1 ( x ) d x = 1 34 711 992 h 8 , ∫ x 4 x 2   φ e ( x ) d x = 1459 34992 h 2 ∫ x 4 x 2   φ 1 ( x ) d x = 486 + 73 73 34992 h 3 , ∫ x 3 z 4   φ 4 ( x ) d x = 3 54 h 2 . {:[int_(A_(0))^(x_(1))|varphi_(1)(x)|dx=(1)/(34)(711)/(992)h^(8)","quadint_(x_(4))^(x_(2))|varphi_(e)(x)|dx=(1459)/(34992)h^(2)],[int_(x_(4))^(x_(2))|varphi_(1)(x)|dx=(486+73sqrt73)/(34992)h^(3)","quadint_(x_(3))^(z_(4))|varphi_(4)(x)|dx=(sqrt3)/(54)h^(2).]:}\begin{aligned} & \int_{A_{0}}^{x_{1}}\left|\varphi_{1}(x)\right| \mathrm{d} x=\frac{1}{34} \frac{711}{992} h^{8}, \quad \int_{x_{4}}^{x_{2}}\left|\varphi_{\mathrm{e}}(x)\right| \mathrm{d} x=\frac{1459}{34992} h^{2} \\ & \int_{x_{4}}^{x_{2}}\left|\varphi_{1}(x)\right| \mathrm{d} x=\frac{486+73 \sqrt{73}}{34992} h^{3}, \quad \int_{x_{3}}^{z_{4}}\left|\varphi_{4}(x)\right| \mathrm{d} x=\frac{\sqrt{3}}{54} h^{2} . \end{aligned}∫A0x1|f1(x)|dx=134711992h8,∫x4x2|fand(x)|dx=145934992h2∫x4x2|f1(x)|dx=486+737334992h3,∫x3With4|f4(x)|dx=354h2.
Taking into account the symmetry of the curve y = φ ( x ) y = φ ( x ) y=varphi(x)y=\varphi(x)and=f(x)right side x = x = x=x=x= got x 2 x 2 x_(2)x_{2}x2, we have
∫ π ∗ π ∗ | φ ( x ) | d x = [ 2 1711 34992 + 2 ′ 1459 24 θ 92 + 2 − 486 + 73 / 7 3 ― 34982 + ( n − 6 ) 3 54 | A ∫ π ∗ π ∗   | φ ( x ) | d x = 2 1711 34992 + 2 ′ 1459 24 θ 92 + 2 − 486 + 73 / 7 3 ¯ 34982 + ( n − 6 ) 3 54 A int_(pi_(**))^(pi_(**))|varphi(x)|dx=[2(1711)/(34992)+2^(')(1459)/(24 theta92)+2-(486+73//7 bar(3))/(34982)+(n-6)(sqrt3)/(54)|_(A):}\int_{\pi_{*}}^{\pi_{*}}|\varphi(x)| \mathrm{d} x=\left[2 \frac{1711}{34992}+2^{\prime} \frac{1459}{24 \theta 92}+2-\frac{486+73 / 7 \overline{3}}{34982}+\left.(n-6) \frac{\sqrt{3}}{54}\right|_{\mathrm{A}}\right.∫p∗p∗|f(x)|dx=[2171134992+2′145924i92+2−486+73/73―34982+(n−6)354|Aand therefore
(b) K c = 3650 + 73 73 − 1944 3 17496 + n ψ 3 ― 54 (b) K c = 3650 + 73 73 − 1944 3 17496 + n ψ 3 ¯ 54 {:(b)K_(c)=(3650+73sqrt73-1944sqrt3)/(17496)+(n psi bar(3))/(54):}\begin{equation*} K_{c}=\frac{3650+73 \sqrt{73}-1944 \sqrt{3}}{17496}+\frac{n \psi \overline{3}}{54} \tag{b} \end{equation*}(b)Kc=3650+7373−1944317496+nψ3―54
  1. HESTUL IN FORMULA (D)
Assuming that the function f ( x ) f ( x ) f(x)f(x)f(x)it is annoying C 2 [ x 01 x n ] C 2 x 01 x n C^(2)[x_(01)x_(n)]C^{2}\left[x_{01} x_{n}\right]C2[x01xn]the rest is deducted R D R D R_(D)R_{D}RDhas the form
( θ ) R D = ∫ n 0 x m ψ ( x ) f ′ ′ ( x ) d x , ( θ ) R D = ∫ n 0 x m   ψ ( x ) f ′ ′ ( x ) d x , {:(theta")"R_(D)=int_(n_(0))^(x_(m))psi(x)f^('')(x)dx",":}\begin{equation*} R_{D}=\int_{n_{0}}^{x_{m}} \psi(x) f^{\prime \prime}(x) \mathrm{d} x, \tag{$\theta$} \end{equation*}(i)RD=∫n0xmψ(x)f′′(x)dx,
where the function ψ ( x ) ψ ( x ) psi(x)\psi(x)ψ(x)coincides on the intervals [ x 0 , x 1 ] , [ x , x 1 ] … x 0 , x 1 , x , x 1 … [x_(0),x_(1)],[x,x_(1)]dots\left[x_{0}, x_{1}\right],\left[x, x_{1}\right] \ldots[x0,x1],[x,x1]…I polynomials ψ 1 ( x ) , ψ 2 ( x ) , … ψ 1 ( x ) , ψ 2 ( x ) , … psi_(1)(x),psi_(2)(x),dots\psi_{1}(x), \psi_{2}(x), \ldotsψ1(x),ψ2(x),…given by the equations
ψ 1 ( x ) = ( x − x 0 ) 1 2 − 5 h 12 ( x − x 0 ) , ψ k ( x ) = ( x − x k − 1 ) 1 2 − h 2 ( x − x 1 − 1 ) + h 1 12 ψ 1 ( x ) = x − x 0 1 2 − 5 h 12 x − x 0 , ψ k ( x ) = x − x k − 1 1 2 − h 2 x − x 1 − 1 + h 1 12 psi_(1)(x)=((x-x_(0))^(1))/(2)-(5h)/(12)(x-x_(0)),quadpsi_(k)(x)=((x-x_(k-1))^(1))/(2)-(h)/(2)(x-x_(1-1))+(h^(1))/(12)\psi_{1}(x)=\frac{\left(x-x_{0}\right)^{1}}{2}-\frac{5 h}{12}\left(x-x_{0}\right), \quad \psi_{k}(x)=\frac{\left(x-x_{k-1}\right)^{1}}{2}-\frac{h}{2}\left(x-x_{1-1}\right)+\frac{h^{1}}{12}ψ1(x)=(x−x0)12−5h12(x−x0),ψk(x)=(x−xk−1)12−h2(x−x1−1)+h112
( k = 2 , 3 , ⋯ ) . ( k = 2 , 3 , ⋯ ) . (k=2,3,cdots).(k=2,3, \cdots) .(k=2,3,⋯).
It is found that
∫ x 1 x 1 ψ 1 ( x ) d x = − h 3 24 , ∫ k + 1 x k ψ 1 ( x ) d x = 0 ( k = 2 , 3 , … ) ∫ x 1 x 1   ψ 1 ( x ) d x = − h 3 24 , ∫ k + 1 x k   ψ 1 ( x ) d x = 0 ( k = 2 , 3 , … ) int_(x_(1))^(x_(1))psi_(1)(x)dx=-(h^(3))/(24),quadint_(k+1)^(x_(k))psi_(1)(x)dx=0quad(k=2,3,dots)\int_{x_{1}}^{x_{1}} \psi_{1}(x) \mathrm{d} x=-\frac{h^{3}}{24}, \quad \int_{k+1}^{x_{k}} \psi_{1}(x) \mathrm{d} x=0 \quad(k=2,3, \ldots)∫x1x1ψ1(x)dx=−h324,∫k+1xkψ1(x)dx=0(k=2,3,…)
and therefore the result is
∫ x 0 x n ψ ( x ) d x = − h 3 12 , ∫ x 0 x n   ψ ( x ) d x = − h 3 12 , int_(x_(0))^(x_(n))psi(x)dx=-(h^(3))/(12),\int_{x_{0}}^{x_{n}} \psi(x) \mathrm{d} x=-\frac{h^{3}}{12},∫x0xnψ(x)dx=−h312,
which means that the degree of accuracy of formula (D) is 1.
From formula (0) it is deduced that
(T) | R D | ⩽ K D M B h 3 (T) R D ⩽ K D M B h 3 {:(T)|R_(D)| <= K_(D)M_(B)h^(3):}\begin{equation*} \left|R_{D}\right| \leqslant K_{D} M_{B} h^{3} \tag{T} \end{equation*}(T)|RD|⩽KDMBh3
ande K b A a = ∫ M 0 x 0 | Ψ ( x ) | d x .  ande  K b A a = ∫ M 0 x 0   | Ψ ( x ) | d x . " ande "K_(b)A^(a)=int_(M_(0))^(x_(0))|Psi(x)|dx.\text { ande } K_{b} A^{\mathrm{a}}=\int_{M_{0}}^{x_{0}}|\Psi(x)| d x . Spirit KbAa=∫M0x0|P(x)|dx.
Doing the calculations, we find
∫ x 0 π 1 | ψ 1 ( x ) | d x = 71 k 8 1296 , ∫ α 1 π 2 | ψ 8 ( x ) | d x = V 3 ― 54 h 0 ∫ x 0 π 1   ψ 1 ( x ) d x = 71 k 8 1296 , ∫ α 1 π 2   ψ 8 ( x ) d x = V 3 ¯ 54 h 0 int_(x_(0))^(pi_(1))|psi_(1)(x)|dx=(71k^(8))/(1296),quadint_(alpha_(1))^(pi_(2))|psi_(8)(x)|dx=(V bar(3))/(54)h^(0)\int_{x_{0}}^{\pi_{1}}\left|\psi_{1}(x)\right| d x=\frac{71 k^{8}}{1296}, \quad \int_{\alpha_{1}}^{\pi_{2}}\left|\psi_{8}(x)\right| d x=\frac{V \overline{3}}{54} h^{0}∫x0p1|ψ1(x)|dx=71k81296,∫a1p2|ψ8(x)|dx=In3―54h0
and taking into account the symmetry of the curve y = ψ ( x ) y = ψ ( x ) y=psi(x)y=\psi(x)and=ψ(x)right side x = x 0 + x 0 2 x = x 0 + x 0 2 x=(x_(0)+x_(0))/(2)x=\frac{x_{0}+x_{0}}{2}x=x0+x02HAVE
∫ x 0 x 0 [ ψ ( x ) | d x = [ 2 71 1296 + ( π − 2 ) V 3 ― 54 ] A 3 = ( 71 − 12 V 3 648 + n V 3 ― 54 ) h 3 ∫ x 0 x 0   ψ ( x ) d x = 2 71 1296 + ( π − 2 ) V 3 ¯ 54 A 3 = 71 − 12 V 3 648 + n V 3 ¯ 54 h 3 int_(x_(0))^(x_(0))[psi(x)|dx=[2(71)/(1296)+(pi-2)(V bar(3))/(54)]A^(3)=((71-12 V3)/(648)+(nV bar(3))/(54))h^(3):}\int_{x_{0}}^{x_{0}}\left[\psi(x) \left\lvert\, \mathrm{d} x=\left[2 \frac{71}{1296}+(\pi-2) \frac{V \overline{3}}{54}\right] A^{3}=\left(\frac{71-12 V 3}{648}+\frac{n V \overline{3}}{54}\right) h^{3}\right.\right.∫x0x0[ψ(x)|dx=[2711296+(p−2)In3―54]A3=(71−12In3648+nIn3―54)h3
(8) K D = 71 − 12 3 648 + n 3 54 . (8) K D = 71 − 12 3 648 + n 3 54 . {:(8)K_(D)=(71-12sqrt3)/(648)+(nsqrt3)/(54).:}\begin{equation*} K_{D}=\frac{71-12 \sqrt{3}}{648}+\frac{n \sqrt{3}}{54} . \tag{8} \end{equation*}(8)KD=71−123648+n354.

1. THE REMAINDER IN THE FORMULA (L)

Presuming that it works f ( x ) f ( x ) f(x)f(x)f(x)it is classy C 2 [ x 01 x p ] C 2 x 01 x p C^(2)[x_(01)x_(p)]C^{2}\left[x_{01} x_{p}\right]C2[x01xp]we can serialize her crime in the form
(9) R L = ∫ x v x n θ ( x ) ∫ ′ ′ ( x ) d x (9) R L = ∫ x v x n   θ ( x ) ∫ ′ ′ ( x ) d x {:(9)R_(L)=int_(x_(v))^(x_(n))theta(x)int^('')(x)dx:}\begin{equation*} R_{L}=\int_{x_{v}}^{x_{n}} \theta(x) \int^{\prime \prime}(x) \mathrm{d} x \tag{9} \end{equation*}(9)RL=∫xinxni(x)∫′′(x)dx
where the function θ ( x ) θ ( x ) theta(x)\theta(x)i(x)coincides on the intervals [ x 0 , x 1 ] [ x 1 , x 1 ] , … x 0 , x 1 x 1 , x 1 , … [x_(0),x_(1)][x_(1),x_(1)],dots\left[x_{0}, x_{1}\right]\left[x_{1}, x_{1}\right], \ldots[x0,x1][x1,x1],…I polynomials
θ 1 ( x ) = ( x − x 0 ) 2 2 − 3 h 8 ( x − x 0 ) , θ 3 ( x ) = ( x − x 1 ) 2 2 − 13 h 24 ( x − x 1 ) + h 2 8 θ k ( x ) = ( x − x k − 1 ) 2 2 − h 2 ( x − x k − 1 ) + h 2 12 ( k = 3 , 4 , … ) θ 1 ( x ) = x − x 0 2 2 − 3 h 8 x − x 0 , θ 3 ( x ) = x − x 1 2 2 − 13 h 24 x − x 1 + h 2 8 θ k ( x ) = x − x k − 1 2 2 − h 2 x − x k − 1 + h 2 12 ( k = 3 , 4 , … ) {:[theta_(1)(x)=((x-x_(0))^(2))/(2)-(3h)/(8)(x-x_(0))","quadtheta_(3)(x)=((x-x_(1))^(2))/(2)-(13 h)/(24)(x-x_(1))+(h^(2))/(8)],[theta_(k)(x)=((x-x_(k-1))^(2))/(2)-(h)/(2)(x-x_(k-1))+(h^(2))/(12)quad(k=3","4","dots)]:}\begin{gathered} \theta_{1}(x)=\frac{\left(x-x_{0}\right)^{2}}{2}-\frac{3 h}{8}\left(x-x_{0}\right), \quad \theta_{3}(x)=\frac{\left(x-x_{1}\right)^{2}}{2}-\frac{13 h}{24}\left(x-x_{1}\right)+\frac{h^{2}}{8} \\ \theta_{k}(x)=\frac{\left(x-x_{k-1}\right)^{2}}{2}-\frac{h}{2}\left(x-x_{k-1}\right)+\frac{h^{2}}{12} \quad(k=3,4, \ldots) \end{gathered}i1(x)=(x−x0)22−3h8(x−x0),i3(x)=(x−x1)22−13h24(x−x1)+h28ik(x)=(x−xk−1)22−h2(x−xk−1)+h212(k=3,4,…)
we
∫ x 0 x 1 θ 1 ( x ) d x = − h 8 48 , ∫ x 1 x 1 θ 2 ( x ) d x = h 8 48 ⋅ ∫ x 1 x 1 θ 3 ( x ) d x = 0 . ∫ x 0 x 1   θ 1 ( x ) d x = − h 8 48 , ∫ x 1 x 1   θ 2 ( x ) d x = h 8 48 ⋅ ∫ x 1 x 1   θ 3 ( x ) d x = 0 . int_(x_(0))^(x_(1))theta_(1)(x)dx=-(h^(8))/(48),quadint_(x_(1))^(x_(1))theta_(2)(x)dx=(h^(8))/(48)*int_(x_(1))^(x_(1))theta_(3)(x)dx=0.\int_{x_{0}}^{x_{1}} \theta_{1}(x) \mathrm{d} x=-\frac{h^{8}}{48}, \quad \int_{x_{1}}^{x_{1}} \theta_{2}(x) \mathrm{d} x=\frac{h^{8}}{48} \cdot \int_{x_{1}}^{x_{1}} \theta_{3}(x) \mathrm{d} x=0 .∫x0x1i1(x)dx=−h848,∫x1x1i2(x)dx=h848⋅∫x1x1i3(x)dx=0.
Result horse
∫ x 0 z u θ ( x ) d x = 0 ∫ x 0 z u   θ ( x ) d x = 0 int_(x_(0))^(z_(u))theta(x)dx=0\int_{x_{0}}^{z_{u}} \theta(x) \mathrm{d} x=0∫x0Withini(x)dx=0
which means that the degree of enactity of the formula (L) is much greater than 1. It is easily seen that the degree of enactity of the formula (L) is ib.
SOME PRACTICAL FORMULAS OF SQUAREDNESS
603
From formula (9) it follows that
| R t | ⩽ K L M 2 A ∗ K t h 3 = ∫ x 0 x ∗ | θ ( x ) | d x R t ⩽ K L M 2 A ∗ K t h 3 = ∫ x 0 x ∗   | θ ( x ) | d x {:[|R_(t)| <= K_(L)M_(2)A^(**)],[K_(t)h^(3)=int_(x_(0))^(x_(**))|theta(x)|dx]:}\begin{gathered} \left|R_{t}\right| \leqslant K_{L} M_{2} A^{*} \\ K_{t} h^{3}=\int_{x_{0}}^{x_{*}}|\theta(x)| d x \end{gathered}|Rt|⩽KLM2A∗Kth3=∫x0x∗|i(x)|dx
Doing the calculations we find
∫ x 1 x 1 | 0 1 ( x ) | d x = 19 384 h 8 , ∫ α 1 x 1 | θ 2 ( x ) | d x = 341 10368 h 2 , ∫ α 0 α 1 | θ 1 ( x ) | d x = 3 54 h 1 ∫ x 1 x 1   0 1 ( x ) d x = 19 384 h 8 , ∫ α 1 x 1   θ 2 ( x ) d x = 341 10368 h 2 , ∫ α 0 α 1   θ 1 ( x ) d x = 3 54 h 1 int_(x_(1))^(x_(1))|0_(1)(x)|dx=(19)/(384)h^(8),quadint_(alpha_(1))^(x_(1))|theta_(2)(x)|dx=(341)/(10368)h^(2),int_(alpha_(0))^(alpha_(1))|theta_(1)(x)|dx=(sqrt3)/(54)h^(1)\int_{x_{1}}^{x_{1}}\left|0_{1}(x)\right| \mathrm{d} x=\frac{19}{384} h^{8}, \quad \int_{\alpha_{1}}^{x_{1}}\left|\theta_{2}(x)\right| \mathrm{d} x=\frac{341}{10368} h^{2}, \int_{\alpha_{0}}^{\alpha_{1}}\left|\theta_{1}(x)\right| \mathrm{d} x=\frac{\sqrt{3}}{54} h^{1}∫x1x1|01(x)|dx=19384h8,∫a1x1|i2(x)|dx=34110368h2,∫a0a1|i1(x)|dx=354h1.
It is a rule that
∫ x n x n | θ ( x ) | d x = 2 ( 19 384 + 341 10368 ) h 3 + ( n − 4 ) 3 54 h 2 = 427 2592 h 2 + + ( n − 4 ) 3 54 h 2 ∫ x n x n   | θ ( x ) | d x = 2 19 384 + 341 10368 h 3 + ( n − 4 ) 3 54 h 2 = 427 2592 h 2 + + ( n − 4 ) 3 54 h 2 {:[int_(x_(n))^(x_(n))|theta(x)|dx=2((19)/(384)+(341)/(10368))h^(3)+(n-4)(sqrt3)/(54)h^(2)=(427)/(2592)h^(2)+],[+(n-4)(sqrt3)/(54)h^(2)]:}\begin{gathered} \int_{x_{n}}^{x_{n}}|\theta(x)| d x=2\left(\frac{19}{384}+\frac{341}{10368}\right) h^{3}+(n-4) \frac{\sqrt{3}}{54} h^{2}=\frac{427}{2592} h^{2}+ \\ +(n-4) \frac{\sqrt{3}}{54} h^{2} \end{gathered}∫xnxn|i(x)|dx=2(19384+34110368)h3+(n−4)354h2=4272592h2++(n−4)354h2
I declare
(11) K L = 427 − 192 3 2 ― 592 + n 3 54 . (11) K L = 427 − 192 3 2 ¯ 592 + n 3 54 . {:(11)K_(L)=(427-192sqrt3)/( bar(2)592)+n(sqrt3)/(54).:}\begin{equation*} K_{L}=\frac{427-192 \sqrt{3}}{\overline{2} 592}+n \frac{\sqrt{3}}{54} . \tag{11} \end{equation*}(11)KL=427−19232―592+n354.
  1. COMPARISON OF NUMBERS K C , K D , K L K C , K D , K L K_(C),K_(D),K_(L)K_{C}, K_{D}, K_{L}KC,KD,KL, WITH THE NUMBER K , K ,  K_(", ")K_{\text {, }}K, CORRESPONDING TO THE FORMULA (T)
As is known in formula (T), assuming that the function aff(x) is of class C 1 [ x 0 , x n ] C 1 x 0 , x n C^(1)[x_(0),x_(n)]C^{1}\left[x_{0}, x_{n}\right]C1[x0,xn], evem
(12) R T = ∫ x 0 a x x ( x ) f ′ ′ ( x ) d x (12) R T = ∫ x 0 a x   x ( x ) f ′ ′ ( x ) d x {:(12)R_(T)=int_(x_(0))^(a_(x))x(x)f^('')(x)dx:}\begin{equation*} R_{T}=\int_{x_{0}}^{a_{x}} x(x) f^{\prime \prime}(x) d x \tag{12} \end{equation*}(12)RT=∫x0axx(x)f′′(x)dx
where the function χ ( x ) χ ( x ) chi(x)\chi(x)x(x)coincides on the intervals [ x 0 , x 1 ] , [ x 1 , x 0 ] , … x 0 , x 1 , x 1 , x 0 , … [x_(0),x_(1)],[x_(1),x_(0)],dots\left[x_{0}, x_{1}\right],\left[x_{1}, x_{0}\right], \ldots[x0,x1],[x1,x0],…with polynomials
χ 1 ( x ) = ( x − x 0 ) ( x − x 1 ) 2 , χ 1 ( x ) = ( x − x 1 ) ( x − x 1 ) 2 , … χ 1 ( x ) = x − x 0 x − x 1 2 , χ 1 ( x ) = x − x 1 x − x 1 2 , … chi_(1)(x)=((x-x_(0))(x-x_(1)))/(2),quadchi_(1)(x)=((x-x_(1))(x-x_(1)))/(2),dots\chi_{1}(x)=\frac{\left(x-x_{0}\right)\left(x-x_{1}\right)}{2}, \quad \chi_{1}(x)=\frac{\left(x-x_{1}\right)\left(x-x_{1}\right)}{2}, \ldotsx1(x)=(x−x0)(x−x1)2,x1(x)=(x−x1)(x−x1)2,…
It is noted that
(13) ( R T ) ⩽ K T M 2 h a (13) R T ⩽ K T M 2 h a {:(13)(R_(T)) <= K_(T)M_(2)h^(a):}\begin{equation*} \left(R_{T}\right) \leqslant K_{T} M_{2} h^{\mathrm{a}} \tag{13} \end{equation*}(13)(RT)⩽KTM2ha
where
((14) K r = n 12 . ((14) K r = n 12 . {:((14)K_(r)=(n)/( 12).:}\begin{equation*} K_{r}=\frac{n}{12} . \tag{(14} \end{equation*}((14)Kr=n12.
Let's see what we have.
K C − K L = 1095 + 292 | 73 − 2592 | 3 69984 = 0 , 0157 … , K p − K c = 1620 + 3 − 73 / 73 − 1739 = 0 , 0253 . K C − K L = 1095 + 292 | 73 − 2592 | 3 69984 = 0 , 0157 … , K p − K c = 1620 + 3 − 73 / 73 − 1739 = 0 , 0253 . {:[K_(C)-K_(L)=(1095+292|73-2592|3)/(69984)=0","0157 dots","],[K_(p)-K_(c)=1620+3-73//73-1739=0","0253.]:}\begin{gathered} K_{C}-K_{L}=\frac{1095+292|73-2592| 3}{69984}=0,0157 \ldots, \\ K_{p}-K_{c}=1620+3-73 / 73-1739=0,0253 . \end{gathered}KC−KL=1095+292|73−2592|369984=0,0157…,Kp−Kc=1620+3−73/73−1739=0,0253.
It is also found that we have
K r − K D = ( 54 − 12 ∣ 3 ) ( n − 2 ) + ( 37 − 12 ∣ 3 ― ) 648 > 0 ( n ⩾ 2 ) K r − K D = ( 54 − 12 ∣ 3 ) ( n − 2 ) + ( 37 − 12 ∣ 3 ¯ ) 648 > 0 ( n ⩾ 2 ) K_(r)-K_(D)=((54-12∣3)(n-2)+(37-12∣ bar(3)))/(648) > 0quad(n >= 2)K_{r}-K_{D}=\frac{(54-12 \mid 3)(n-2)+(37-12 \mid \overline{3})}{648}>0 \quad(n \geqslant 2)Kr−KD=(54−12∣3)(n−2)+(37−12∣3―)648>0(n⩾2)
from which it follows that we have
(15) K t < K r < K D < K T . (15) K t < K r < K D < K T . {:(15)K_(t) < K_(r) < K_(D) < K_(T).:}\begin{equation*} \boldsymbol{K}_{t}<\boldsymbol{K}_{r}<\boldsymbol{K}_{\mathrm{D}}<\boldsymbol{K}_{T} . \tag{15} \end{equation*}(15)Kt<Kr<KD<KT.
In conclusion, we have determined and evaluated the residues in formulas (C), (D), (L). From inequalities (15) it follows that the more advantageous is formula (L), which immediately follows, with a very small difference for the coefficient K K KKK, formula (C). The most disadvantageous of all the formulas considered in this note is formula (T).
We will now introduce a method for constructing formulas (C), (D), (L) and other analogous ones.
RPR Academy Cluj Branch,
Computing Institute
SUMMARY
The remainders of the quadrature formulas (C), (D), (L), (T) are determined by the formulas (1), (6), (9), (12), and the evaluations (2), (7), (10), (13), characterized by the numbers K G , K D , K L , K T K G , K D , K L , K T K_(G),K_(D),K_(L),K_(T)K_{G}, K_{D}, K_{L}, K_{T}KG,KD,KL,KTgiven by the formulas
Cateva formula practice of cuadratura
895
(1), (14). We demonstrate the inequalities, (15), from which it follows that (1). Ja (1) is the most advantageous, followed by formula (C), and that formula (T) is the most disadvantageous.
BIBLIOGRAPHY
G Goblarl Operalions sur les commes experimeniales, CF Aced. Sel. Fatii, 1058, e46, 1799-1800.
Minkull, Technique de Caleut namérique, Libralite Polylechnique Ch. Béranger, Paris, 1052, 244.
3. DV IONFBCT, CDOdraturi numerice, Edit. Lehnica, Bucurest, 1957.
(() => { window.MathpixRenderConfig = {"width":1200,"outMath":{"output_format":"svg","include_svg":true,"include_asciimath":true,"include_latex":true,"include_mathml":true,"include_mathml_word":true,"include_tsv":true,"include_table_markdown":true,"include_smiles":true},"accessibility":{"assistive_mml":true,"include_speech":false}}; const autoScript = document.createElement('script'); autoScript.src = "https://cdn.jsdelivr.net/npm/mathpix-markdown-it@2.0.37/es5/browser/auto-render.js"; autoScript.onload = () => { console.log("[Mathpix] auto-render loaded"); }; autoScript.onerror = () => { console.error("[Mathpix] Failed to load auto-render script:", autoScript.src); }; document.head.append(autoScript); })(); (() => { const ctxScript = document.createElement('script'); ctxScript.src = "https://cdn.jsdelivr.net/npm/mathpix-markdown-it@2.0.37/es5/context-menu.js"; ctxScript.onload = () => { const init = () => { window.addListenerContextMenuEvents?.(); }; if (document.readyState === "loading") { document.addEventListener("DOMContentLoaded", init, { once: true }); } else { init(); } if (!window.__mathpixContextMenuKeyHandlerInstalled) { window.__mathpixContextMenuKeyHandlerInstalled = true; document.addEventListener("keydown", (e) => { if (e.key !== "ArrowUp" && e.key !== "ArrowDown") return; const el = document.querySelector(".mmd-menu"); if (el) e.preventDefault(); }); } }; ctxScript.onerror = () => { console.error("[Mathpix] Failed to load context-menu script:", ctxScript.src); }; document.head.appendChild(ctxScript); })();
1963

Related Posts