Stability of Runge-Kutta Methods for Neutral Delay Differential Equations

<div>The norm of difference <svg height="12.6917pt" id="M184" style="vertical-align:-3.40337pt" version="1.1" viewbox="-0.0498162 -9.28833 48.5993 12.6917" width="48.5993pt" 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="M563 395C563 428 543 451 519 451C500 451 480 438 472 423C468 415 465 403 470 394C477 381 482 368 482 347C482 272 386 136 342 74H340C334 163 323 257 312 337C302 410 290 451 257 451C214 451 160 385 127 338L145 310C169 341 203 373 209 373C219 373 225 369 232 325C249 215 269 62 271 -15C220 -80 135 -165 18 -231L26 -257L142 -230C254 -105 279 -74 342 13C394 85 484 211 526 286C553 334 563 365 563 395Z" id="g50-122"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.506,3.001)"><path d="M504 89L488 119C453 85 428 69 416 69C408 69 405 75 412 102L462 304C492 438 460 451 436 451S388 441 356 423C311 397 229 336 170 256H168L189 340C208 418 200 451 177 451C144 451 82 413 24 352L40 322C65 345 98 369 108 369C114 369 119 363 112 335L28 -3L36 -12C54 -3 83 4 112 10C125 73 138 122 153 174C205 258 328 382 376 382C395 382 399 369 384 305L330 93C312 25 323 -12 351 -12C378 -12 439 21 504 89Z" id="g50-111"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="512" vert-adv-y="512"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,20.621,0)"><path d="M556 236V289H56V236H556Z" id="g54-33"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.504,0)"><path d="M481 445L456 468C426 437 403 427 388 427C372 427 356 428 302 438C256 447 221 451 187 451C159 451 130 428 105 384C93 363 80 340 65 305L96 290C122 342 149 381 193 381C216 381 267 377 294 374S346 372 365 376L75 105C37 70 24 35 24 16C24 2 30 -5 35 -8C42 -8 46 -6 48 -2C58 14 75 35 90 41C102 46 116 46 147 32C200 8 267 -21 298 -21C347 -21 390 12 448 126L420 151C369 77 338 56 311 56C287 56 249 68 212 81C180 92 149 100 127 98C212 172 353 308 413 372L481 445Z" id="g50-123"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="505" vert-adv-y="505"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,37.744,3.001)"><path d="M504 89L488 119C453 85 428 69 416 69C408 69 405 75 412 102L462 304C492 438 460 451 436 451S388 441 356 423C311 397 229 336 170 256H168L189 340C208 418 200 451 177 451C144 451 82 413 24 352L40 322C65 345 98 369 108 369C114 369 119 363 112 335L28 -3L36 -12C54 -3 83 4 112 10C125 73 138 122 153 174C205 258 328 382 376 382C395 382 399 369 384 305L330 93C312 25 323 -12 351 -12C378 -12 439 21 504 89Z" id="g50-111"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="512" vert-adv-y="512"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,42.919,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></svg> with step size <svg height="9.49473pt" id="M185" style="vertical-align:-0.2063999pt" version="1.1" viewbox="-0.0498162 -9.28833 57.1661 9.49473" width="57.1661pt" 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="M509 90L493 118C465 91 427 70 418 70S409 77 415 103L469 327C494 429 466 451 436 451C412 451 387 444 356 427C307 400 222 342 160 264H158L246 667C253 701 252 710 243 710C228 710 169 680 87 673L83 644H117C155 644 160 640 149 588L24 -4L32 -12C57 -4 86 4 108 9L143 187C198 275 324 386 375 386C394 386 399 367 385 311L331 87C315 22 322 -12 351 -12C387 -12 452 24 509 90Z" id="g50-105"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="527" vert-adv-y="527"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,10.527,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.145,0)"><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(.013,0,0,-0.013,28.476,0)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.557,0)"><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(.013,0,0,-0.013,37.888,0)"><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(.013,0,0,-0.013,44.219,0)"><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><g transform="matrix(.013,0,0,-0.013,50.55,0)"><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>

Discrete Dynamics in Nature and Society

fig2

Figure 2

Figure 2: Stability of Runge-Kutta Methods for Neutral Delay Differential Equations 