On some Aitken-Steffensen-Halley-type method for approximating the roots of scalar equations

Abstract

We extend the Aitken-Steffensen method to the Halley transformation. Under some rather simple assumptions we obtain error bounds for each iteration step; moreover, the convergence order of the iterates is 3, i.e. higher than for the Aitken-Steffensen case.

Author

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

Keywords

nonlinear equations in R; Aitken-Steffensen-Halley method.

PDF

PDF-LaTeX file on the journal website.

Cite this paper as:

I. Păvăloiu, On some Aitken-Steffensen-Halley-type method for approximating the roots of scalar equations, Rev. Anal. Numér. Théor. Approx., 30 (2001) no. 2, pp.,  207-212, https://ictp.acad.ro/jnaat/journal/article/view/2001-vol30-no2-art8

About this paper

Print ISSN

1222-9024

Online ISSN

2457-8126

References

[1] Argyros, I.K. On the convergence of some projection methods with perturbations, J. Comp. Appl. Math. 36, (1991), 255–258.

[2] Argyros, I.K. On an application of the Zincenko method to the approximation of implicit functions, Z.A.A. 10, 3, (1991), 391– 396.

[3] Argyros, I.K. and Szidarovszky, F. The Theory and Application of Iteration Methods, C.R.C. Press, Inc., Boca Raton, Florida, 1993.

[4] Catinas, E. On some iterative methods for solving nonlinear equations, Revue d’analyse Numerique et de theorie de l’approximation, 23, 1, (1994), 47–53.

[5] Kantorovich, L.V. The method of successive approximation for functional equations, Acta Math. 71 (1939), 63–97.

[6] Pavaloiu, I. Sur une generalisation de la methode de Steffensen, Revue d’analyse Numerique et de theorie de l’approximation, 21, 1, (1992), 59–65.

[7] Pavaloiu, I. Bilateral approximations for the solutions of scalar equations, Revue d’analyse numerique et de theorie de l’approximation, 23, 1, (1994), 95–100.

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; } ON SOME AITKEN-STEFFENSEN-HALLEY-TYPE METHODS FOR APPROXIMATING THE ROOTS OF SCALAR EQUATIONS

ON SOME AITKEN-STEFFENSEN-HALLEY-TYPE METHODS
FOR APPROXIMATING THE ROOTS OF SCALAR EQUATIONS

ION PĂVĂLOIU∗
Abstract.

In this note we extend the Aitken-Steffensen method to the Halley transformation. Under some rather simple assumptions we obtain error bounds for each iteration step; moreover, the convergence order of the iterates is 3, i.e. higher than for the Aitken-Steffensen case.

1991 Mathematics Subject Classification:
65H05
∗”T. Popoviciu” Institute of Numerical Analysis, P.O. Box 68-1, 3400 Cluj–Napoca, Romania, e-mail: pavaloiu@ictp-acad.math.ubbcluj.ro.

1. Introduction

Let f:[a,b]→ℝ, a,b∈ℝ, a<b and suppose that f∈C4⁢[a,b], and f′⁢(x)>0, ∀x∈[a,b]. Consider the function h:[a,b]→ℝ given by

h⁢(x)=f⁢(x)f′⁢(x).

As it was shown in [2], the Halley method for solving:

(1) f⁢(x)=0,

is given by

(2) xn+1=xn−h⁢(xn)h′⁢(xn),n=0,1,…,x0∈[a,b].

This sequence is in fact generated by the Newton method for solving h⁢(x)=0.

The first and second order derivatives of h are given by

(3) h′⁢(x) =2⁢(f′⁢(x))2−f′′⁢(x)⋅f⁢(x)2⁢(f′⁢(x))3/2
(4) h′′⁢(x) =3⁢(f′′⁢(x))2−2⁢f′′′⁢(x)⁢f′⁢(x)4⁢(f′⁢(x))5/2⋅f⁢(x)

which yield the following equalities for a solution x¯ of (1):

(5) h′⁢(x¯) =(f′⁢(x¯))1/2,and
(6) h′′⁢(x¯) =0.

Relation (6) ensures the convergence order 3 for the sequence (xk)k≥0.

In the papers [2]–[8] and [12] there are studied the convergence and the convergence order of some sequences generated by some interpolatory methods applied to equation h⁢(x)=0.

We shall consider other two equations equivalent to (1) of the form:

(7) x−φ1⁢(x) =0,and
(8) x−φ2⁢(x) =0

The Aitken method for solving h⁢(x)=0 is given by the iteration

(9) xn+1=φ1⁢(xn)−h⁢(φ1⁢(xn))[φ1⁢(xn),φ2⁢(xn);h],n=0,1,…,x0∈[a,b].

In this note we shall study the convergence of these iterates. We shall show that the functions φ1 and φ2 may be chosen in order to obtain bilateral approximations at each iteration step; this fact allows the control of the errors. On the other hand, the convergence order of (xn)n≥0 given by (9) is at least equal to 3.

Hypotheses f∈C4⁢[a,b] and f′⁢(x)>0, ∀x∈[a,b] imply, taking into account (5), that there exist α,β∈ℝ, a≤α<x¯<β≤b such that h′⁢(x)>0,∀x∈[α,β].

