Let \(X,Y\) be normed spaces, \(G:X\rightarrow Y\) a nonlinear operator, and the nonlinar equation \(G\left( x\right) =0\). We define the divided differences of \(G\) and we give some examples of nonlinear operators on Banach spaces for which we construct the divided differences of different orders. We construct the Lagrange interpolation polynomial in the Newton form for \(G\). The solution of equation is a approximated by using the inverse interpolation Lagrange polynomial.
Ion Păvăloiu
(Tiberiu Popoviciu Institute of Numerical Analysis)
Original title (in French)
La résolution des equations par intérpolation
English translation of the title
Solving equations by interpolation
divided difference in normed spaces, Lagrange inverse interpolation; multistep iterative method of Newton type
I. Păvăloiu, La résolution des equations par intérpolation, Mathématica, 23(46) (1981) no. 1, pp. 61-72 (in French).
[1] A.M. Ostrowski, Resenie uravnenii i sistem uravnenii, Izd. inostr. lit. Moskva (1963).
[2] I. Pavaloiu, Interpolation dans les espaces lineaires normes et applications, Mathematica, Cluj, nr.12 (35), 2, (1970), pp. 309–324.
[3] I. Pavaloiu, Introducere in teoria aproximarii solutiilor ecuatiilor, Editura Dacia, 1976.
[4] T. Popoviciu, Introduction a la theorie des differences divisees, Bulletin mathematique de la societe Roumaine de Sciences, 42, 1, (1940) pp. 65–78.
[5] A.S. Sergeev, O metode hord, Sibirski mat. Jurnal, XI, (2), (1961), pp. 282–289.
[6] J.F. Traub, Iterative Methods for the Solution of Equations, Prentice Hall, Series in Automatic Computation, 1964.
[7] S. Ul’m, Ob obobscennyh razdelennîh raznostiah, II, Izv. Nauk Estonskoi SSR, 16, 2, (1967), 146–155.
/*======================================================================
Core CSS for LaTeXML documents converted to (X)HTML */
/* Generic Page layout */
.ltx_page_header,
.ltx_page_footer { font-size:0.8em; }
.ltx_page_header *[rel~="prev"],
.ltx_page_footer *[rel~="prev"] { float:left; }
.ltx_page_header *[rel~="up"],
.ltx_page_footer *[rel~="up"] { display:block; text-align:center; }
.ltx_page_header *[rel~="next"],
.ltx_page_footer *[rel~="next"] { float:right; }
/* What was I trying for here; need more selective rule!
.ltx_page_header .ltx_ref,
.ltx_page_footer .ltx_ref {
margin:0 1em; }
*/
.ltx_page_header li {
padding:0.1em 0.2em 0.1em 1em;}
/* Main content */
.ltx_page_content { clear:both; }
.ltx_page_header { border-bottom:1px solid; margin-bottom:5px; }
.ltx_page_footer { clear:both; border-top:1px solid; margin-top:5px; }
.ltx_page_header:after,
.ltx_page_footer:after,
.ltx_page_content:after {
content:"."; display:block; height:0; clear:both; visibility:hidden; }
.ltx_page_footer:before {
content:"."; display:block; height:0; clear:both; visibility:hidden; }
.ltx_page_logo { font-size:80%; margin-top: 5px; clear:both; float:right; }
.ltx_page_logo a { font-variant: small-caps; }
.ltx_page_logo img { vertical-align:-3px; }
/* if shown */
.ltx_page_navbar li { white-space:nowrap; display:block; overflow:hidden; }
/* If ref got turned into span, it's "this section"*/
.ltx_page_navbar li span.ltx_ref { white-space:normal; overflow:visible; }
/* Ought to be easily removable/overridable? */
.ltx_pagination.ltx_role_newpage { height:2em; }
/*======================================================================
Document Structure; Titles & Frontmatter */
/* undo bold here to remove the browser's native h# styling,
at let all other styles override it (with more specific rules)*/
.ltx_title { font-size:100%; font-weight:normal; }
/* Hack to simulate run-in! put class="ltx_runin" on a title or tag
for it to run-into the following text. */
.ltx_runin { display:inline; }
.ltx_runin:after { content:" "; }
.ltx_runin + .ltx_para,
.ltx_runin + .ltx_para p,
.ltx_runin + p { display:inline; }
.ltx_outdent { margin-left: -2em; }
/* .ltx_chapter_title, etc should be in ltx-article.css etc.
*/
.ltx_page_main { margin:0px; padding:1em 3em 1em 2em; }
.ltx_tocentry { list-style-type:none; }
/* support for common author block layouts.*/
/* add class ltx_authors_1line to get authors in single line
with pop-up affiliation, etc. */
.ltx_authors_1line .ltx_creator,
.ltx_authors_1line .ltx_author_before,
.ltx_authors_1line .ltx_author_after { display:inline;}
.ltx_authors_1line .ltx_author_notes { display:inline-block; }
.ltx_authors_1line .ltx_author_notes:before { content:"*"; color:blue;}
.ltx_authors_1line .ltx_author_notes span { display:none; }
.ltx_authors_1line .ltx_author_notes:hover span {
display:block; position:absolute; z-index:10;
background:white; text-align:left;
border: 1px solid black; border-radius: 0 5px 5px 5px; box-shadow: 5px 5px 10px gray; }
/* add class=ltx_authors_multiline to get authors & affliations on separate lines*/
.ltx_authors_multiline .ltx_creator,
.ltx_authors_multiline .ltx_author_before,
.ltx_authors_multiline .ltx_author_after,
.ltx_authors_multiline .ltx_author_notes,
.ltx_authors_multiline .ltx_author_notes .ltx_contact {
display:block; }
/*======================================================================
Para level */
.ltx_float {
margin: 1ex 3em 1ex 3em; }
td.ltx_subfigure,
td.ltx_subtable,
td.ltx_subfloat { width:50%; }
/* theorems, figure, tables, floats captions.. */
/*======================================================================
Blocks, Lists, Floats */
.ltx_p,
.ltx_quote,
.ltx_block,
.ltx_para {
display: block; }
/* alignment within blocks */
.ltx_align_left { text-align:left; }
.ltx_align_right { text-align:right; }
.ltx_align_center { text-align:center; }
.ltx_align_justify { text-align:justify; }
.ltx_align_top { vertical-align:top; }
.ltx_align_bottom { vertical-align:bottom; }
.ltx_align_middle { vertical-align:middle; }
.ltx_align_baseline { vertical-align:baseline; }
.ltx_align_floatleft { float:left; }
.ltx_align_floatright { float:right; }
.ltx_td.ltx_align_left, .ltx_th.ltx_align_left,
.ltx_td.ltx_align_right, .ltx_th.ltx_align_right,
.ltx_td.ltx_align_center, .ltx_th.ltx_align_center { white-space:nowrap; }
.ltx_td.ltx_align_left.ltx_wrap, .ltx_th.ltx_align_left.ltx_wrap,
.ltx_td.ltx_align_right.ltx_wrap, .ltx_th.ltx_align_right.ltx_wrap,
.ltx_td.ltx_align_center.ltx_wrap, .ltx_th.ltx_align_center.ltx_wrap,
.ltx_td.ltx_align_justify, .ltx_th.ltx_align_justify { white-space:normal; }
.ltx_tabular .ltx_tabular { width:100%; }
.ltx_inline-block { display:inline-block; }
/* equations in non-aligned mode (not normally used) */
.ltx_eqn_div { display:block; width:95%; text-align:center; }
/* equations in aligned mode (aligning tags, etc as well as equations) */
.ltx_eqn_table { display:table; width:100%; border-collapse:collapse; }
.ltx_eqn_row { display:table-row; }
.ltx_eqn_cell { display:table-cell; width:auto; }
/* Padding between column pairs in ams align */
table.ltx_eqn_align tr.ltx_equation td.ltx_align_left + td.ltx_align_right,
table.ltx_eqn_align tr.ltx_equation td.ltx_align_left + td.ltx_align_center,
table.ltx_eqn_align tr.ltx_equation td.ltx_align_center + td.ltx_align_right,
table.ltx_eqn_align tr.ltx_equation td.ltx_align_center + td.ltx_align_center { padding-left:3em; }
table.ltx_eqn_eqnarray tr.ltx_eqn_lefteqn + tr td.ltx_align_right { min-width:2em; }
.ltx_eqn_eqno { max-width:0em; overflow:visible; white-space: nowrap; }
.ltx_eqn_eqno.ltx_align_right .ltx_tag { float:right; }
.ltx_eqn_center_padleft,
.ltx_eqn_center_padright { width:50%; min-width:2em;}
.ltx_eqn_left_padleft,
.ltx_eqn_right_padright { min-width:2em; }
.ltx_eqn_left_padright,
.ltx_eqn_right_padleft { width:100%; }
/* Various lists */
.ltx_itemize,
.ltx_enumerate,
.ltx_description {
display:block; }
.ltx_itemize .ltx_item,
.ltx_enumerate .ltx_item {
display: list-item; }
/* Position the tag to look like a normal item bullet. */
li.ltx_item > .ltx_tag {
display:inline-block; margin-left:-1.5em; min-width:1.5em;
text-align:right; }
.ltx_item .ltx_tag + .ltx_para,
.ltx_item .ltx_tag + .ltx_para .ltx_p { display:inline; }
/* NOTE: Need to try harder to get runin appearance? */
dl.ltx_description dt { margin-right:0.5em; float:left;
font-weight:bold; font-size:95%; }
dl.ltx_description dd { margin-left:5em; }
dl.ltx_description dl.ltx_description dd { margin-left:3em; }
/* Theorems */
.ltx_theorem {margin:1em 0em 1em 0em; }
.ltx_title_theorem { font-size:100%; }
/* Bibliographies */
.ltx_bibliography dt { margin-right:0.5em; float:left; }
.ltx_bibliography dd { margin-left:3em; }
/*.ltx_biblist { list-style-type:none; }*/
.ltx_bibitem { list-style-type:none; }
.ltx_bibitem .ltx_tag { font-weight:bold; margin-left:-2em; width:3em; }
/*.bibitem-tag + div { display:inline; }*/
.ltx_bib_title { font-style:italic; }
.ltx_bib_article .bib-title { font-style:normal !important; }
.ltx_bib_journal { font-style:italic; }
.ltx_bib_volume { font-weight:bold; }
/* Indices */
.ltx_indexlist li { list-style-type:none; }
.ltx_indexlist { margin-left:1em; padding-left:1em;}
/* Listings */
.ltx_listing {
display:block;
margin: 1ex 3em 1ex 0em;
overflow-x:auto;
text-align: left; }
.ltx_float .ltx_listing {
margin: 0; }
.ltx_listingline { white-space:nowrap; min-height:1em; }
.ltx_lst_numbers_left .ltx_listingline .ltx_tag {
background-color:transparent;
margin-left:-3em; width:2.5em;
position:absolute;
text-align:right; }
.ltx_lst_numbers_right .ltx_listingline .ltx_tag {
background-color:transparent;
width:2.5em;
position:absolute; right:3em;
text-align:right; }
/*
position:absolute; left:0em;
max-width:0em; text-align:right; }
*/
.ltx_parbox {text-indent:0em; }
/* NOTE that it is CRITICAL to put position:relative outside & absolute inside!!
I wish I understood why!
Outer box establishes resulting size, neutralizes any outer positioning, etc;
inner establishes position of stuff to be rotated */
.ltx_transformed_outer {
position:relative; bottom:0pt;left:0pt;
overflow:visible; }
.ltx_transformed_inner {
display:block;
position:absolute;bottom:0pt;left:0pt; }
.ltx_transformed_inner > .ltx_p {text-indent:0em; margin:0; padding:0; }
/* If simulating a table (html5), try to get rowspan to work...sorta? */
span.ltx_rowspan { position:absolute; top:0; bottom:0; }
/* by default, p doesn't indent */
.ltx_p { text-indent:0em; white-space:normal; }
/* explicit control of indentation (on ltx_para) */
.ltx_indent > .ltx_p:first-child { text-indent:2em!important; }
.ltx_noindent > .ltx_p:first-child { text-indent:0em!important; }
/*======================================================================
Columns */
.ltx_page_column1 {
width:44%; float:left; } /* IE uses % of wrong container*/
.ltx_page_column2 {
width:44%; float:right; }
.ltx_page_columns > .ltx_page_column1 {
width:48%; float:left; }
.ltx_page_columns > .ltx_page_column2 {
width:48%; float:right; }
.ltx_page_columns:after {
content:"."; display:block; height:0; clear:both; visibility:hidden; }
/*======================================================================
Borders and such */
.ltx_tabular { display:inline-table; border-collapse:collapse; }
.ltx_tabular.ltx_centering { display:table; }
.ltx_thead,
.ltx_tfoot,
.ltx_tbody { display:table-row-group; }
.ltx_tr { display:table-row; }
.ltx_td,
.ltx_th { display:table-cell; }
.ltx_framed { border:1px solid black;}
.ltx_tabular .ltx_td,
.ltx_tabular .ltx_th { padding:0.1em 0.5em; }
/* regular lines */
.ltx_border_t { border-top:1px solid black; }
.ltx_border_r { border-right:1px solid black; }
.ltx_border_b { border-bottom:1px solid black; }
.ltx_border_l { border-left:1px solid black; }
/* double lines */
.ltx_border_tt { border-top:3px double black; }
.ltx_border_rr { border-right:3px double black; }
.ltx_border_bb { border-bottom:3px double black; }
.ltx_border_ll { border-left:3px double black; }
/* Light lines */
.ltx_border_T { border-top:1px solid gray; }
.ltx_border_R { border-right:1px solid gray; }
.ltx_border_B { border-bottom:1px solid gray; }
.ltx_border_L { border-left:1px solid gray; }
/* Framing */
.ltx_framed_rectangle { border-style:solid; border-width:1px; }
.ltx_framed_top { border-top-style:solid; border-top-width:1px; }
.ltx_framed_left { border-left-style:solid; border-left-width:1px; }
.ltx_framed_right { border-right-style:solid; border-right-width:1px; }
.ltx_framed_bottom,
.ltx_framed_underline { border-bottom-style:solid; border-bottom-width:1px; }
.ltx_framed_topbottom { border-top-style:solid; border-top-width:1px;
border-bottom-style:solid; border-bottom-width:1px; }
.ltx_framed_leftright { border-left-style:solid; border-left-width:1px;
border-right-style:solid; border-right-width:1px; }
/*======================================================================
Misc */
/* .ltx_verbatim*/
.ltx_verbatim { text-align:left; }
/*======================================================================
Meta stuff, footnotes */
.ltx_note_content { display:none; }
/*right:5%; */
.ltx_note_content {
max-width: 70%; font-size:90%; left:15%;
text-align:left;
background-color: white;
padding: 0.5em 1em 0.5em 1.5em;
border: 1px solid black; border-radius: 0 5px 5px 5px; box-shadow: 5px 5px 10px gray; }
.ltx_note_mark { color:blue; }
.ltx_note_type { font-weight: bold; }
.ltx_note { display:inline-block; text-indent:0; } /* So we establish containing block */
.ltx_note_content .ltx_note_mark { position:absolute; left:0.2em; top:-0.1em; }
.ltx_note:hover .ltx_note_content,
.ltx_note .ltx_note_content:hover {
display:block; position:absolute; z-index:10; }
.ltx_ERROR { color:red; }
.ltx_rdf { display:none; }
.ltx_missing { color:red;}
.ltx_nounicode { color:red; }
/*======================================================================
SVG (pgf/tikz ?) basics */
/* Stuff appearing in svg:foreignObject */
.ltx_svg_fog foreignObject { margin:0; padding:0; overflow:visible; }
.ltx_svg_fog foreignObject > p { margin:0; padding:0; display:block; }
/*.ltx_svg_fog foreignObject > p { margin:0; padding:0; display:block; white-space:nowrap; }*/
/*======================================================================
Low-level Basics */
/* Note that LaTeX(ML)'s font model doesn't map quite exactly to CSS's */
/* Font Families => font-family */
.ltx_font_serif { font-family: serif; }
.ltx_font_sansserif { font-family: sans-serif; }
.ltx_font_typewriter { font-family: monospace; }
/* dingbats should be converted to unicode? */
/* Math font families handled within math: script, symbol, fraktur, blackboard ? */
/* Font Series => font-weight */
.ltx_font_bold { font-weight: bold; }
.ltx_font_medium { font-weight: normal; }
/* Font Shapes => font-style or font-variant */
.ltx_font_italic { font-style: italic; font-variant:normal; }
.ltx_font_upright { font-style: normal; font-variant:normal; }
.ltx_font_slanted { font-style: oblique; font-variant:normal; }
.ltx_font_smallcaps { font-variant: small-caps; font-style:normal; }
.ltx_font_oldstyle { font-variant: oldstyle-nums; /* experimental css3 ? Doesn't seem to work!*/
font-style:normal;
-moz-font-feature-settings: "onum";
-ms-font-feature-settings: "onum";
-webkit-font-feature-settings: "onum";
font-variant-numeric: oldstyle-nums; }
.ltx_font_mathcaligraphic { font-family: "Lucida Calligraphy", "Zapf Chancery","URW Chancery L"; }
/*
.ltx_font_mathscript { ? }
*/
cite { font-style: normal; }
.ltx_red { color:red; }
/*.ltx_centering { text-align:center; margin:auto; }*/
/*.ltx_inline-block.ltx_centering,*/
/* Hmm.... is this right in general? */
.ltx_centering { display:block; margin:auto; text-align:center; }
/* Dubious stuff */
.ltx_hflipped {
display:inline-block;
-moz-transform: scaleX(-1);
-o-transform: scaleX(-1);
-webkit-transform: scaleX(-1);
transform: scaleX(-1);
filter: FlipH;
-ms-fliter: "FlipH"; }
.ltx_vflipped {
display:inline-block;
-moz-transform: scaleY(-1);
-o-transform: scaleY(-1);
-webkit-transform: scaleY(-1);
transform: scaleY(-1);
filter: FlipV;
-ms-fliter: "FlipV"; }
/* .ltx_phantom handled in xslt */
/* A somewhat amsart-like style */
.ltx_title_document {
text-transform:uppercase; font-size: 150%; text-align:center; }
.ltx_title_section {
font-family:serif; font-weight:normal; font-variant:small-caps; font-size: 120%;
text-align:center; }
.ltx_authors { text-align:center; }
.ltx_author_before,.ltx_author_after,
.ltx_role_author .ltx_personname {
text-transform:uppercase; }
.ltx_abstract { margin-bottom:1em; }
.ltx_title_abstract { font-family:serif; font-weight:normal; font-variant:small-caps;}
.ltx_title_keywords,
.ltx_title_classification { display:inline; font-size:100%; font-weight:normal; font-style:italic; }
.ltx_role_dedicatory { font-size:100%; font-style:italic; text-align:center; margin:1em;}
.ltx_figure .ltx_tag,
.ltx_table .ltx_tag {
font-family:serif; font-variant:small-caps;}
.ltx_title_subsection,
.ltx_title_subsubsection {
font-family:serif; font-weight:bold;}
.ltx_title_subsection:after,
.ltx_title_subsubsection:after,
.ltx_title_paragraph:after { content:"."; }
/* abstract, subsection, subsubsection and paragraph titles are runin */
.ltx_abstract .ltx_title,
.ltx_abstract .ltx_title + .ltx_p { display:inline; }
.ltx_subsection .ltx_title,
.ltx_subsection .ltx_title + .ltx_para,
.ltx_subsection .ltx_title + .ltx_para > .ltx_p { display:inline; }
.ltx_subsubsection .ltx_title,
.ltx_subsubsection .ltx_title + .ltx_para,
.ltx_subsubsection .ltx_title + .ltx_para > .ltx_p { display:inline; }
.ltx_paragraph .ltx_title,
.ltx_paragraph .ltx_title + .ltx_para,
.ltx_paragraph .ltx_title + .ltx_para > .ltx_p { display:inline; }
.ltx_title +.ltx_para > .ltx_p {text-indent:2em; }
.ltx_itemize,
.ltx_enumerate,
.ltx_description
{ margin-left: clamp(1em, 3%, 2em); }
.ltx_title_document { font-size:170%; text-align:center; margin:0.5em 0 0.5em 0; }
.ltx_authors,
.ltx_role_author { text-align:center; margin:0.5em 0 0.5em 0; }
.ltx_role_author .ltx_personname { font-size: 120%; }
.ltx_date { text-align:center; font-size: 120%; margin:0.5em 0 0.5em 0; }
.ltx_subtitle { text-align:center; font-size: 120%; padding-left:0.2em; margin-left:-0.5em; }
.ltx_title_abstract { text-align:center; font-size: 100%; font-weight:bold; }
.ltx_abstract { margin-left:4em; margin-right:4em; }
.ltx_title_acknowledgements,
.ltx_title_keywords,
.ltx_title_classification {
text-align:left; font-size: 100%; font-weight:bold; margin:0.5 0 0 0; }
.ltx_appendix,
.ltx_section,
.ltx_subsection,
.ltx_subsubsection { margin-top:1.5em; }
.ltx_paragraph,
.ltx_subparagraph { margin-top:1.0em; }
.ltx_title_appendix,
.ltx_title_section,
.ltx_title_bibliography { font-size:140%; font-weight:bold; margin-bottom:1em; }
.ltx_title_subsection { font-size:120%; font-weight:bold; margin-bottom:1em; }
.ltx_title_subsubsection { font-size:100%; font-weight:bold; margin-bottom:1em; }
/* Paragraph & Subparagraph titles should be runin! */
.ltx_title_paragraph { font-size:100%; font-weight:bold; display:inline;
margin-right:1em; }
.ltx_paragraph .ltx_title,
.ltx_paragraph .ltx_title + .ltx_para,
.ltx_paragraph .ltx_title + .ltx_para > .ltx_p { display:inline; }
.ltx_title_subparagraph { font-size:100%; font-weight:bold; display:inline;
margin-left:2em; margin-right:1em; }
.ltx_subparagraph .ltx_title,
.ltx_subparagraph .ltx_title + .ltx_para,
.ltx_subparagraph .ltx_title + .ltx_para > .ltx_p { display:inline; }
.ltx_figure { text-align:center; margin:auto; margin:0.5em; }
.ltx_table { text-align:center; margin:auto; margin:0.5em; }
/* experimental: flex model, to center short captions, justify long ones */
.ltx_table,
.ltx_figure { display:flex; flex-direction:column; align-items:center; }
.ltx_table .ltx_caption,
.ltx_figure .ltx_caption { text-align:justify; }
/* first p in para gets indented */
.ltx_para > .ltx_p:first-child { text-indent:2em; }
/* except the initial in a section */
section > .ltx_title +.ltx_para > .ltx_p,
section > .ltx_title +.ltx_date +.ltx_para > .ltx_p {text-indent:0em; }
.ltx_title_abstract + .ltx_p {text-indent:2em; }
.ltx_itemize,
.ltx_enumerate { margin-left:1em; }
.ltx_theorem { margin-top:0.5em; margin-bottom:0.5em; }
.ltx_theorem .ltx_title { margin-bottom:0.1em; }
.ltx_theorem .ltx_title + .ltx_para,
.ltx_theorem .ltx_title + .ltx_para .ltx_p,
.ltx_theorem .ltx_title + .ltx_p { margin-top:0em; }
Solving equations by interpolation
by
Ion Păvăloiu
1. Divided difference notation.
LetAndtwo normed linear spaces andan application. We considerdistinct elements of space
We designate bythe set of linear and continuous applications defined onand values inand we define by recurrence the following sets:
|
|
|
Or
In other words, the setsare the spaces of applications-linear and continuous defined onand values in
Definition 1 .
The applicationis called the first-order divided difference of the approximation
on the pointsof the system ( 1 ) if
-
1)
|
|
|
-
2)
if there exists the derivative in the Fréchet sense of the applicationon the pointSO
We have defined above the divided difference of order one of the applicationon two arbitrary elements of the system ( 1 ).
Now considering consecutive points of the system ( 1 ) we can consider the following divided differences:
|
|
|
Using the divided differences defined above we define the notion of divided difference of order of the application
We assume that we have defined the application
with the help of which we can define the divided differences of orderof the applicationonany elements of the system ( 1 ).
LetAndtwo divided differences of orderon the pointsrespectively
of the system ( 1 ).
Definition 2 .
The applicationis called divided difference of orderof the applicationon the pointsif
-
1)
|
|
|
|
|
|
-
2)
if there exists the derivative in the sense of Fréchet of order
of the applicationon the pointSO:
|
|
|
In order to be able to use these divided differences in the construction of interpolation polynomials, it is necessary to introduce the notion of divided difference symmetrical with respect to the points considered.
Definition 3 .
The Divided DifferenceOrare points of the system ( 1 ) is said to be symmetrical with respect to the points, if we have:
| (2) |
|
|
|
for each permutationnumbers
Ifthen, equality ( 2 ) is always verified. Then the divided differences of the real functions of a real variable are symmetric with respect to the points considered.
If we now assume for exampleAnd and if we designate by
two points ofWe have
| (3) |
|
|
|
Or
|
|
|
We show that the divided difference ( 3 ) satisfies conditions ( 1 ) and ( 2 ) of definition 1 , but that it is not symmetric with respect to the pointsAnd
But we can construct a divided difference on these points which is symmetric. For example if we consider the divided difference defined by:
|
|
|
|
|
|
|
|
so we obviously have
Be it nowwhere by
we designate the set of functions defined and contained on the intervalNow let us consider the application
given by the following equality:
|
|
|
Oris a continuous function over its entire domain of definition.
Be it now-functions
from spacewhich have the property that
Forand for each
The application's split differenceon functions
of the system ( 4 ) can be defined using the following equality:
|
|
|
Or
It is immediately shown that the linear application defined above has the properties ( 1 ) and ( 2 ) of definition 1 .
If we designate bythe divided difference of orderof the function, then the divided difference of orderof the applicationis defined by the following equality:
|
|
|
|
|
|
Or
2. Interpolation
In all that follows we will assume that the divided differences are symmetrical with respect to the points considered.
We now designate, byan application defined by the equality:
| (5) |
|
|
|
|
|
|
|
|
So we have the following theorem:
Theorem 1 .
The applicationdefined by ( 5 ) has the following properties
-
1)
-
2)
For
Demonstration..
We prove equality 2) by induction. Forequality 2) is obvious because
which is none other than equality 1) of definition 1. We assume that equality 2) is true byand we will show that it is also true for
Indeed, if the following equality takes place
| (6) |
|
|
|
|
|
|
|
|
|
|
|
|
then, taking into account equality
|
|
|
which results from the definition of the divided difference of order
The following equality results:
|
|
|
|
|
|
|
|
|
|
|
|
which shows that equality 2) holds for
It remains to demonstrate the equalities 1). From the definition of the applicationit results:
| (7) |
|
|
|
|
|
|
|
|
If we now consider the last term on the right of equality ( 7 ) and take into account the equality:
|
|
|
|
|
|
|
|
|
|
|
|
we get
|
|
|
|
|
|
|
|
Proceeding in the same way after a finite number of steps we arrive at the following equality:
|
|
|
The applicationis called the Lagrange interpolation polynomial of the applicationon the knotsof the system ( 1 ) ∎
3. Solving equations by interpolation
We consider the equation
is the zero element of the spaceAnd
is a nonlinear operator.
To solve equation ( 8 ) we can proceed as follows.
Assuming that the points of the system ( 1 ) are in a neighborhood of the solutionfrom equation ( 8 ), from equality 2) (theorem 1 ) the following equality results:
|
|
|
from which it results
|
|
|
If we now assume that we have the inequalities Or
is a real and positive number, then we have:
|
|
|
from which it follows that to find an approximation of the solution
you just have to solve the following equation:
The method outlined above has some drawbacks, including:
First, solving an equation in the form ( 9 ) even in the caseis a difficult mathematical problem and the difficulty of the problem increases considerably if
Even in the casethat is to say, whenis an algebraic polynomial with real coefficients, the problem can be very difficult because of the fact that an algebraic equation of degreehasroots, and then it is necessary to choose the root of equation ( 9 ) which approaches the root
of equation ( 8 ).
Second, it is well known that the position of the roots of an algebraic equation in the plane and the nature of these roots can change considerably when the coefficients of the equation are affected by non-essential errors.
In the following we will present another method by which some of the difficulties mentioned above can be eliminated.
We assume that linear operators
|
|
|
are reversible.
It follows that
Be it nowany set of elements of spaceWe designate bythe restriction of the applicationto the whole
We assume that the applicationis a homeomorphism of setsAndWhere bywe have designated the image of the wholeby
We have the following theorem:
Theorem 2 .
If the linear operator
Oris invertible, then the equality holds:
Or
Demonstration..
From the definition of divided difference of order 1 it follows that
|
|
|
Or
|
|
|
That's to say
|
|
|
from which results the equality of the statement of the theorem.
Let us now - elements of the set
We designate byapplication values
on these points, that is to say
We consider the application
given by
| (10) |
|
|
|
|
|
|
|
|
which verifies equalities
| (11) |
|
|
|
And
| (12) |
|
|
|
If we now assume that the solutionof equation ( 8 ) belongs to the setso obviously
and from ( 11 ) we deduce:
| (13) |
|
|
|
If we designate by
we deduce:
| (15) |
|
|
|
From equality ( 15 ) it follows that if the elements
are chosen in such a way thatOris a sufficiently small real number, thenis a good approximation for
For clarity we refer tothe operatorgiven by the relation ( 10 ). Now let
-given initial elements and either
the sequence generated by the recurrence relation:
| (16) |
|
|
|
We will now present some special cases of the iterative process ( 16 ).
A. LetAndtwo arbitrary elements of space
If we consider the first two terms of the inverse interpolation polynomial ( 10 ) taking into account Theorem
2 , we have
| (17) |
|
|
|
which is nothing other than the extension of the rope method.
B. In case we consider three initial elements
and we limit ourselves to the first three terms of the relation ( 10 ) we have:
| (18) |
|
|
|
|
|
|
|
|
|
|
|
|
This method is an extension to linear spaceof a method analogous to the well-known Tchebycheff method. From ( 15 ) and ( 16 ) we deduce
| (19) |
|
|
|
|
Regarding the method ( 16 ) the following two questions arise:
1. What are the conditions for the convergence of the method ( 16 ).
2. In the case of convergence of method ( 16 ), what is the speed of convergence?
Concerning the convergence speed of an iterative method of type ( 16 ), AM Ostrowski [ 1 ] showed that the order of convergence of this method does not exceed 2, for each given number of interpolation points.
In the following we will show that if at each iteration step we choose the interpolation points appropriately, then the order of convergence increases considerably.
As is well known for the resolution of an equation of the form ( 8 ) by an iterative method it is necessary to highlight an applicationwhich has the property that each solution of equation ( 8 ) is a fixed point for the application
An operatorwhich has the above property is said to be an iterative operator attached to equation ( 8 ).
For the operators attached to equation ( 8 ) we can define the notion of order.
Definition 4 .
Eithera set of space elements
Anda real number.
We say that the iterative operatorto order on the wholewith respect to equation ( 8 ) if for eachwe have
| (20) |
|
|
|
Let us now-iterative operators attached to equation ( 8 ) , respectively of order
We consider an arbitrary elementand designate bythe following expressions:
| (21) |
|
|
|
If the initial elements in the iterative method ( 16 ) are:
so if we writewe have:
| (23) |
|
|
|
We now assume that we have built the first
elements of the suitethen the elementis obtained in the following way:
We write:
| (24) |
|
|
|
And
|
|
|
So the elementhas the following form:
| (25) |
|
|
|
From ( 15 ) and ( 25 ) we have:
| (26) |
|
|
|
But since we assumed that iterative operators
have the orders respectivelywith the constantsthe following relationships result:
And
| (28) |
|
|
|
|
|
|
|
|
|
|
|
|
For
If we write now
| (29) |
|
|
|
And
| (31) |
|
|
|
and if we assume that the divided differencesare strongly bounded by the same constantthen we deduce from ( 26 );
| (32) |
|
|
|
If we now assume that the first-order divided differences of the applicationare strongly bounded by the same constant
we have:
| (33) |
|
|
|
Taking into account inequalities ( 33 ) and ( 32 ) we obtain:
| (34) |
|
|
|
By multiplying the inequalities ( 34 ) byand writing
|
|
|
we obtain from ( 34 )
We will now admit that we can choose the constants
and the initial pointin such a way that
From inequalities ( 35 ) we deduce:
hence taking into account the conditionwe deduce
That's to say
We now consider some special cases of the iterative method ( 25 ).
1. If we consider a single iterative operatorof order
then the iterative method ( 25 ) returns to the well-known Aitken-Steffensen method. In this case we have:
2. IfOr
SOAnd
3. If AndSO
And
Bibliography
-
[1]
AM Ostrowski, Reşenie uravnenii i sistem uravnenii, Izd. inostr. bed. Moskva (1963).
-
[2]
I. Păvăloiu,
††margin:
clickable
Interpolation in normed linear spaces and applications , Mathematica, Cluj, nr. 12 (35) , 2, (1970), pp. 309–324.
-
[3]
I. Păvăloiu,
††margin:
clickable
Introducere în teoria aproximării soluţiilor ecuaţiilor , Editura Dacia, 1976.
-
[4]
T. Popoviciu,
††margin:
clickable
Introduction to the theory of divided differences , Mathematical Bulletin of the Romanian Society of Sciences, 42 , 1, (1940) pp. 65–78. JSTOR
-
[5]
AS Sergeev, O metode hord, Sibirski mat. Jurnal, XI, (2), (1961), pp. 282–289.
-
[6]
JF Traub, Iterative Methods for the Solution of Equations, Prentice Hall, Series in Automatic Computation, 1964.
-
[7]
S. Ul'm, Ob obobscennyh razdelennîh raznostiah , II, Izv. Nauk Estonskoi SSR, 16 , 2, (1967), 146–155.