Two Embedded Pairs of Runge-Kutta Type Methods for Direct Solution of Special Fourth-Order Ordinary Differential Equations

<div>The graph of <svg height="15.4397pt" id="M248" style="vertical-align:-3.051001pt" version="1.1" viewbox="-0.0498162 -12.3887 37.3021 15.4397" width="37.3021pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M167 -163V704H103V-163H167ZM317 -163V704H253V-163H317Z" id="g50-10"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="420" vert-adv-y="420"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,5.46,0)"><path d="M487 458L455 460C436 434 425 433 397 433C331 433 279 437 225 437C103 437 52 374 24 338L49 312C100 363 152 379 229 374C210 308 184 213 163 117C158 94 156 74 156 58C156 6 180 -12 211 -12C261 -12 313 27 348 68L332 99C307 76 285 64 267 64C252 64 236 86 250 148L301 373C320 373 344 369 364 369C399 369 435 370 450 372C458 386 470 411 487 458Z" id="g50-241"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="499" vert-adv-y="499"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,11.95,-5.872)"><path d="M310 541L304 571C290 586 211 619 185 610L80 76L131 52L310 541Z" id="g50-31"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="361" vert-adv-y="361"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,15.235,-5.872)"><path d="M309 -142C211 -61 151 99 151 271S210 601 309 683L288 710C141 611 75 451 75 272V271C75 90 141 -72 288 -169L309 -142Z" id="g50-41"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,18.493,-5.872)"><path d="M141 347C171 480 264 541 327 570C364 587 403 599 436 605L429 641C389 634 339 622 300 609C191 570 37 458 37 239C37 86 126 -12 245 -12C367 -12 454 89 454 210C454 314 380 397 272 397C252 397 230 390 208 380L141 347ZM231 338C323 338 366 257 366 174C366 108 340 27 263 27C180 27 129 123 129 244C129 268 130 292 135 310C158 323 193 338 231 338Z" id="g50-55"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,22.925,-5.872)"><path d="M283 271C283 451 220 610 70 710L48 683C146 603 207 443 207 271S149 -59 48 -142L70 -169C220 -67 283 90 283 271Z" id="g50-42"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,26.734,0)"><path d="M167 -163V704H103V-163H167ZM317 -163V704H253V-163H317Z" id="g50-10"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="420" vert-adv-y="420"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,32.194,3.001)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg> versus <svg height="16.4288pt" id="M249" style="vertical-align:-3.852201pt" version="1.1" viewbox="-0.0498162 -12.5766 10.1742 16.4288" width="10.1742pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,-.535,-3.51)"><path d="M490 563C409 598 330 631 258 673H233C161 631 81 597 0 563L14 533C89 551 168 576 245 598C320 575 400 551 476 533L490 563Z" id="g56-17"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="490" vert-adv-y="490"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M460 334C460 396 434 451 378 451C330 451 241 408 148 300H146L237 679C241 697 241 710 232 710C213 710 153 684 67 675L66 646H95C141 646 145 642 134 595L39 170C23 97 31 54 46 33C64 8 100 -12 137 -12C178 -12 234 3 298 43C391 101 460 222 460 334ZM371 320C371 204 316 89 248 51C230 41 208 37 192 37C143 37 102 72 119 166C124 194 129 215 135 235C202 323 298 392 335 392C353 392 371 372 371 320Z" id="g50-99"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="478" vert-adv-y="478"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.218,-5.872)"><path d="M310 541L304 571C290 586 211 619 185 610L80 76L131 52L310 541Z" id="g50-31"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="361" vert-adv-y="361"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,5.109,3.693)"><path d="M290 377C321 398 342 415 358 430C378 450 389 473 389 502C389 578 329 635 238 635H237C184 635 137 610 109 578L64 515L88 493C112 529 154 573 208 573S303 542 303 482C303 409 233 370 141 341L149 308C165 313 190 319 215 319C272 319 341 283 341 193C342 98 292 43 222 43C163 43 122 72 96 94C88 101 79 100 70 94C61 87 47 73 46 60C44 47 48 37 62 23C76 10 118 -12 165 -12C238 -12 430 62 430 223C430 297 379 359 290 375V377Z" id="g50-52"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>.</div>

Mathematical Problems in Engineering

fig2

Figure 2

Figure 2: Two Embedded Pairs of Runge-Kutta Type Methods for Direct Solution of Special Fourth-Order Ordinary Differential Equations 