The problem considered in this paper is to approximate the simple zeros of the function <svg style="vertical-align:-3.56265pt;width:32.299999px;" id="M1" height="16.6625" version="1.1" viewBox="0 0 32.299999 16.6625" width="32.299999" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns="http://www.w3.org/2000/svg"> <g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path id="x1D453" d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" /></g><g transform="matrix(.017,-0,0,-.017,10.976,12.162)"><path id="x28" d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" /></g><g transform="matrix(.017,-0,0,-.017,16.857,12.162)"><path id="x1D465" d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" /></g><g transform="matrix(.017,-0,0,-.017,26.36,12.162)"><path id="x29" d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" /></g> </svg> by iterative processes. An optimal 16th order class is constructed. The class is built by considering any of the optimal three-step derivative-involved methods in the first three steps of a four-step cycle in which the first derivative of the function at the fourth step is estimated by a combination of already known values. Per iteration, each method of the class reaches the efficiency index <svg style="vertical-align:-1.13977pt;width:85.974998px;" id="M2" height="17.125" version="1.1" viewBox="0 0 85.974998 17.125" width="85.974998" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns="http://www.w3.org/2000/svg"> <g transform="matrix(.008,-0,0,-.008,3.525,6.012)"><path id="x35" d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" /></g> <g transform="matrix(.017,-0,0,-.017,.663,12.862)"><path id="x221A" d="M782 703l-309 -950l-231 505l-91 -47l-18 35l147 78l194 -431h2l264 810h42z" /></g> <g transform="matrix(1.25,0,0,-1.25,90,90)"> <path d="m -61.08,70.93 13.06,0" style="fill:none;stroke:#000000;stroke-width:0.68000001;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:10;stroke-opacity:1;stroke-dasharray:none" /> </g> <g transform="matrix(.017,-0,0,-.017,13.65,15.637)"><path id="x31" d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" /></g><g transform="matrix(.017,-0,0,-.017,21.809,15.637)"><path id="x36" d="M137 343l67 33q37 17 63 17q79 0 129.5 -53t50.5 -131q0 -92 -58 -156.5t-147 -64.5t-147 68t-58 182q0 63 17 119t43 95.5t61.5 72t69 52t67.5 31.5q62 22 128 33l6 -32q-56 -11 -108 -35q-149 -71 -184 -231zM227 337q-47 0 -95 -27q-6 -23 -6 -70q0 -93 36 -155.5
t96 -62.5q53 0 78 45.5t25 105.5q0 68 -35 116t-99 48z" /></g><g transform="matrix(.017,-0,0,-.017,34.694,15.637)"><path id="x2248" d="M537 353q-54 -70 -120 -70q-47 0 -137 43q-78 38 -108 38q-55 0 -91 -57l-31 33q51 74 122 74q45 0 137 -44q76 -37 108 -37q55 0 91 57zM537 168q-54 -70 -120 -70q-47 0 -137 43q-78 38 -108 38q-55 0 -91 -57l-31 33q51 74 122 74q45 0 137 -44q76 -37 108 -37
q55 0 91 57z" /></g><g transform="matrix(.017,-0,0,-.017,49.381,15.637)"><use xlink:href="#x31"/></g><g transform="matrix(.017,-0,0,-.017,57.541,15.637)"><path id="x2E" d="M113 -12q-24 0 -39.5 16t-15.5 42q0 24 16 40.5t40 16.5t40 -16.5t16 -40.5q0 -26 -16 -42t-41 -16z" /></g><g transform="matrix(.017,-0,0,-.017,61.416,15.637)"><path id="x37" d="M447 623l8 -12l-283 -613l-74 -10l-7 11q174 283 297 551h-216q-48 0 -62.5 -12t-33.5 -63h-29q10 60 18 148h382z" /></g><g transform="matrix(.017,-0,0,-.017,69.576,15.637)"><path id="x34" d="M456 178h-96v-72q0 -51 12.5 -62.5t72.5 -16.5v-27h-256v27q65 5 78 17t13 62v72h-260v28q182 271 300 426h40v-407h96v-47zM280 225v295h-2q-107 -148 -196 -295h198z" /></g><g transform="matrix(.017,-0,0,-.017,77.735,15.637)"><use xlink:href="#x31"/></g> </svg>, by carrying out four evaluations of the function and one evaluation of the first derivative. The error equation for one technique of the class is furnished analytically. Some methods of the class are tested by challenging the existing high-order methods. The interval Newton's method is given as a tool for extracting enough accurate initial approximations to start such high-order methods. The obtained numerical results show that the derived methods are accurate and efficient.