2. ERROR EVALUATION AND LOCAL CONVERGENCE

Consider the interval [α,β] given above. We shall make the following assumptions on φ1 and φ2:

  • i.

    the function f∈C4⁢[a,b];

  • ii.

    equation (1) has the solution x¯∈[a,b];

  • iii.

    the inequality f′⁢(x)>0 holds for x∈[α,β];

  • iv.

    the function φ1verifies the relation 0<[x,y;φ1]<1 for all x,y∈[α,β], where [x,y;φ] denotes the first order divided difference of φ on x and y;

  • v.

    the function φ2 verifies the relations −1<[x,y;φ2]<0 for all x,y∈[α,β].

We can state the following result:

Theorem 2.1.

Assume that i–v hold, and for some x0∈[α,β] sufficiently close to x¯ we have φ1⁢(x0),φ2⁢(x0)∈[α,β]. Then the following relations hold:

  • j.

    the sequences (xn)n≥0,(φ1⁢(xn))n≥0 and (φ2⁢(xn))n≥0 converge to x¯;

  • jj.

    for any n=0,1,…, one has

    |x¯−xn+1|≤max⁡{|xn+1−φ1⁢(xn)|,|xn+1−φ2⁢(xn)|};
  • jjj.

    there exists k∈ℝ, k>0,which does not depend on n∈ℕ such that

    |xn+1−x¯|≤k⁢|xn−x¯|3,n=0,1,…
Proof.

By φ1⁢(x0),φ2⁢(x0)∈[α,β] it obviously follows that h′⁢(φ1⁢(x0))>0 and h′⁢(φ2⁢(x0))>0. Denote by I0 the interval having the extremities φ1⁢(x0) and φ2⁢(x0). We notice, taking into account the mean formula, that

[φ1⁢(x0),φ2⁢(x0);h]>0.

When x0<x¯, by iv. and x¯=φ1⁢(x¯) it follows φ1⁢(x0)<x¯. Analogously, φ1⁢(x0)>x¯ for x0>x¯. Taking into account v. and x¯=φ2⁢(x¯) we get φ2⁢(x0)>x¯ for x0<x¯ and φ2⁢(x0)<x¯ for x0>x¯. It is obvious that in both situations x¯∈I. It can be easily seen that for all n=0,1,… we have

φ1⁢(xn)−h⁢(φ1⁢(xn))[φ1⁢(xn),φ2⁢(xn);h]=φ2⁢(xn)−h⁢(φ2⁢(xn))[φ1⁢(xn),φ2⁢(xn);h],

which, for n=0 imply x1>φ1⁢(x0) and x1<φ2⁢(x0) if x0<x¯, respectively x1<φ1⁢(x0) and x1>φ2⁢(x0) if x0>x¯, i.e., x1∈intI0. It is clear now that, analogously, x1<φ1⁢(x1)<x¯<φ2⁢(x1) if x1<x¯ or x1>φ1⁢(x1)>x¯>φ2⁢(x1) if x1>x¯. Denoting by I1 the interval determined by φ1⁢(x1) and φ2⁢(x1) then

I1⊂I0,

and the element x2 constructed by (9) satisfies x2,x¯∈I1.

Repeating the above reason we get xn+1,x¯∈In, the interval being determined by φ1⁢(xn),φ2⁢(xn), and also that

In+1⊂In.

and x¯∈In+1. It is clear now that jj. holds. In order to obtain jjj. we shall use the identity

h⁢(x¯)= h⁢(φ1⁢(xn))+[φ1⁢(xn),φ2⁢(xn);h]⁢(x¯−φ1⁢(xn))
+[x¯,φ1⁢(xn),φ2⁢(xn);h]⁢(x¯−φ1⁢(xn))⁢(x¯−φ2⁢(xn)),

which, together with (9) and h⁢(x¯)=0, imply

x¯−xn+1=−[x¯,φ1⁢(xn),φ2⁢(xn);h][φ1⁢(xn),φ2⁢(xn);h]⁢(x¯−φ1⁢(xn))⁢(x¯−φ2⁢(xn)).

For the difference x¯−φ1⁢(xn), by iv. one gets

x¯−φ1⁢(xn)=[x¯,xn;φ1]⁢(x¯−xn),

i.e.,

x¯−φ1⁢(xn)<|x¯−xn|.

Analogously, by v. we get

|x¯−φ2⁢(xn)|<|x¯−xn|.

The mean formula for divided differences implies

[x¯,φ1⁢(xn),φ2⁢(xn);h] =12⁢h′′⁢(ξn),with ⁢ξn∈In,and
[φ1⁢(xn),φ2⁢(xn);h] =h′⁢(ηn),ηn∈In.

For h′′⁢(ξn) we have

|h′′⁢(ξn)|=|h′′⁢(ξn)−h′′⁢(x¯)|=|h′′′⁢(θn)|⁢|x¯−ξn|.

Since ξn∈In it follows

|x¯−ξn|<|x¯−xn|.

Denoting m1=infx∈[α,β]|h′⁢(x)|, M3=supx∈[α,β]|h′′′⁢(x)|, the above relations lead to

