Bilateral approximations of the roots of scalar equations by Lagrange-Aitken-Steffensen method of order three

Abstract

We study the monotone convergence of two general methods of Aitken-Steffenssen type. These methods are obtained from the Lagrange inverse interpolation polynomial of degree two, having controlled nodes. The obtained results provide information on controlling the errors at each iteration step.

Author

Ion Păvăloiu
(Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy)

Keywords

Aitken-Steffenssen methods; Lagrange inverse interpolation

PDF

PDF-LaTeX file (on the journal website).

Cite this paper as:

I. Păvăloiu, Bilateral approximations of the roots of scalar equations by Lagrange-Aitken-Steffensen method of order three, Rev. Anal. Numér. Théor. Approx., 35 (2006) no. 2, pp. 173-182. https://doi.org/10.33993/jnaat352-843

About this paper

Print ISSN

1222-9024

Online ISSN

2457-8126

References

[1] Balazs, M., A bilateral approximating method for finding the real roots of real equations, Rev. Anal. Numer. Theor. Approx., 21 (2), pp. 111-117, 1992.

[2] Casulli, V., Trigiante, D.,  The convergence order for iterative multipoint procedures, Calcolo, 13 (1), pp. 25-44, 1997.

[3] Costabile, F., Gualtieri, I.M., Luceri, R., A new iterative method for the computation of the solution of nonlinear equations, Numer. Algorithms, 28, pp. 87-100, 2001.

[4] Frontini, M., Hermite interpolation and a new iterative method for the computation of the roots of non-linear equations, Calcolo, 40, pp. 109-119, 2003.

[5] Grau, M., An improvement to the computing of nonlinear equation solutions, Numer. Algorithms., 34, pp. 1-12, 2003.

[6] Ostrowski, A., Solution of Equations in Euclidian and Banach Spaces, Academic Press, New York and London, 1973.

[7] Păvăloiu I., Optimal efficiency index for iterative methods of interpolatory type, Computer Science Journal of Moldova, 1 (5), pp. 20-43, 1997.

[8] Păvăloiu I.,  Approximation of the roots of equations by Aitken-Steffensen-type monotonic sequences, Calcolo, 32 (1-2), pp. 69-82, 1995.

[9] Păvăloiu I., Optimal problems concerning interpolation methods of solution of equations, Publications de L’Institut Mathematique, 52 (66), pp. 113-126, 1992.

[10] Păvăloiu I., Optimal effiency index of a class of Hermite iterative methods, with two steps, Rev. Anal. Numer. Theor. Approx., 29 (1), pp. 83-89, 2000.

[11] Păvăloiu I., Local convergence of general Steffensen type methods, Rev. Anal. Numer. Theor. Approx., 33 (1), pp. 79-86, 2004.

[12] Păvăloiu I.,  and Pop, N., Interpolation and Applications, Risoprint, Cluj-Napoca, 2005 (in Romanian).

[13] Păvăloiu I., On a Steffensen-Hermite-type Method for approximating the solution of nonlinear equations, Rev. Anal. Numer. Theor. Approx., 25 1, pp. 87-94, 2006.

[14] Păvăloiu I.,  Bilateral approximation of solutions of equations by order-three Steffensen type methods, Studia Univ. “Babeș-Bolyai”, Mathematica, Vol. LI, no. 3, pp. 105-114, 2006.

[15] Traub, J.F., Iterative Methods for Solutions of Equations, Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1964.

[16] Turowicz, B. A., Sur les derivees d’ordre superieur d’une function inverse, Ann. Polon. Math., 8, pp. 265-269, 1960.

Paper (preprint) in HTML form

