The Numerical Method for Solving Differential Equations of Lane-Emden Type by Padé Approximation

<table>Values of <svg height="16.625" id="M97" style="vertical-align:-3.56265pt" version="1.1" viewbox="0 0 31.237499 16.625" width="31.237499" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.138)"><path d="M556 393q0 -39 -36 -106q-42 -78 -185 -279q-47 -66 -81 -108t-117 -135l-112 -26l-8 22q150 90 251 219q-6 136 -39 340q-8 53 -21 53q-6 0 -27 -19.5t-38 -42.5l-16 26q80 111 127 111q23 0 35 -28t20 -90q18 -137 27 -263h2q142 200 142 279q0 24 -14 48q-4 7 5 26
q13 28 43 28q18 0 30 -15.5t12 -39.5z" id="x1D466"></path></g><g transform="matrix(.017,-0,0,-.017,9.905,12.138)"><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,15.786,12.138)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,25.289,12.138)"><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> and its <svg height="14.85" id="M98" style="vertical-align:-2.22495pt" version="1.1" viewbox="0 0 45.987499 14.85" width="45.987499" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.012)"><path d="M619 482q0 -69 -40.5 -119t-96 -72.5t-118.5 -26.5h-44l-70 20l-31 -151q-14 -67 0.5 -83t89.5 -22l-5 -28h-287l8 28q65 7 81.5 22t29.5 83l79 398q12 56 0.5 70.5t-78.5 20.5l7 28h255q108 0 164 -43t56 -125zM524 478q0 141 -146 141q-25 0 -47 -8q-16 -6 -20.5 -13.5
t-10.5 -39.5l-43 -241q37 -13 83 -13q67 0 125.5 45t58.5 129z" id="x1D443"></path></g><g transform="matrix(.017,-0,0,-.017,10.874,12.012)"><path d="M290 -163h-170v866h170v-28q-79 -7 -94 -19.5t-15 -72.5v-627q0 -59 14.5 -71.5t94.5 -19.5v-28z" id="x5B"></path></g><g transform="matrix(.017,-0,0,-.017,16.738,12.012)"><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,24.898,12.012)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,31.901,12.012)"><path d="M456 178h-96v-72q0 -51 12.5 -62.5t72.5 -16.5v-27h-256v27q65 5 78 17t13 62v72h-260v28q182 271 300 426h40v-407h96v-47zM280 225v295h-2q-107 -148 -196 -295h198z" id="x34"></path></g><g transform="matrix(.017,-0,0,-.017,40.061,12.012)"><path d="M226 -163h-170v27q79 7 94 20t15 73v627q0 59 -15 72t-94 20v27h170v-866z" id="x5D"></path></g>
</svg> Padé approximant.</table>

Discrete Dynamics in Nature and Society

fig2

Figure 2

Figure 2: The Numerical Method for Solving Differential Equations of Lane-Emden Type by Padé Approximation 