About the Use of the HdHr Algorithm Group in Integrating the Movement Equation with Nonlinear Terms

<table>Logarithm for the relative error within the range <svg height="11.175" id="M174" style="vertical-align:-0.0pt" version="1.1" viewbox="0 0 13.5875 11.175" width="13.5875" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.113)"><path d="M775 650l-6 -28q-60 -6 -81.5 -16t-61.5 -54l-175 -191l125 -243q30 -58 48.5 -71t82.5 -19l-5 -28h-275l7 28l35 4q31 4 37 12t-6 34l-108 216q-140 -165 -177 -219q-16 -22 -10.5 -30.5t41.5 -13.5l22 -3l-7 -28h-244l8 28q52 4 75 15.5t67 52.5q48 46 206 231
l-110 215q-26 51 -44 63t-72 17l6 28h250l-6 -28l-27 -4q-30 -5 -35 -10t3 -27q17 -43 95 -190q70 78 154 185q15 21 10 29.5t-33 12.5l-30 4l5 28h236z" id="x1D44B"></path></g>
</svg>  step estimate <svg height="16.450001" id="M175" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 79.474998 16.450001" width="79.474998" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path 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" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,5.944,12.162)"><path d="M475 507q0 -83 -20 -172t-56 -167.5t-93.5 -129t-125.5 -50.5q-157 0 -157 227q0 78 21.5 164t59 161t96.5 123.5t126 48.5q79 0 114 -58t35 -147zM391 522q0 155 -81 155q-62 0 -111 -82.5t-73 -200.5h253q12 81 12 128zM373 346h-255q-12 -91 -12 -150q0 -72 20 -123
t63 -51q34 0 64 28.5t52.5 77t39 103t28.5 115.5z" id="x1D703"></path></g>
<g transform="matrix(.012,-0,0,-.012,14.413,16.25)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
<g transform="matrix(.017,-0,0,-.017,25.462,12.162)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,40.149,12.162)"><path 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" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,48.309,12.162)"><path d="M574 449q-23 -58 -38 -79q-16 -4 -79 -4q-27 -100 -46 -219q-14 -87 7 -87t74 36l13 -27q-74 -81 -139 -81q-48 0 -48 62q0 22 8 59l60 258l-154 4q-21 -103 -54.5 -209t-64.5 -151q-38 -23 -83 -23l-7 15q46 36 94 152.5t64 216.5q-81 0 -138 -54l-18 23q22 27 39.5 43
t61.5 33.5t100 17.5q43 0 131.5 -2.5t129.5 -2.5q23 0 33.5 5t24.5 25z" id="x1D70B"></path></g><g transform="matrix(.017,-0,0,-.017,58.355,12.162)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,65.359,12.162)"><path 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" id="x36"></path></g><g transform="matrix(.017,-0,0,-.017,73.518,12.162)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
</svg>.</table>

Mathematical Problems in Engineering

fig3

Figure 3

Figure 3: About the Use of the HdHr Algorithm Group in Integrating the Movement Equation with Nonlinear Terms 