/*====================================================================== 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; } BILATERAL APPROXIMATIONS OF THE ROOTS OF SCALAR EQUATIONS BY LAGRANGE-AITKEN-STEFFENSEN METHOD OF ORDER THREE∗

BILATERAL APPROXIMATIONS OF THE ROOTS OF SCALAR EQUATIONS BY LAGRANGE-AITKEN-STEFFENSEN METHOD
OF ORDER THREE∗

Ion Păvăloiu
(Date: January 11, 2006.)
Abstract.

We study the monotone convergence of two general methods of Aitken-Steffenssen type. These methods are obtained from the Lagrange inverse interpolation polynomial of degree two, having controlled nodes. The obtained results provide information on controlling the errors at each iteration step.

Key words and phrases:
Aitken-Steffenssen methods, Lagrange inverse interpolation.
1991 Mathematics Subject Classification:
65H05.
∗This work has been supported by MEdC-ANCS under grant 2CEEX-06-11-96.
†“Tiberiu Popoviciu” Institute of Numerical Analysis, P.O. Box. 68-1, Cluj-Napoca,Romania, e-mail: pavaloiu@ictp.acad.ro.

1. Introduction

It is well known that the Steffensen, Aitken, and Aitken-Steffensen methods are obtained from the inverse Lagrange interpolation polynomial of degree one, with controlled nodes [9]–[13], [7]. Consider the equation

(1.1) f⁢(x)=0

where

f:[a,b]→R,a,b∈R,a<b.

We also consider the following three equations, each of them equivalent with equation (1.1):

(1.2) x=g⁢(x),g:[a,b]→[a,b]

and

(1.3) x = g1⁢(x),g1:[a,b]→[a,b],
x = g2⁢(x),g2:[a,b]→[a,b].

The Steffensen method in given by relations

(1.4) xn+1=xn−f⁢(xn)[xn,g⁢(xn);f],x0∈[a,b],n=0,1,2,…,

and, analogously, the Aitken method is of the following form:

(1.5) xn+1=g1⁢(xn)−f⁢(g1⁢(xn))[g1⁢(xn),g2⁢(xn);f],x0∈[a,b],n=0,1,2,….

Finally, the Aitken-Steffensen method is given by the relations:

(1.6) xn+1=g1⁢(xn)−f⁢(g1⁢(xn))[g1⁢(xn),g2⁢(g1⁢(xn));f],x0∈[a,b],n=0,1,2,….

The order of convergence of all three methods (1.4)–(1.6) is at least two, and this order depends on the functions g and g1, g2 respectively. Essentially, the methods (1.4)–(1.6) are obtained from the method of chord where the interpolation nodes depend on the functions g and respectively g1,g2. In papers [2], [9], [10], [13] some conditions had to be considered in order that all the three methods (1.4)–(1.6) generate two sequences (un)n≥0 and (vn)n≥0 with the properties:

α) The sequence (un)n≥0 is increasing and the sequence (vn)n≥1 is decreasing:

β)limun=limvn=x¯, where x¯ is the root of equation (1.1), x¯ ∈[a,b].

Practically, such sequences are very interesting, because by inequalities

max⁡{x¯−un,vn−x¯}≤vn−un,n=0,1,2,…,

the errors of approximation at every step of iteration may be controlled.

Let a1,a2, a3∈[a,b] be three nodes of interpolation, and let b1, b2, b3 the values of the function f, i.e.

b1=f⁢(a1),b2=f⁢(a2),b3=f⁢(a3).

Suppose that the function f:[a,b]→F, is bijective where F=f⁢([a,b]).

Then there exists f−1:F→[a,b] and the following equality holds [13], [15]:

f−1⁢(y) = a1+[b1,b2;f−1]⁢(y−b1)
+[b1,b2,b3;f−1]⁢(y−b1)⁢(y−b2)
+[y,b1,b2,b3;f−1]⁢(y−b1)⁢(y−b2)⁢(y−b3),

for every y∈F.

If x¯∈[a,b] is the root of equation (1.1), then x¯=f−1⁢(0), and by (1) one obtains the following representation of x¯:

x¯ = a1−[b1,b2;f−1]⁢b1+[b1,b2,b3;f−1]⁢b1⁢b2
−[0,b1,b2,b3;f−1]⁢b1⁢b2⁢b3.

By relations (see [13])

[b1,b2;f−1]=1[a1,a2;f]

and

[b1,b2,b3;f−1]=−[a1,a2,a3;f][a1,a2;f]⁢[a1,a3;f]⁢[a2,a3;f],

using (1), one obtains the following approximation for x¯:

(1.9) a4=a1−f⁢(a1)[a1,a2;f]−[a1,a2,a3;f]⁢f⁢(a1)⁢f⁢(a2)[a1,a2;f]⁢[a1,a3;f]⁢[a2,a3;f]

and the error:

(1.10) x¯−a4=−[0,b1,b2,b3;f−1]⁢f⁢(a1)⁢f⁢(a2)⁢f⁢(a3).

If f∈C3⁢([a,b]) and f′⁢(x)≠0, ∀x∈[a,b], then f−1∈C3⁢(F) and the following equality holds (see [13], [17]):

(1.11) [f−1⁢(y)]′′′=3⁢[f′′⁢(x)]2−f′⁢(x)⁢f′′′⁢(x)[f′⁢(x)]5,

where y=f⁢(x).

Using the mean value formula for divided differences (see [13]) it follows that there exists η∈F such that

(1.12) [0,b1,b2,b3;f−1]=(f−1)′′′⁢(η)6.

Because f is bijective, it follows that there exists ξ∈[a,b] such that η=f⁢(ξ), and by (1.11) and (1.12) one obtains:

(1.13) [0,b1,b2,b3;f−1]=3⁢[f′′⁢(ξ)]2−f′⁢(ξ)⁢f′′′⁢(ξ)6⁢[f′⁢(ξ)]5.

By (1.9), if one considers particular nodes a1,a2,a3 it is possible to obtain different methods of Steffensen type, of Aitken type or of Aitken-Steffensen type.

Let xn∈[a,b] be an approximation of the root x¯ of equation (1.1).

If one considers a1=xn, a2=g⁢(xn), a3=g⁢(g⁢(xn)), then it follows the following method of Steffensen type [9], [10], [13]:

xn+1 = xn−f⁢(xn)[xn,g⁢(xn);f]
−[xn,g⁢(xn),g⁢(g⁢(xn));f]⁢f⁢(xn)⁢f⁢(g⁢(xn))[xn,g⁢(xn);f]⁢[xn,g⁢(g⁢(xn));f]⁢[g⁢(xn),g⁢(g⁢(xn));f],

x0∈[a,b],n=0,1,2,….

If a1=xn, a2=g1⁢(xn), a3=g2⁢(g1⁢(xn)), then one obtains the following method of Aitken-Steffensen type:

xn+1 = xn−f⁢(xn)(xn,g1⁢(xn);f]
−[xn,g1⁢(xn),g2⁢(g1⁢(xn));f]⁢f⁢(xn)⁢f⁢(g1⁢(xn))[xn,g1⁢(xn);f]⁢[xn,g2⁢(g1⁢(xn));f]⁢[g1⁢(xn),g2⁢(g1⁢(xn));f],

n=0,1,…, x0∈[a,b], and finally, for a1=xn,a2=g1⁢(xn), a3=g2⁢(xn) one obtains the following method of Aitken type:

xn+1 = xn−f⁢(xn)[xn,g1⁢(xn);f]
−[xn,g1⁢(xn),g2⁢(xn);f]⁢f⁢(xn)⁢f⁢(g1⁢(xn))[xn,g1⁢(xn);f]⁢[xn,g2⁢(xn);f]⁢[g1⁢(xn),g2⁢(xn);f],
n = 0,1,…,x0∈[a,b].

Using the symmetry of Lagrange polynomial with respect to nodes, by permutations of a1, a2, a3 in methods (1)–(1), one obtains the same results for xn+1.

In [15] conditions are given in order that method (1) generates sequences approximating the root of equation (1.1) bilaterally. In this paper we study methods (1) and (1) and we obtain same conditions in order that the sequences generated by there methods are bilateral approximations of the root x¯ of equations (1.1).

2. The convergence of Aitken-Steffensen method

In the following we study the method (1) and we search the conditions on the sequences (xn)n≥0, (g1⁢(xn))n≥0 and (g2(g1(xn))n≥0 generated from this method in order that they are monotonic sequences, bilaterally approximating the root x¯ of equation (1.1).

We need the following hypothesis:

  1. a)

    g1 is increasing on [a,b];

  2. b)

    g2 is continuous and decreasing on [a,b];

  3. c)

    equation (1.1) has a solution x¯∈[a,b] and g1⁢(x¯)=g2⁢(x¯)=x¯,

  4. d)

    function f is in C3⁢([a,b]), and for every x∈[a,b] the following relation is fulfilled:

    (2.1) 3⁢[f′′⁢(x)]2−f′⁢(x)⁢f′′′⁢(x)<0;
  5. e)

    function g1 satisfies inequality

    |g1⁢(x)−g1⁢(x¯)|≤L⁢|x−x¯|,∀x∈[a,b],

    where 0<L<1.

Concerning the convergence of sequence (xn)n≥0 generated by (1), the following Theorem holds:

Theorem 2.1.

Let x0∈[a,b] and f,g1,g2 verify the following conditions:

  1. i)1

    f is increasing on [a,b];

  2. ii)1

    f is convex on [a,b];

  3. iii)1

    functions g1,g2 and f verify hypotheses a)−e);

  4. iv)1

    x0>x¯ and g2⁢(g1⁢(x0))≥a.

Then the sequences (xn)n≥0, (g1⁢(xn))n≥0 and (g2(g1(xn))n≥0 generated by (1) verify the properties:

  1. j)1

    sequences (xn)n≥0 and (g1⁢(xn))n≥0 are decreasing and bounded from below by x¯ ;

  2. jj)1

    sequence (g2(g1(xn))n≥0 is increasing and bounded from above by x¯;

  3. jjj)1

    at every iteration step the following inequalities hold:

    (2.2) xn−x¯≤xn−g2⁢(g1⁢(xn)),n=0,1,…;
  4. jv)1

    limxn=limg1(xn)=limg2(g1(xn)=x¯.

Proof.

By hypotheses a) and x0>x¯ it follows g1⁢(x0)>g1⁢(x¯), i.e g1⁢(x0)>x¯. Hypothesis e) implies g1⁢(x0)<x0. From 𝐢𝐯1) it follows g2(g1(x0)≥a and from b) and g1⁢(x0)>x¯ one obtains g2⁢(g1⁢(x1))<x¯. Now using 𝐢1) it follows:

f⁢(x0)>0,f⁢(g1⁢(x0))>0

and

f(g2(g1(x0))<0,

and by considering 𝐢1) and 𝐢𝐢1) for n=0 in (1) it follows that x1<x0. Using (1.13), by hypotheses d) and (1.10) for a4=x1 and b1=f(x0),b2=f(g1(x0)),b3=f(g2(g1(x0)) it follows x¯−x1<0, i.e x1>x¯. By hypotheses a) and x1<x0 it follows g1⁢(x1)<g1⁢(x0) and then, by b) one obtains g2⁢(g1⁢(x1))≥g2⁢(g1⁢(x0)), and g2⁢(g1⁢(x1))<x¯. Let xm, m∈ ℕ be an element of sequence (xn)n≥0 generated by (1) and suppose that xm>x¯. Then one obtains the relations:

(2.3) g2⁢(g1⁢(xm))<g2⁢(g1⁢(xm+1))<x¯<g1⁢(xm+1)<xm+1<g1⁢(xm)<xm.

Inequality xm+1<g1⁢(xm) follows from equality:

xm+1 = xm−f⁢(xm)[xm,g1⁢(xm);f]
−[xm,g1(xm),g2(g1(xm));f]f(xm)f(g1(xm)[xm,g1⁢(xm);f]⁢[xm,g2⁢(g1⁢(xm));f]⁢[g1⁢(xm),g2⁢(g1⁢(xm));f]
= g1⁢(xm)−f(g1(xm)[xm,g1⁢(xm);f]
−[xm,g1⁢(xm),g2⁢(g1⁢(xm));f]⁢f⁢(g1⁢(xm))⁢f⁢(xm)[xm,g1⁢(xm);f]⁢[xm,g2⁢(g1⁢(xm));f]⁢[g1⁢(xm),g2⁢(g1⁢(xm));f]

and from the hypothesis of the theorem.

From relations (2.3) it follows (2.2) and conclusions 𝐣1) and 𝐣𝐣1). Conclusion 𝐣𝐯1) is obvious. The theorem is proved. ∎

Remark 2.2.

If function f is concave and decreasing on [a,b], then function h:[a,b]→ℝ defined by h⁢(x)=−f⁢(x) in convex and increasing. The relation (2.1) is also verified for h. Consequently the sequence generated by (1) for function h, verifies all the conditions of Theorem 2.1, and then the conclusions of this theorem hold.

An analogous proof with that from Theorem 2.1 is valid for the following:

Theorem 2.3.

Let x0∈[a,b] and f,g1,g2 verify the following conditions:

  1. i)2

    function f is increasing on [a,b];

  2. ii)2

    function f is concave on [a,b];

  3. iii)2

    functions g1,g2 and f verify the hypotheses a)−e);

  4. iv)2

    x0<x¯ and g2(g1(x0)≤b.

Then sequences (xn)n≥0, (g1⁢(xn))n≥0 and (g2⁢(g1⁢(xn≥0))) generated by (1) have the following properties:

  1. j)2

    sequences (xn)n≥0 and (g1⁢(xn))n≥0 are increasing and bounded from above by x¯;

  2. jj)2

    sequence g2⁢(g1⁢(xn)) is decreasing and bounded from below by x¯;

  3. jjj)2

    the following relations hold:

    (2.4) x¯−xn≤g2⁢(g1⁢(xn))−xn,n=0,1,2,…,
  4. jv)2

    limxn=limg1⁢(xn)=limg2⁢(g1⁢(xn))=x¯.

Remark 2.4.

If function f is decreasing and convex, then h1:[a,b]→ℝ defined by h1⁢(x)=−f⁢(x) is increasing and concave, h1verifies all hypotheses of Theorem 2.3 and consequently all the conclusions 𝐣2)−𝐣𝐯2) hold.

3. Convergence of Aitken type method

In order to study the convergence of the sequences generated by (1) we must suppose that functions f,g1 and g2 verify hypotheses a)-e) of section 2.

The following theorem holds:

Theorem 3.1.

Let x0∈[a,b] and the functions f,g1 and g2 verify the conditions:

  1. i)3

    f is increasing an [a,b];

  2. ii)3

    f is convex an [a,b];

  3. iii)3

    f, g1 and g2 verifies the hypotheses a)−e);

  4. iv)3

    x0>x¯ and g2⁢(x0)>a.

Then sequences (xn)n≥0, (g1⁢(xn))n≥0 and (g2⁢(xn))n≥0 generated by (1) have the properties:

  1. j)3

    sequences (xn)n≥0 and (g1⁢(xn))n≥0 are decreasing and bounded from below by x¯ ;

  2. jj)3

    sequence (g2⁢(xn))n≥0 is increasing and bounded from above by x¯;

  3. jjj)3

    the following relations hold,

    xn−x¯≤xn−g2⁢(xn),n=0,1,…;
  4. jv)3

    limxn=limg1⁢(xn)=limg2⁢(xn)=x¯.

Proof.

Let be xm∈[a,b], xm>x¯ and g2⁢(xm)>a, where m∈ℕ. By a) it follows that g1⁢(xm)>x¯ and by b) g2⁢(xm)<x¯. Using e) one obtains g1⁢(xm)<xm. By hypothesis 𝐢3) and 𝐢𝐢3), and by the above relations and (1), for n=m it follows xm+1<xm.

For a1=xm, a2=g1⁢(xm) and a3=g2⁢(xm) in (1.10), by (1.13) and hypothesis d) it follows that xm+1>x¯. Observe that xm+1 in (1) may by expressed in the following way:

xm+1 = g1⁢(xm)−f⁢(g1⁢(xm))[g1⁢(xm),xm;f]−
−[xm,g1⁢(xm),g2⁢(xm);f]⁢f⁢(g1⁢(xm))⁢f⁢(xm)[xm,g1⁢(xm);f]⁢[xm,g2⁢(xm);f]⁢[g1⁢(xm),g2⁢(xm);f]

and then xm+1<g1⁢(xm). By relation xm+1<xm it follows that g2⁢(xm+1)>g2⁢(xm). Consequently it follows that for every m∈ℕ, the following relations hold:

g2⁢(xm)<g2⁢(xm+1)<x¯<g1⁢(xm+1)<xm+1<g1⁢(xm)<xm.

By these relations conclusions 𝐣𝐣𝐣3) and 𝐣𝐯3) also follow. ∎

Remark 3.2.

If function f is decreasing and concave, then function h⁢(x)=−f⁢(x), x∈[a,b] verifies all the hypothesis of Theorem 3.1 and consequently the sequence (xn)n≥1 generated by (1) verifies all the conclusions in Theorem 3.1.

Analogously, the following result can be proved

Theorem 3.3.

Let x0∈[a,b] and functions f,g1,g2 have the following properties:

  1. i4)

    f is increasing on [a,b];

  2. ii)4

    f is concave on [a,b];

  3. iii)4

    f,g1, and g2 verify hypothesis a)-e);

  4. iv)4

    x0<x¯ and g2⁢(x0)≤b.

Then sequences (xn)n≥0, (g1⁢(xn))n≥0 and (g2⁢(xn))n≥0 generated by (1) verify the properties:

  1. j)4

    (xn)n≥0 and (g1⁢(xn))n≥0 are increasing and bounded from above by x¯;

  2. jj)4

    sequence (g2⁢(xn))n≥0 is decreasing and bounded from below by x¯;

  3. jjj)4

    the following relations hold:

    x¯−xn≤g2⁢(xn)−xn,n=0,1,2,…,
  4. jv)4

    limxn=limg1⁢(xn)=limg2⁢(xn)=x¯.

Remark 3.4.

If function f is decreasing and convex, then function h⁢(x) =−f⁢(x), x∈[a,b] is increasing and concave. It follows that sequence (xn)n≥0 generated by (1) verifies the hypotheses and all the conclusions of Theorem 3.3.

4. The determination of the auxiliary functions

In the following, for every situation concerning the monotonicity and convexity of function f, functions g1 and g2 can be determined such that conditions a), b), c) and e) should be verified. In the following we thoroughly present the case in which function f is increasing and convex.

Supposing that f′⁢(x)>0 for every x∈[a,b], one considers the functions:

(4.1) g1⁢(x)=x−f⁢(x)f′⁢(b),

and

(4.2) g2⁢(x)=x−f⁢(x)f′⁢(a).

Then

g1′⁢(x)=1−f′⁢(x)f′⁢(b)≥0

for every x∈[a,b] and, consequently g1 is increasing. Analogously

g2′⁢(x)=1−f′⁢(x)f′⁢(a)≤0

for every x∈[a,b] and then g2 is a decreasing function. It follows that g1 and g2 verify hypotheses a) and b). The hypothesis e) is also verified because g1′⁢(x)<1, for every x∈[a,b]. In iv1) the inequality g2(g1(x0)≥a holds. Because g1⁢(x0)<x0 and x0>x¯, the above inequality is verified if g2⁢(x0)≥a. This follows, because g2⁢(g1⁢(x0))>g2⁢(x0) (g2 is a decreasing functions). This means that the following relation must hold:

x0−f⁢(x0)f′⁢(a)≥a

i.e.

x0≥a+f(x0))f′⁢(a).

Because f⁢(x0)>0 and f′⁢(a)>0, for the veridicity of last inequality it is sufficient that

a+f⁢(x0)f′⁢(a)≤x¯,

where x¯ is the root of equation (1.1). This last inequality may be realized if x0 is sufficiently close to x¯.

Consequently, the hypotheses of theorems 2.1 and 3.1 are realized for g1 and g2 considered above. The other cases may be similarly analyzed.

5. Order of convergence

In the following we prove that every method (1) and (1) have the order of convergence three. The order of convergence for method (1) was treated in [15], and this order is at least three. Assume the following hypotheses:

α) function g2 verifies relation

|g2⁢(x)−g2⁢(x¯)|≤p⁢|x−x¯|,

for every x∈[a,b], where p>0, p∈ℝ;

β) m≤|f′⁢(x)|≤M, for every x∈[a,b], where m>0 and M>0 are real numbers.

γ)|3[f′(x)]2−f′(x)f′′′(x)|≤q, for every x∈[a,b], where q>0, q∈ℝ .

In hypotheses α),β),γ) and e), by (1.10) and (1.13), for sequence (xn)n≥0 generated by (1), one obtains:

|x¯−xn+1|≤q⁢M3⁢L2⁢p6⁢m5⁢|x¯−xn|3,n=0,1,2,…

and this means that the order of convergence for method (1) is at least three.

With the same hypotheses, for tree sequence (xn)n≥0 generated by (1) it follows:

|x¯−xn+1|≤q⁢M3⁢L⁢p6⁢m5⁢|x¯−xn|3,n=0,1,2,…

i.e. the order of convergence of (1) is also at least three.

References

  • [1]
  • [2] Balázs, M., A bilateral approximating method for finding the real roots of real equations, Rev. Anal. Numér. Théor. Approx., 21 (2), pp. 111–117, 1992.
  • [3] Casulli, V., Trigiante, D., The convergence order for iterative multipoint procedures, Calcolo, 13 (1), pp. 25–44, 1997.
  • [4] Costabile, F., Gualtieri, I. M., Luceri, R., A new iterative method for the computation of the solution of nonlinear equations, Numer. Algorithms, 28, pp. 87–100, 2001.
  • [5] Frontini, M., Hermite interpolation and a new iterative method for the computation of the roots of non-linear equations, Calcolo, 40, pp. 109–119, 2003.
  • [6] Grau, M., An improvement to the computing of nonlinear equation solutions, Numer. Algorithms., 34, pp. 1–12, 2003.
  • [7] Ostrowski, A., Solution of Equations in Euclidian and Banach Spaces, Academic Press, New York and London, 1973.
  • [8] Păvăloiu, I., Optimal efficiency index for iterative methods of interpolatory type, Computer Science Journal of Moldova, 1 (5), pp. 20–43, 1997.
  • [9] Păvăloiu, I., Approximation of the roots of equations by Aitken-Steffensen-type monotonic sequences, Calcolo, 32 (1–2), pp. 69–82, 1995.
  • [10] Păvăloiu, I., Optimal problems concerning interpolation methods of solution of equations, Publications de L’Institut Mathématique, 52 (66), pp. 113–126, 1992.
  • [11] Păvăloiu, I., Optimal effiency index of a class of Hermite iterative methods, with two steps, Rev. Anal. Numér. Théor. Approx., 29 (1), pp. 83–89, 2000.
  • [12] Păvăloiu, I., Local convergence of general Steffensen type methods, Rev. Anal. Numér. Théor. Approx., 33 (1), pp. 79–86, 2004.
  • [13] Păvăloiu, I. and Pop, N., Interpolation and Applications, Risoprint, Cluj-Napoca, 2005 (in Romanian).
  • [14] Păvăloiu, I., On a Steffensen-Hermite-type Method for approximating the solution of nonlinear equations, Rev. Anal. Numér. Théor. Approx., 25 1, pp. 87–94, 2006.
  • [15] Păvăloiu, I., Bilateral approximation of solutions of equations by order-three Steffensen type methods, Studia Univ. “Babeş-Bolyai”, Mathematica, Vol. LI, no. 3, pp. 105–114, 2006.
  • [16] Traub, J. F., Iterative Methods for Solutions of Equations, Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1964.
  • [17] Turowicz, B. A., Sur les derivées d’ordre supérieur d’une function inverse, Ann. Polon. Math., 8, pp. 265–269, 1960.
  • [18]
2006

Related Posts