A Zero-Dissipative Phase-Fitted Fourth Order Diagonally Implicit Runge-Kutta-Nyström Method for Solving Oscillatory Problems

<table>Graph for <svg height="21.8375" id="M138" style="vertical-align:-3.29552pt" version="1.1" viewbox="0 0 100.2 21.8375" width="100.2" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.675)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,10.041,16.675)"><path d="M238 0h-223v26q48 5 60 17t12 61v473q0 50 -10.5 60.5t-58.5 15.5v24q73 8 148 35v-608q0 -49 11.5 -61t60.5 -17v-26z" id="x6C"></path></g><g transform="matrix(.017,-0,0,-.017,14.341,16.675)"><path d="M257 449q92 0 154 -65t62 -158q0 -112 -67 -175t-150 -63q-98 0 -158.5 66.5t-60.5 154.5q0 59 21 106.5t54.5 75.5t71 43t73.5 15zM244 416q-48 0 -81 -47t-33 -128q0 -96 38 -158t99 -62q51 0 82 43.5t31 139.5q0 91 -36 151.5t-100 60.5z" id="x6F"></path></g><g transform="matrix(.017,-0,0,-.017,23.011,16.675)"><path d="M463 437l5 -8q-21 -38 -45 -51l-57 3q30 -35 30 -88q0 -76 -52.5 -118t-122.5 -42q-24 0 -51 5q-30 -21 -30 -42q0 -16 16 -27t47 -11q19 0 54 0.5t52 0.5q36 0 65.5 -8t53.5 -36t24 -74q0 -77 -75 -137.5t-174 -60.5q-79 0 -126 37.5t-48 84.5q0 32 24 56q21 21 83 70
q-49 15 -74 55q-13 20 -11 42q56 21 88 61q-37 16 -62.5 51t-25.5 83q0 78 56.5 122t122.5 44h1q51 0 96 -27q71 4 136 15zM219 418q-34 0 -58 -32t-24 -87q0 -59 26.5 -96.5t67.5 -37.5q35 1 59 32t24 87q0 59 -26 96.5t-69 37.5zM241 -11q-62 0 -87 -21q-46 -39 -46 -78
q0 -41 38 -68.5t94 -27.5q69 0 106.5 31t37.5 80q0 52 -45 71q-33 13 -98 13z" id="x67"></path></g>
<g transform="matrix(.012,-0,0,-.012,30.963,21.638)"><path 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" id="x31"></path></g><g transform="matrix(.012,-0,0,-.012,36.674,21.638)"><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,43.012,16.675)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g><g transform="matrix(.017,-0,0,-.017,49.999,16.675)"><path d="M471 456q-22 -60 -37 -85q-20 -3 -84 -3q-9 0 -29 2t-29 2q-38 -145 -54 -234q-15 -80 14 -80q23 0 68 38l14 -25q-78 -83 -137 -83q-50 0 -50 64q0 25 7 57l67 266q-57 4 -98 -9t-80 -49l-20 21q14 18 27 31t36.5 30t56.5 26t73 9q26 0 82.5 -2t89.5 -2q21 0 32 5t24 24
z" id="x1D70F"></path></g>
<g transform="matrix(.012,-0,0,-.012,58.212,8.513)"><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(.012,-0,0,-.012,62.329,8.513)"><path d="M570 304q0 -108 -87 -199q-40 -42 -94.5 -74t-105.5 -43q-41 0 -65 11l-29 -141q-9 -45 -1.5 -58t45.5 -16l26 -2l-5 -29l-241 -10l4 26q51 10 67.5 24t26.5 60l113 520q-54 -20 -89 -41l-7 26q38 28 105 53l11 49q20 25 77 58l8 -7l-17 -77q39 14 102 14q82 0 119 -36
t37 -108zM482 289q0 114 -113 114q-26 0 -66 -7l-70 -327q12 -14 32 -25t39 -11q59 0 118.5 81.5t59.5 174.5z" id="x1D45D"></path></g><g transform="matrix(.012,-0,0,-.012,69.385,8.513)"><path d="M535 230h-212v-233h-58v233h-213v50h213v210h58v-210h212v-50z" id="x2B"></path></g><g transform="matrix(.012,-0,0,-.012,76.37,8.513)"><path 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" id="x31"></path></g><g transform="matrix(.012,-0,0,-.012,82.081,8.513)"><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>
<g transform="matrix(.017,-0,0,-.017,86.812,16.675)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(.012,-0,0,-.012,93.8,20.75)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg> and <svg height="22.8125" id="M139" style="vertical-align:-3.29552pt" version="1.1" viewbox="0 0 103.7375 22.8125" width="103.7375" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,17.65)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,10.041,17.65)"><path d="M238 0h-223v26q48 5 60 17t12 61v473q0 50 -10.5 60.5t-58.5 15.5v24q73 8 148 35v-608q0 -49 11.5 -61t60.5 -17v-26z" id="x6C"></path></g><g transform="matrix(.017,-0,0,-.017,14.341,17.65)"><path d="M257 449q92 0 154 -65t62 -158q0 -112 -67 -175t-150 -63q-98 0 -158.5 66.5t-60.5 154.5q0 59 21 106.5t54.5 75.5t71 43t73.5 15zM244 416q-48 0 -81 -47t-33 -128q0 -96 38 -158t99 -62q51 0 82 43.5t31 139.5q0 91 -36 151.5t-100 60.5z" id="x6F"></path></g><g transform="matrix(.017,-0,0,-.017,23.011,17.65)"><path d="M463 437l5 -8q-21 -38 -45 -51l-57 3q30 -35 30 -88q0 -76 -52.5 -118t-122.5 -42q-24 0 -51 5q-30 -21 -30 -42q0 -16 16 -27t47 -11q19 0 54 0.5t52 0.5q36 0 65.5 -8t53.5 -36t24 -74q0 -77 -75 -137.5t-174 -60.5q-79 0 -126 37.5t-48 84.5q0 32 24 56q21 21 83 70
q-49 15 -74 55q-13 20 -11 42q56 21 88 61q-37 16 -62.5 51t-25.5 83q0 78 56.5 122t122.5 44h1q51 0 96 -27q71 4 136 15zM219 418q-34 0 -58 -32t-24 -87q0 -59 26.5 -96.5t67.5 -37.5q35 1 59 32t24 87q0 59 -26 96.5t-69 37.5zM241 -11q-62 0 -87 -21q-46 -39 -46 -78
q0 -41 38 -68.5t94 -27.5q69 0 106.5 31t37.5 80q0 52 -45 71q-33 13 -98 13z" id="x67"></path></g>
<g transform="matrix(.012,-0,0,-.012,30.963,22.613)"><path 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" id="x31"></path></g><g transform="matrix(.012,-0,0,-.012,36.674,22.613)"><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,43.012,17.65)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g><g transform="matrix(.017,-0,0,-.017,49.999,17.65)"><path d="M471 456q-22 -60 -37 -85q-20 -3 -84 -3q-9 0 -29 2t-29 2q-38 -145 -54 -234q-15 -80 14 -80q23 0 68 38l14 -25q-78 -83 -137 -83q-50 0 -50 64q0 25 7 57l67 266q-57 4 -98 -9t-80 -49l-20 21q14 18 27 31t36.5 30t56.5 26t73 9q26 0 82.5 -2t89.5 -2q21 0 32 5t24 24
z" id="x1D70F"></path></g>
<g transform="matrix(.012,-0,0,-.012,58.212,9.488)"><path d="M227 744l-123 -338l-31 15l73 368q12 3 41.5 -8t36.5 -20z" id="x2032"></path></g><g transform="matrix(.012,-0,0,-.012,61.746,9.488)"><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(.012,-0,0,-.012,65.863,9.488)"><path d="M570 304q0 -108 -87 -199q-40 -42 -94.5 -74t-105.5 -43q-41 0 -65 11l-29 -141q-9 -45 -1.5 -58t45.5 -16l26 -2l-5 -29l-241 -10l4 26q51 10 67.5 24t26.5 60l113 520q-54 -20 -89 -41l-7 26q38 28 105 53l11 49q20 25 77 58l8 -7l-17 -77q39 14 102 14q82 0 119 -36
t37 -108zM482 289q0 114 -113 114q-26 0 -66 -7l-70 -327q12 -14 32 -25t39 -11q59 0 118.5 81.5t59.5 174.5z" id="x1D45D"></path></g><g transform="matrix(.012,-0,0,-.012,72.919,9.488)"><path d="M535 230h-212v-233h-58v233h-213v50h213v210h58v-210h212v-50z" id="x2B"></path></g><g transform="matrix(.012,-0,0,-.012,79.904,9.488)"><path 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" id="x31"></path></g><g transform="matrix(.012,-0,0,-.012,85.615,9.488)"><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>
<g transform="matrix(.017,-0,0,-.017,90.35,17.65)"><path d="M162 -163h-60v866h60v-866zM309 -163h-60v866h60v-866z" id="x2016"></path></g>
<g transform="matrix(.012,-0,0,-.012,97.338,21.725)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg> against <svg height="8.4250002" id="M140" style="vertical-align:-0.17555pt" version="1.1" viewbox="0 0 8.5874996 8.4250002" width="8.5874996" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,8.013)"><path d="M475 445l-64 -69q-119 -128 -289 -280q27 6 81 -15q66 -26 102 -26q48 0 111 93l23 -19q-41 -84 -75.5 -116.5t-73.5 -32.5q-37 0 -147 53q-43 21 -63 7q-22 -15 -38 -41q-4 -7 -11 -7q-8 10 -8 26q0 38 47 83l292 273q-32 -5 -67 -1q-63 7 -105 7q-32 0 -53.5 -22
t-47.5 -67l-25 13q15 39 38 78q40 66 83 66q57 0 120 -12q52 -10 80 -10q29 0 67 42z" id="x1D467"></path></g>
</svg> value.</table>

Mathematical Problems in Engineering

fig2

Figure 2

Figure 2: A Zero-Dissipative Phase-Fitted Fourth Order Diagonally Implicit Runge-Kutta-Nyström Method for Solving Oscillatory Problems 