|x¯−xn+1|≤M32⁢m1⁢|x¯−xn|3,n=0,1,…,

i.e., jjj. for k=M32⁢m1.

Since the initial approximation x0 was supposed sufficiently close to the solution x¯, then

M32⁢m1⁢|x¯−x0|<1

implies, together with jjj., statement j. ∎

Remark 2.1.

Supposing that f′′⁢(x)≥0 for all x∈[a,b], and if instead of iii. we assume that f′⁢(x)>0, then obviously f⁢(x)<0 for a≤x<x¯, and so in Theorem 2.1 we may take α=a. From the above conditions it follows that h′⁢(x)>0 for x∈[a,x¯], and since f⁢(x¯)=0, one gets β>x¯. ∎

3. DETERMINING THE FUNCTIONS φ1 AND φ2

Under reasonable hypotheses on f, we shall show that there exist two classes of functions among we can choose the functions φ1 and φ2 such that hypotheses iv. and v. to be satisfied.

Besides assumptions i.–iii. on f, we shall suppose that f is strictly convex, i.e., f′′⁢(x)>0, ∀x∈[a,b]. This condition implies that f′ is increasing on [a,b]. If, moreover, f′⁢(x)<2⁢λ, with 0<λ≤fr′⁢(a), then we may consider the functions

(10) φ1⁢(x) =x−f⁢(x)μand
(11) φ2⁢(x) =x−f⁢(x)λ

where μ may be taken as any real number greater than the left derivative of f at b, fl′⁢(b).

In the following we shall show that the functions φ1⁢(x) and φ2⁢(x) chosen above obey iv. and v. The derivatives of these functions are given by

φ1′⁢(x)=1−f′⁢(x)μ, resp. ⁢φ1′⁢(x)=1−f′⁢(x)λ.

Obviously, 0≤φ1′⁢(x)<1. Also, the monotonicity of f′ implies φ2′⁢(x)<0, while f′⁢(x)<2⁢λ implies −1<φ2′⁢(x).

Taking into account Remark 2.1, under the above hypotheses it is obvious that if x0<x¯, then condition φ1⁢(x0)∈[α,β] from Theorem 2.1 is obviously satisfied. Indeed, this fact follows from φ1′⁢(x)<1, since φ1⁢(x0)−x¯=φ1⁢(x0)−φ1⁢(x¯)=φ1′⁢(ξ)⁢(x0−x¯)<0,x0<θ<x¯, i.e., φ1⁢(x0)<x¯.

On the other hand, |φ1⁢(x0)−x¯|<|x0−x¯|, and so the relations x0<φ2⁢(x0)<x¯ hold. The hypothesis φ2⁢(x0)∈[α,β] must be kept.

References

  • [1] A. Ben-Israel, Newton’s method with modified functions, Contemp. Math., 204, pp. 39–50, 1997.
  • [2] G.H. Brown, Jr., On Halley’s variation of Newton’s method, Amer. Math. Monthly, 84, pp. 726–728, 1977.
  • [3] V. Candela and A. Marquina, Recurrence relations for rational cubic methods I: The Halley’s method, Computing, 44, pp. 169–184, 1990.
  • [4] G. Deslauries and S. Dubuc, Le calcul de la racine cubique selon Héron, El. Math., 51, pp. 28–34, 1996.
  • [5] W.F. Ford and J.A. Pennline, Accelerated convergence in Newton method, SIAM Rev., 38, pp. 658–659, 1996.
  • [6] J. Gerlach, Accelerated convergence in Newton’s method, SIAM Rev., 36, pp. 272–276, 1994.
  • [7] ††margin: clickable → D. Luca and I. Păvăloiu, On the Heron’s method for the approximation of the cubic root of a real number, Rev. Anal. Numér. Théor. Approx., 28, pp. 103–108, 1997.
  • [8] A. Melman, Geometry and convergence of Euler’s and Halley’s methods, SIAM Rev., 39, pp. 728–735, 1997.
  • [9] A.M. Ostrowski, The Solution of Equations and Systems of Equations, Academic Press, New York–London, 1960.
  • [10] ††margin: clickable → I. Păvăloiu, On the monotonicity of the sequences of approximations obtained by Stef- fensen method, Mathematica (Cluj), 35 (58), pp. 171–76, 1993.
  • [11] ††margin: clickable → I. Păvăloiu, Approximation of the roots of equations by Aitken–Steffensen-type mono- tonic sequences, Calcolo, 32, pp. 69–82, 1995.
  • [12] ††margin: clickable → I. Păvăloiu, On a Halley-Steffensen method for approximating the solutions of scalar equations, Rev. Anal. Numér. Théor. Approx., 30, 2001 no. 1, pp. 69–74.
  • [13] T. Popoviciu, ††margin: available soon,
    refresh and click here →
    Sur la délimitation de l’erreur dans l’approximation des racines d’une équation par interpolation linéaire ou quadratique, Rev. Roumaine Math. Pures Appl., 13, pp. 75–78, 1968.
  • [14]

Received June 13, 2000.

2001

Related Posts