Phase-Fitted and Amplification-Fitted Higher Order Two-Derivative Runge-Kutta Method for the Numerical Solution of Orbital and Related Periodical IVPs

<div>Energy conservation. The logarithm error of energy (MAXERR) at each integration point when solving the harmonic oscillator (Problem <a href="https://static.hindawi.com/articles/mpe/volume-2017/1871278/figures/#problem1">1</a>) for <svg height="8.55521pt" id="M167" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -8.34882 29.6114 8.55521" width="29.6114pt" 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="M615 290C615 392 563 448 529 448C511 448 491 441 484 434L483 425C518 400 546 349 546 273C546 158 504 46 417 46C381 46 356 77 356 142C356 184 380 293 392 334L387 340L307 323C305 288 296 228 287 188C265 94 221 52 180 52C144 52 112 89 112 171C112 278 164 372 241 430L226 454C106 394 23 288 23 155C23 40 79 -12 139 -12C184 -12 245 32 285 85C292 31 322 -12 376 -12C423 -12 470 18 501 45C570 104 615 184 615 290Z" id="g113-246"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="638" vert-adv-y="638"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.927,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,23.19,0)"><path d="M249 635C141 635 70 555 70 471C70 401 114 353 179 316C143 294 106 267 90 252C68 231 45 202 45 157C45 50 130 -12 237 -12C322 -12 435 52 435 169C435 256 372 304 303 343C349 374 375 398 383 407C401 429 411 458 411 487C411 569 344 635 249 635ZM238 603C285 603 337 567 337 482C337 422 310 385 276 358C205 393 145 426 145 500C145 552 179 603 238 603ZM248 20C183 20 125 70 125 163C125 218 158 268 206 300C284 261 355 217 355 143C355 66 308 20 248 20Z" id="g113-57"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg> with initial condition <svg height="13.6311pt" id="M168" style="vertical-align:-5.28228pt" version="1.1" viewbox="-0.0498162 -8.34882 36.2957 13.6311" width="36.2957pt" 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="M474 429L457 433C435 440 389 448 367 448C348 448 323 446 309 443C266 434 195 406 148 366C78 307 23 210 23 101C23 35 55 -12 92 -12C118 -12 146 1 196 35C247 70 311 130 346 173H348L281 -148C268 -211 256 -221 208 -229L191 -232L187 -257L433 -245L437 -219L411 -216C357 -210 350 -205 362 -140L427 207C447 315 461 381 474 429ZM387 387C379 337 363 262 355 236C318 180 201 57 142 57C126 57 112 81 112 128C112 205 150 321 220 376C244 395 280 403 312 403C345 403 370 396 387 387Z" id="g113-114"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="491" vert-adv-y="491"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.318,3.132)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0065,0,0,-0.0065,10.75,5.233)"><path d="M394 0V37C301 42 295 50 295 122V635C236 613 178 596 109 583V553L165 551C205 549 209 544 209 495V122C209 50 203 42 111 37V0H394Z" id="g176-50"></path><glyph.data ascent="3439" descent="-2858" horiz-adv-x="494" vert-adv-y="494"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,18.618,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.88,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, <svg height="13.6311pt" id="M169" style="vertical-align:-5.28228pt" version="1.1" viewbox="-0.0498162 -8.34882 43.9553 13.6311" width="43.9553pt" 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="M474 429L457 433C435 440 389 448 367 448C348 448 323 446 309 443C266 434 195 406 148 366C78 307 23 210 23 101C23 35 55 -12 92 -12C118 -12 146 1 196 35C247 70 311 130 346 173H348L281 -148C268 -211 256 -221 208 -229L191 -232L187 -257L433 -245L437 -219L411 -216C357 -210 350 -205 362 -140L427 207C447 315 461 381 474 429ZM387 387C379 337 363 262 355 236C318 180 201 57 142 57C126 57 112 81 112 128C112 205 150 321 220 376C244 395 280 403 312 403C345 403 370 396 387 387Z" id="g113-114"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="491" vert-adv-y="491"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.318,3.132)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0065,0,0,-0.0065,10.75,5.233)"><path d="M452 139L416 148C386 81 373 77 319 77H141L282 220C372 312 414 373 414 464C414 569 332 635 240 635C181 635 132 609 101 574L42 493L70 467C98 506 143 562 209 562C280 562 324 515 324 435C324 347 261 272 197 197C148 139 89 84 34 29V0H417L452 139Z" id="g176-51"></path><glyph.data ascent="3439" descent="-2858" horiz-adv-x="494" vert-adv-y="494"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,18.618,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.88,0)"><path d="M535 230V280H52V230H535Z" id="g117-33"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,37.54,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g113-51"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, and <svg height="11.439pt" id="M170" style="vertical-align:-2.15067pt" version="1.1" viewbox="-0.0498162 -9.28833 45.7397 11.439" width="45.7397pt" 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="M600 0V24L339 665L311 657L43 24V0H600ZM497 50H107L303 539L497 50Z" id="g113-133"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="641" vert-adv-y="641"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,8.327,0)"><path d="M324 430H196L233 583L223 592L145 529L120 430H54L29 396L31 388H111L56 126C33 15 54 -12 77 -12C137 -12 214 57 250 95L233 119C208 92 155 59 138 59C126 59 120 70 131 125L186 390L298 394L324 430Z" id="g113-117"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="341" vert-adv-y="341"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,16.392,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,27.655,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.895,0)"><path d="M368 703H309L44 -163H104L368 703Z" id="g113-48"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="412" vert-adv-y="412"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,39.256,0)"><path d="M153 550H386L412 615L406 623H120L82 318C104 327 142 338 184 338C294 338 347 275 347 187C347 112 305 39 221 39C160 39 119 71 97 89C88 97 80 96 71 90C59 80 50 67 49 57C48 45 52 36 66 23C80 9 123 -12 169 -12C221 -11 288 15 342 59C403 109 431 165 431 225C431 308 366 395 238 395C212 395 165 379 127 364L153 550Z" id="g113-54"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg> with <svg height="11.6445pt" id="M171" style="vertical-align:-3.29568pt" version="1.1" viewbox="-0.0498162 -8.34882 58.8977 11.6445" width="58.8977pt" 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="M324 430H196L233 583L223 592L145 529L120 430H54L29 396L31 388H111L56 126C33 15 54 -12 77 -12C137 -12 214 57 250 95L233 119C208 92 155 59 138 59C126 59 120 70 131 125L186 390L298 394L324 430Z" id="g113-117"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="341" vert-adv-y="341"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,4.368,3.132)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,8.317,3.132)"><path d="M532 0V30C473 38 467 38 467 107V298C467 395 416 451 336 451C307 451 279 439 251 420C226 404 205 389 185 374V453C142 434 94 422 43 413V386C98 378 104 373 104 309V107C104 41 100 38 31 30V0H251V30C191 35 185 40 185 107V337C216 364 255 389 295 389C360 389 386 347 386 276V110C386 44 380 35 321 30V0H532Z" id="g58-108"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="557" vert-adv-y="557"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,13.359,3.132)"><path d="M526 56L493 58C457 61 448 68 448 119V710C411 697 344 682 288 675V646C362 641 367 639 367 573V439C344 448 315 451 300 451C164 451 40 342 40 202C40 61 147 -12 227 -12C239 -12 266 -7 305 16L367 53V-12C427 10 505 21 526 26V56ZM367 88C342 69 306 54 271 54C207 54 133 111 133 229C133 373 219 409 264 409C301 409 342 394 367 359V88Z" id="g58-98"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="538" vert-adv-y="538"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.437,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.7,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,39.94,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,46.18,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,52.42,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>.</div>

Mathematical Problems in Engineering

fig2

Figure 2

Figure 2: Phase-Fitted and Amplification-Fitted Higher Order Two-Derivative Runge-Kutta Method for the Numerical Solution of Orbital and Related Periodical IVPs 

Figure 2 | Phase-Fitted and Amplification-Fitted Higher Order Two-Derivative Runge-Kutta Method for the Numerical Solution of Orbital and Related Periodical IVPs 