Metode Frobenius

Batas sumber. Bagian ini mengikat special_functions.tex baris 1-673 pada sumber beku.

Tinjauan Beberapa Fungsi Khusus

Dalam bab ini, kami menyajikan secara sangat singkat beberapa fungsi khusus yang digunakan untuk menyatakan solusi umum persamaan diferensial parsial linear orde dua, seperti persamaan panas, persamaan gelombang, dan persamaan potensial. Rujukan yang baik untuk materi bab ini adalah Simm; RenRog.

“Polinom Legendre” yang akan kita perkenalkan kemudian didefinisikan sebagai solusi persamaan diferensial tertentu yang diperoleh menggunakan “metode Frobenius” di sekitar “titik reguler.” Demikian pula, fungsi Bessel yang juga akan kita perkenalkan kemudian didefinisikan sebagai solusi persamaan diferensial tertentu yang diperoleh menggunakan “metode Frobenius” di sekitar “titik singular reguler.” Kita mulai dengan menyajikan metode Frobenius.

Metode Frobenius

Dalam kondisi apa solusi yy dari suatu persamaan diferensial biasa bersifat analitik di sekitar titik aa? Dengan kata lain, dalam kondisi apa kita memperoleh y=n=0an(xa)n\displaystyle y = \sum_{n=0}^\infty a_n (x-a)^n, yakni deret dengan jari-jari kekonvergenan positif, sebagai solusi suatu persamaan diferensial biasa? Contoh berikut menunjukkan bahwa kita tidak dapat selalu mengasumsikan bahwa hal ini mungkin.

Contoh. Tinjau persamaan diferensial biasa

x2y+(x1)y+1=0.x^2 y' + (x-1)y + 1 = 0 \ . Misalkan kita mempunyai solusi berbentuk y=n=0anxn\displaystyle y = \sum_{n=0}^\infty a_n x^n dengan jari-jari kekonvergenan positif. Maka y=n=1nanxn1\displaystyle y' = \sum_{n=1}^\infty n a_n x^{n-1} pada interval konvergensi. Jika deret-deret ini kita substitusikan sebagai yy dan yy' ke dalam (5.1.1), maka kita memperoleh

0=n=1nanxn+1+(x1)n=0anxn+1=n=1nanxn+1+n=0anxn+1n=0anxn+1=n=2(n1)an1xn+n=1an1xnn=0anxn+1=n=2(nan1an)xn+(a0a1)xa0+1.\begin{aligned} 0 &= \sum_{n=1}^\infty n a_n x^{n+1} +(x-1) \sum_{n=0}^\infty a_n x^n + 1 = \sum_{n=1}^\infty n a_n x^{n+1} + \sum_{n=0}^\infty a_n x^{n+1} - \sum_{n=0}^\infty a_n x^n + 1 \\ &= \sum_{n=2}^\infty (n-1) a_{n-1} x^n + \sum_{n=1}^\infty a_{n-1} x^n - \sum_{n=0}^\infty a_n x^n + 1 \\ &= \sum_{n=2}^\infty \left( n a_{n-1} - a_n\right) x^n + (a_0 - a_1)x - a_0 + 1 \ . \end{aligned} Dengan demikian, kita memperoleh a0=1a_0 = 1 dari 1a0=01-a_0 = 0, a1=a0=1a_1 = a_0 = 1 dari a0a1=0a_0 - a_1 = 0, a2=2a1=2a_2 = 2a_1 = 2 dari 2a1a2=02 a_1 - a_2 = 0, a3=3a2=32=3!a_3 = 3a_2 = 3\cdot 2 = 3! dari 3a2a3=03 a_2 - a_3 = 0, dan seterusnya. Jika kita mengasumsikan bahwa an1=(n1)!a_{n-1} = (n-1)!, maka kita memperoleh an=n!a_n = n! dari nan1an=0n a_{n-1} - a_n = 0. Oleh induksi, kita mempunyai an=n!a_n = n! untuk setiap bilangan bulat tak negatif. Kita memperoleh y=n=0n!xn\displaystyle y = \sum_{n=0}^\infty n! x^n. Uji rasio menunjukkan bahwa deret ini divergen untuk x0x\neq 0. Oleh karena itu, kita tidak dapat mengasumsikan bahwa terdapat solusi analitik dari (5.1.1) di lingkungan titik asal.

Tinjau persamaan diferensial biasa berbentuk

y+P(x)y+Q(x)y=0.y'' + P(x) y' + Q(x) y = 0 \ . Jika P(x)P(x) dan Q(x)Q(x) merupakan fungsi analitik pada x=0x=0, kita mengatakan bahwa 00 adalah suatu titik reguler. Jika xP(x)x P(x) dan x2Q(x)\displaystyle x^2 Q(x) merupakan fungsi analitik pada x=0x=0, kita mengatakan bahwa 00 adalah suatu titik singular reguler.

Andaikan bahwa 00 merupakan titik singular reguler; misalkan xP(x)=j=0αjxj\displaystyle x P(x) = \sum_{j=0}^\infty \alpha_j x^j dan x2Q(x)=j=0βjxj\displaystyle x^2 Q(x) = \sum_{j=0}^\infty \beta_j x^j untuk |x|<R|x| < R. Untuk cabang Frobenius real yang digunakan di bawah, kita andaikan bahwa kedua akar indisial real dan bekerja pada 0<x<R0<x<R. Untuk akar kompleks atau domain kompleks, cabang logaritma yang menentukan xλx^\lambda harus dipilih secara eksplisit. Selanjutnya, andaikan λ1\lambda_1 dan λ2\lambda_2 merupakan akar-akar dari persamaan indisial F(λ)=λ2+(α01)λ+β0=0F(\lambda) = \lambda^2 + (\alpha_0 - 1) \lambda + \beta_0 = 0 dengan λ1λ2\lambda_1 \geq \lambda_2.

Salah satu solusi persamaan diferensial biasa tersebut diberikan oleh

y1(x)=xλ1(a0+a1x+a2x2+),y_1(x) = x^{\lambda_1}(a_0 + a_1 x + a_2 x^2 + \ldots) \ , dengan deret yang konvergen secara mutlak untuk 0<x<R0<x<R pada cabang real tersebut. Solusi kedua yang bebas secara linear dari persamaan diferensial biasa tersebut diberikan oleh

y2(x)={xλ2(b0+b1x+b2x2+) jika λ1λ2y1(x)ln(x)+xλ2(b1x+b2x2+) jika λ1=λ2ky1(x)ln(x)+xλ2(b0+b1x+b2x2+) jika 0<λ1λ2y_2(x) = \begin{cases} x^{\lambda_2} ( b_0 + b_1 x + b_2 x^2 + \ldots ) & \text{ jika } \lambda_1 - \lambda_2 \not\in \mathbb{N}\\ y_1(x) \ln(x) + x^{\lambda_2} ( b_1 x + b_2 x^2 + \ldots ) & \text{ jika } \lambda_1 = \lambda_2 \\ k y_1(x) \ln(x) + x^{\lambda_2} ( b_0 + b_1 x + b_2 x^2 + \ldots ) & \text{ jika } 0<\lambda_1 - \lambda_2 \in \mathbb{N} \end{cases} dengan deret yang juga konvergen secara mutlak untuk 0<x<R0<x<R pada cabang real tersebut dan kk suatu konstanta.

Dalam beberapa keadaan, kita dapat menghindari manipulasi rumus untuk y2y_2 yang memuat suku ln(x)\ln(x). Jika y1y_1 merupakan fungsi yang diketahui, kita dapat menggunakan metode reduksi orde untuk memperoleh solusi kedua yang bebas secara linear; cukup dengan mengganti yy dengan y2(x)=y1(x)u(x)y_2(x) = y_1(x) u(x) dalam persamaan diferensial biasa tersebut lalu menyelesaikannya untuk uu.

Jika 0<s=λ1λ20<s = \lambda_1-\lambda_2 \in \mathbb{N}, solusi y2y_2 dengan k=0k=0 memberikan solusi kedua yang bebas secara linear ketika b0b_0 dan bsb_s bebas.

Terdapat pula metode untuk memperoleh y2y_2 ketika λ1=λ2\lambda_1=\lambda_2 yang tidak memerlukan bentuk dengan ln(x)\ln(x). Untuk memperoleh solusi deret pertama y1y_1, kita substitusikan deret

y(x,λ)=xλn=0cnxny(x,\lambda) = x^\lambda \sum_{n=0}^\infty c_n x^n ke dalam (5.1.2) untuk memperoleh bentuk

c0F(λ)xλ+xλn=1(cnF(n+λ)+G(λ,cnp,cnp+1,,cn1))xn=0,c_0 F(\lambda) x^\lambda + x^\lambda \sum_{n=1}^\infty \left(c_n F(n+\lambda) + G(\lambda,c_{n-p},c_{n-p+1},\ldots,c_{n-1})\right) x^n = 0 \ , dengan GG suatu fungsi, p>0p>0 suatu bilangan bulat tetap, dan cj=0c_j = 0 untuk j<0j<0. Jika λ=λ1\lambda=\lambda_1, maka c0c_0 bebas dan relasi rekurensi

cn=G(λ1,cnp,cnp+1,,cn1)F(n+λ1)c_n = -\frac{G(\lambda_1,c_{n-p},c_{n-p+1},\ldots,c_{n-1})} {F(n+\lambda_1)} untuk n1n\geq 1 memberikan semua koefisien solusi lainnya. Ketika λ1=λ2\lambda_1=\lambda_2, persamaan indisialnya adalah F(λ)=(λλ1)2\displaystyle F(\lambda) = (\lambda-\lambda_1)^2. Jika kita substitusikan (5.1.4) dengan

cn=G(λ,cnp,cnp+1,,cn1)F(n+λ)c_n = -\frac{G(\lambda,c_{n-p},c_{n-p+1},\ldots,c_{n-1})}{F(n+\lambda)} untuk n1n \geq 1 ke dalam (5.1.2), maka dari (5.1.5) kita memperoleh

x2(d2ydx2(x,λ)+P(x)dydx(x,λ)+Q(x)y(x,λ))=c0(λλ1)2xλ.x^2\left( \frac{\mathrm{d}^{2}{y}}{\mathrm{d}{x}^{2}}(x,\lambda) + P(x) \frac{\mathrm{d}y}{\mathrm{d}{x}}(x,\lambda) + Q(x) y(x,\lambda)\right) = c_0(\lambda-\lambda_1)^2x^{\lambda} \ . Jika kita mendiferensialkan hubungan ini terhadap λ\lambda, maka kita memperoleh

x2(d2dx2(yλ(x,λ))+P(x)d(yλ(x,λ))dx+Q(x)(yλ(x,λ)))=c0(2(λλ1)+(λλ1)2ln(x))xλ.\begin{aligned} & x^2\left( \frac{\mathrm{d}^{2}}{\mathrm{d}{x}^{2}}{\left( \frac{\partial{y}}{\partial \lambda}(x,\lambda) \right)} + P(x) \frac{\mathrm{d}\left( \frac{\partial{y}}{\partial \lambda}(x,\lambda) \right)}{\mathrm{d}{x}} + Q(x) \left( \frac{\partial{y}}{\partial \lambda}(x,\lambda) \right)\right) \\ &\qquad = c_0\left(2(\lambda-\lambda_1) + (\lambda-\lambda_1)^2\ln(x)\right) x^{\lambda} \ . \end{aligned} Dengan demikian, yλ(x,λ)|λ=λ1\displaystyle\frac{\partial{y}}{\partial \lambda}(x,\lambda)\big|_{\lambda=\lambda_1} merupakan solusi lain dari (5.1.2) yang bebas secara linear dari y1y_1.

Contoh. Selesaikan persamaan diferensial biasa

4xy+2y+y=0.4xy''+ 2 y' + y = 0 \ .

Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=12x\displaystyle P(x) = \frac{1}{2x} dan Q(x)=14x\displaystyle Q(x) = \frac{1}{4x}. Jadi, xP(x)=12\displaystyle x P(x) = \frac{1}{2} dan x2Q(x)=x4\displaystyle x^2 Q(x) = \frac{x}{4}. Kita memperoleh α0=1/2\alpha_0 = 1/2 dan β0=0\beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ212λ=λ(λ12)=0\displaystyle\lambda^2 + (\alpha_0 - 1)\lambda + \beta_0 = \lambda^2 -\frac{1}{2}\,\lambda = \lambda \left(\lambda -\frac{1}{2}\right) = 0. Akar-akarnya adalah λ1=12\displaystyle\lambda_1 = \frac{1}{2} dan λ2=0\lambda_2=0. Karena λ1λ2\lambda_1-\lambda_2 bukan bilangan bulat, kita memperoleh dua solusi yang bebas secara linear, y1=xλ1n=0anxn\displaystyle y_1 = x^{\lambda_1} \sum_{n=0}^\infty a_n x^n dan y2=xλ2n=0bnxn\displaystyle y_2 = x^{\lambda_2} \sum_{n=0}^\infty b_n x^n.

Jika kita substitusikan y=xλn=0cnxn\displaystyle y = x^\lambda \sum_{n=0}^\infty c_n x^n ke dalam 4xy+2y+y=0\displaystyle 4xy''+ 2 y' + y = 0, maka kita memperoleh

0=4xy+2y+y=4n=0(λ+n)(λ+n1)cnxλ+n1+2n=0(λ+n)cnxλ+n1+n=0cnxλ+n=4n=1(λ+n+1)(λ+n)cn+1xλ+n+2n=1(λ+n+1)cn+1xλ+n+n=0cnxλ+n=2λ(2λ1)c0xλ1+n=0(2(λ+n+1)(2λ+2n+1)cn+1+cn)xλ+n.\begin{aligned} 0 & = 4xy''+ 2 y' + y = 4 \sum_{n=0}^\infty (\lambda+n)(\lambda+n-1) c_n x^{\lambda+n-1} + 2 \sum_{n=0}^\infty (\lambda + n )c_n x^{\lambda+n-1} + \sum_{n=0}^\infty c_n x^{\lambda+n} \\ &= 4 \sum_{n=-1}^\infty (\lambda+n+1)(\lambda+n) c_{n+1} x^{\lambda+n} + 2 \sum_{n=-1}^\infty (\lambda + n+1 )c_{n+1} x^{\lambda+n} + \sum_{n=0}^\infty c_n x^{\lambda+n} \\ &= 2\lambda(2\lambda-1) c_0 x^{\lambda-1} + \sum_{n=0}^\infty \big( 2(\lambda+n+1)(2\lambda+2n +1)c_{n+1} +c_n\big) x^{\lambda+n} \ . \end{aligned} Kita memperoleh 2λ(2λ1)c0=02\lambda( 2\lambda -1)c_0 = 0 dan relasi rekurensi cn+1=12(λ+n+1)(2λ+2n+1)cn\displaystyle c_{n+1} = \frac{-1}{2(\lambda+n+1)(2\lambda+2n+1)} c_n untuk n0n\geq 0.

Untuk λ=1/2\lambda = 1/2, kita mempunyai 2λ(2λ1)c0=0\displaystyle 2\lambda( 2\lambda -1)c_0 = 0 untuk setiap nilai c0c_0. Jadi, c0c_0 bebas. Relasi rekurensinya menjadi cn+1=1(2n+3)(2n+2)cn\displaystyle c_{n+1} = \frac{-1}{(2n+3)(2n+2)} c_n untuk n0n \geq 0. Kita memperoleh c1=c03×2=c03!\displaystyle c_1 = -\frac{c_0}{3\times 2} = - \frac{c_0}{3!}, c2=c15×4=c05×4×3×2=c05!\displaystyle c_2 = - \frac{c_1}{5\times 4} = \frac{c_0}{5\times 4 \times 3 \times 2} = \frac{c_0}{5!}, c3=c27×6=c07×6×5×4×3×2=c07!\displaystyle c_3 = - \frac{c_2}{7\times 6} = -\frac{c_0}{7\times 6\times 5\times 4 \times 3 \times 2} = - \frac{c_0}{7!}, dan seterusnya. Oleh induksi, kita memperoleh cn=(1)nc0(2n+1)!\displaystyle c_n = (-1)^n \frac{c_0}{(2n+1)!} untuk n0n \geq 0. Solusi pertamanya adalah

y1=c0x1/2n=0(1)nxn(2n+1)!=c0n=0(1)n(x)2n+1(2n+1)!=c0sin(x).y_1 = c_0 \, x^{1/2} \sum_{n=0}^\infty (-1)^n \frac{x^n}{(2n+1)!} = c_0 \, \sum_{n=0}^\infty (-1)^n \frac{(\sqrt{x})^{2n+1}}{(2n+1)!} = c_0 \sin(\sqrt{x}) \ .

Untuk λ=0\lambda = 0, kita mempunyai 2λ(2λ1)c0=0\displaystyle 2\lambda( 2\lambda -1)c_0 = 0 untuk setiap nilai c0c_0. Jadi, c0c_0 juga bebas. Relasi rekurensinya menjadi cn+1=1(2n+2)(2n+1)cn\displaystyle c_{n+1} = \frac{-1}{(2n+2)(2n+1)} c_n untuk n0n \geq 0. Kita memperoleh c1=c02×1=c02!\displaystyle c_1 = - \frac{c_0}{2\times 1} = - \frac{c_0}{2!}, c2=c14×3=c04×3×2=c04!\displaystyle c_2 = - \frac{c_1}{4\times 3} = \frac{c_0}{4 \times 3 \times 2} = \frac{c_0}{4!}, c3=c26×5=c06×5×4×3×2=c06!\displaystyle c_3 = - \frac{c_2}{6\times 5} = -\frac{c_0}{6\times 5\times 4 \times 3 \times 2} = - \frac{c_0}{6!}, dan seterusnya. Oleh induksi, kita memperoleh cn=(1)nc0(2n)!\displaystyle c_n = (-1)^n \frac{c_0}{(2n)!} untuk n0n\geq 0. Solusi kedua yang bebas secara linear dari solusi pertama adalah

y2=c0n=0(1)nxn(2n)!=c0n=0(1)n(x)2n(2n)!=c0cos(x).y_2 = c_0 \, \sum_{n=0}^\infty (-1)^n \frac{x^n}{(2n)!} = c_0 \, \sum_{n=0}^\infty (-1)^n \frac{(\sqrt{x})^{2n}}{(2n)!} = c_0 \cos(\sqrt{x}) \ .

Solusi umumnya adalah y=Csin(x)+Dcos(x)\displaystyle y = C\, \sin(\sqrt{x}) + D\, \cos(\sqrt{x}) untuk konstanta sembarang CC dan DD.

Contoh. Selesaikan persamaan diferensial biasa

2xy+y+xy=0.2xy''+ y' + x y = 0 \ .

Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=12x\displaystyle P(x) = \frac{1}{2x} dan Q(x)=12\displaystyle Q(x) = \frac{1}{2}. Dengan demikian, xP(x)=12\displaystyle x P(x) = \frac{1}{2} dan x2Q(x)=x22\displaystyle x^2 Q(x) = \frac{x^2}{2}. Kita memperoleh α0=1/2\alpha_0 = 1/2 dan β0=0\beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ212λ=λ(λ12)=0\displaystyle\lambda^2 + (\alpha_0 - 1)\lambda + \beta_0 = \lambda^2 -\frac{1}{2}\,\lambda = \lambda \left(\lambda -\frac{1}{2}\right) = 0. Akar-akarnya adalah λ1=1/2\displaystyle\lambda_1 = 1/2 dan λ2=0\lambda_2=0. Karena λ1λ2\lambda_1-\lambda_2 bukan bilangan bulat, kita mempunyai dua solusi yang bebas secara linear, y1=xλ1n=0anxn\displaystyle y_1 = x^{\lambda_1} \sum_{n=0}^\infty a_n x^n dan y2=xλ2n=0bnxn\displaystyle y_2 = x^{\lambda_2} \sum_{n=0}^\infty b_n x^n.

Jika kita substitusikan y=xλn=0cnxn\displaystyle y = x^\lambda \sum_{n=0}^\infty c_n x^n ke dalam 2xy+y+xy=0\displaystyle 2xy''+ y' + xy = 0, kita memperoleh

0=2xy+y+xy=2n=0(λ+n)(λ+n1)cnxλ+n1+n=0(λ+n)cnxλ+n1+n=0cnxλ+n+1=2n=2(λ+n+2)(λ+n+1)cn+2xλ+n+1+n=2(λ+n+2)cn+2xλ+n+1+n=0cnxλ+n+1=λ(2λ1)c0xλ1+(2λ+1)(λ+1)c1xλ+n=0((λ+n+2)(2λ+2n+3)cn+2+cn)xλ+n+1.\begin{aligned} 0 & = 2xy''+ y' + xy \\ &= 2 \sum_{n=0}^\infty (\lambda+n)(\lambda+n-1) c_n x^{\lambda+n-1} + \sum_{n=0}^\infty (\lambda + n )c_n x^{\lambda+n-1} + \sum_{n=0}^\infty c_n x^{\lambda+n+1} \\ &= 2 \sum_{n=-2}^\infty (\lambda+n+2)(\lambda+n+1) c_{n+2} x^{\lambda+n+1} + \sum_{n=-2}^\infty (\lambda + n+2 )c_{n+2} x^{\lambda+n+1} + \sum_{n=0}^\infty c_n x^{\lambda+n+1} \\ &= \lambda(2\lambda-1)c_0 x^{\lambda-1} + (2\lambda+1)(\lambda+1) c_1 x^{\lambda} \\ &\qquad + \sum_{n=0}^\infty \big( (\lambda+n+2)(2\lambda+2n+3) c_{n+2} +c_n\big) x^{\lambda+n+1} \ . \end{aligned} Kita memperoleh λ(2λ1)c0=0\lambda(2\lambda -1)c_0 = 0, (2λ+1)(λ+1)c1=0(2\lambda+1)(\lambda+1) c_1 = 0 dan hubungan rekurensi

cn+2=1(λ+n+2)(2λ+2n+3)cnc_{n+2} = -\frac{1}{(\lambda+n+2)(2\lambda+2n+3)} c_n untuk n0n \geq 0.

Untuk λ=1/2\lambda = 1/2, kita mempunyai λ(2λ1)c0=0\lambda( 2\lambda -1)c_0 = 0 untuk setiap nilai c0c_0, dan (2λ+1)(λ+1)c1=3c1=0(2\lambda+1)(\lambda+1) c_1 = 3 c_1 = 0. Dengan demikian, c0c_0 bebas dan c1=0c_1=0. Hubungan rekurensinya menjadi cn+2=1(n+2)(2n+5)cn\displaystyle c_{n+2} = -\frac{1}{(n+2)(2n+5)} c_n untuk n0n \geq 0. Kita memperoleh c3=13×7c1=0\displaystyle c_3 = \frac{-1}{3\times 7} c_1 = 0, c5=15×11c3=0\displaystyle c_5 = \frac{-1}{5\times 11} c_3 = 0, dan seterusnya. Dengan demikian, dengan induksi kita memperoleh c2n1=0c_{2n-1} = 0 untuk n>0n>0. Kita juga memperoleh bahwa c2=12×5c0\displaystyle c_2 = -\frac{1}{2\times 5}c_0, c4=14×9c2=c02×4×5×9\displaystyle c_4 = - \frac{1}{4\times 9}c_2 = \frac{c_0}{2\times 4 \times 5 \times 9}, c6=16×13c4=c02×4×6×5×9×13\displaystyle c_6 = - \frac{1}{6\times 13}c_4 = -\frac{c_0}{2\times 4\times 6\times 5 \times 9 \times 13}, dan seterusnya. Dengan induksi, kita memperoleh

c2n=(1)nc02nn!×5×9×13××(4n+1)c_{2n} = (-1)^n \frac{c_0}{2^n n! \times 5 \times 9 \times 13 \times \ldots \times (4n+1)} untuk n>0n >0. Dengan c0=1c_0=1, kita memperoleh solusi pertama

y1=x1/2n=0(1)n12nn!×5×9×13××(4n+1)x2n.y_1 = x^{1/2} \sum_{n=0}^\infty (-1)^n \frac{1}{2^n n! \times 5 \times 9 \times 13 \times \ldots \times (4n+1)}x^{2n} \ .

Untuk λ=0\lambda = 0, kita mempunyai λ(2λ1)c0=0\lambda( 2\lambda -1)c_0 = 0 untuk setiap nilai c0c_0, dan (2λ+1)(λ+1)c1=c1=0(2\lambda+1)(\lambda+1) c_1 = c_1 = 0. Dengan demikian, c0c_0 bebas dan c1=0c_1=0. Hubungan rekurensinya menjadi cn+2=1(n+2)(2n+3)cn\displaystyle c_{n+2} = -\frac{1}{(n+2)(2n+3)} c_n untuk n0n \geq 0. Kita memperoleh c3=13×5c1=0\displaystyle c_3 = \frac{-1}{3\times 5} c_1 = 0, c5=15×9c3=0\displaystyle c_5 = \frac{-1}{5\times 9} c_3 = 0, dan seterusnya. Dengan demikian, dengan induksi kita memperoleh c2n1=0c_{2n-1} = 0 untuk n>0n>0. Kita juga memperoleh bahwa c2=c02×3\displaystyle c_2 = - \frac{c_0}{2\times 3}, c4=14×7c2=c02×4×3×7\displaystyle c_4 = - \frac{1}{4\times 7}c_2 = \frac{c_0}{2\times 4 \times 3 \times 7}, c6=16×11c4=c02×4×6×3×7×11\displaystyle c_6 = - \frac{1}{6\times 11}c_4 = -\frac{c_0}{2\times 4\times 6\times 3 \times 7 \times 11}, dan seterusnya. Dengan induksi, kita memperoleh

c2n=(1)nc02nn!×3×7×11××(4n1)c_{2n} = (-1)^n \frac{c_0}{2^n n! \times 3 \times 7 \times 11 \times \ldots \times (4n-1)} untuk n>0n > 0. Dengan c0=1c_0=1, kita memperoleh solusi kedua yang bebas secara linear dari solusi pertama

y2=1+n=1(1)n12nn!×3×7×11××(4n1)x2n.y_2 = 1+ \sum_{n=1}^\infty (-1)^n \frac{1}{2^n n! \times 3 \times 7 \times 11 \times \ldots \times (4n-1)}x^{2n} \ .

Solusi umumnya adalah y=Cy1(x)+Dy2(x)\displaystyle y = C\, y_1(x) + D\, y_2(x) untuk konstanta sembarang CC dan DD.

Contoh. Selesaikan persamaan diferensial biasa yy=0\displaystyle y''- y = 0.

Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=0\displaystyle P(x) = 0 dan Q(x)=1\displaystyle Q(x) = -1. Dengan demikian, xP(x)=0\displaystyle x P(x) = 0 dan x2Q(x)=x2\displaystyle x^2 Q(x) = -x^2. Kita memperoleh α0=β0=0\alpha_0 = \beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ2λ=λ(λ1)=0\displaystyle\lambda^2 + (\alpha_0 - 1) \lambda + \beta_0 = \lambda^2 - \lambda = \lambda(\lambda -1) = 0. Akar-akarnya adalah λ1=1\displaystyle\lambda_1 = 1 dan λ2=0\lambda_2=0.

Karena λ1λ2\lambda_1-\lambda_2 merupakan bilangan bulat, salah satu solusi berbentuk y1=xn=0anxn=n=0anxn+1\displaystyle y_1 = x \sum_{n=0}^\infty a_n x^n = \sum_{n=0}^\infty a_n x^{n+1} dan solusi lain yang bebas secara linear dari y1y_1 berbentuk y2=ky1ln(x)+n=0bnxn\displaystyle y_2 = k y_1 \ln(x) + \sum_{n=0}^\infty b_n x^n. Namun, seperti akan kita lihat, kita berada dalam keadaan yang memungkinkan kita mengambil k=0k=0 dalam y2y_2 karena b0b_0 dan b1b_1 bebas. Dalam kasus ini, kita mencari solusi berbentuk y=n=0cnxn\displaystyle y = \sum_{n=0}^\infty c_n x^n. Ini akan menghasilkan dua solusi yang bebas secara linear. Jika kita substitusikan bentuk yy ini ke dalam yy=0\displaystyle y''- y = 0, kita memperoleh

yy=n=2n(n1)cnxn2n=0cnxn=n=0(n+1)(n+2)cn+2xnn=0cnxn=n=0((n+1)(n+2)cn+2cn)xn=0.\begin{aligned} y''- y & = \sum_{n=2}^\infty n(n-1) c_n x^{n-2} - \sum_{n=0}^\infty c_n x^n = \sum_{n=0}^\infty (n+1)(n+2) c_{n+2} x^n - \sum_{n=0}^\infty c_n x^n \\ &= \sum_{n=0}^\infty \left( (n+1)(n+2)c_{n+2} - c_n\right) x^n = 0 \ . \end{aligned} Kita memperoleh hubungan rekurensi cn+2=cn(n+1)(n+2)\displaystyle c_{n+2} = \frac{c_n}{(n+1)(n+2)} untuk n0n \geq 0, dengan c0c_0 dan c1c_1 bebas. Kita mempunyai c2=c02=c02!\displaystyle c_2 = \frac{c_0}{2} = \frac{c_0}{2!}, c4=c23×4=c02×3×4=c04!\displaystyle c_4 = \frac{c_2}{3\times 4} = \frac{c_0}{2\times 3 \times 4} = \frac{c_0}{4!}, c6=c45×6=c02×3×4×5×6=c06!\displaystyle c_6 = \frac{c_4}{5\times 6} = \frac{c_0}{2\times 3 \times 4 \times 5 \times 6} = \frac{c_0}{6!}, dan seterusnya. Dengan induksi, kita memperoleh cn=c0n!\displaystyle c_n = \frac{c_0}{n!} untuk nn genap. Selain itu, c3=c12×3=c13!\displaystyle c_3 = \frac{c_1}{2\times 3} = \frac{c_1}{3!}, c5=c34×5=c12×3×4×5=c15!\displaystyle c_5 = \frac{c_3}{4\times 5} = \frac{c_1}{2\times 3\times 4 \times 5} = \frac{c_1}{5!}, c7=c56×7=c12×3×4×5×6×7=c17!\displaystyle c_7 = \frac{c_5}{6\times 7} = \frac{c_1}{2\times 3\times 4 \times 5 \times 6 \times 7} = \frac{c_1}{7!}, dan seterusnya. Dengan induksi, kita memperoleh cn=c1n!\displaystyle c_n = \frac{c_1}{n!} untuk nn ganjil.

Dengan memilih dua pasangan (c0,c1)(c_0,c_1) yang bebas secara linear, kita memperoleh dua solusi yang bebas secara linear. Pilihan (c0,c1)=(0,1)(c_0,c_1) = (0,1) menghasilkan solusi berbentuk y1y_1 di atas, dan pilihan (c0,c1)=(1,0)(c_0,c_1)= (1,0) menghasilkan solusi berbentuk y2y_2 di atas dengan k=0k=0. Namun, kita tidak akan memilih pasangan-pasangan tersebut. Dengan c0=c1=1c_0=c_1 =1, kita memperoleh solusi

y1=n=0x2n(2n)!+n=0x2n+1(2n+1)!=n=0xnn!=ex,y_1 = \sum_{n=0}^\infty \frac{x^{2n}}{(2n)!} + \sum_{n=0}^\infty \frac{x^{2n+1}}{(2n+1)!} = \sum_{n=0}^\infty \frac{x^{n}}{n!} = e^x \ , dan dengan c0=c1=1c_0 = -c_1 = 1, kita memperoleh solusi

y2=n=0x2n(2n)!n=0x2n+1(2n+1)!=n=0(1)nxnn!=ex.y_2 = \sum_{n=0}^\infty \frac{x^{2n}}{(2n)!} - \sum_{n=0}^\infty \frac{x^{2n+1}}{(2n+1)!} = \sum_{n=0}^\infty \frac{(-1)^nx^{n}}{n!} = e^{-x} \ . Kedua solusi tersebut adalah solusi yang biasanya diperoleh dengan metode baku untuk menyelesaikan persamaan diferensial orde dua dengan koefisien konstan. Solusi umumnya adalah y(x)=Aex+Bex\displaystyle y(x) = A e^x + B e^{-x} untuk konstanta sembarang AA dan BB.

Contoh. Selesaikan persamaan diferensial biasa

(1x2)y2xy+2y=0.(1-x^2)y''- 2x y' + 2y = 0 \ .

Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=2x1x2=2xn=0x2n\displaystyle P(x) = \frac{-2x}{1-x^2} = -2x \sum_{n=0}^\infty x^{2n} dan Q(x)=21x2=2n=0x2n\displaystyle Q(x) = \frac{2}{1-x^2} = 2 \sum_{n=0}^\infty x^{2n} untuk |x|<1|x|<1. Dengan demikian, xP(x)=2n=0x2n+2\displaystyle x\,P(x) = -2 \sum_{n=0}^\infty x^{2n+2} dan x2Q(x)=2n=0x2n+2\displaystyle x^2\,Q(x) = 2\sum_{n=0}^\infty x^{2n+2} untuk |x|<1|x|<1. Kita memperoleh α0=β0=0\alpha_0 = \beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ2λ=λ(λ1)=0\displaystyle\lambda^2 + (\alpha_0 - 1)\lambda + \beta_0 = \lambda^2 - \lambda = \lambda (\lambda -1) = 0. Akar-akarnya adalah λ1=1\displaystyle\lambda_1 = 1 dan λ2=0\lambda_2=0. Seperti pada contoh sebelumnya, λ1λ2\lambda_1-\lambda_2 merupakan bilangan bulat dan kita berada dalam keadaan yang memungkinkan kita mengambil k=0k=0 dalam y2y_2 yang diberikan oleh (5.1.3) karena b0b_0 dan b1b_1 bebas, seperti akan kita lihat. Kita mencari solusi berbentuk y=n=0cnxn\displaystyle y = \sum_{n=0}^\infty c_n x^n. Ini akan menghasilkan dua solusi yang bebas secara linear. Jika kita substitusikan bentuk yy ini ke dalam (1x2)y2xy+2y=0\displaystyle(1-x^2)y''- 2x y' + 2y = 0, kita memperoleh

0=(1x2)y2xy+2y=(1x2)n=2n(n1)cnxn22xn=1ncnxn1+2n=0cnxn=n=2n(n1)cnxn2n=2n(n1)cnxn2n=1ncnxn+2n=0cnxn=n=0(n+2)(n+1)cn+2xnn=2n(n1)cnxn2n=1ncnxn+2n=0cnxn=2(c0+c2)+6c3x+n=2((n+2)(n+1)cn+2(n+2)(n1)cn)xn.\begin{aligned} 0 & = (1-x^2)y''- 2x y' + 2y = (1-x^2) \sum_{n=2}^\infty n(n-1) c_n x^{n-2} - 2 x \sum_{n=1}^\infty n c_n x^{n-1} + 2 \sum_{n=0}^\infty c_n x^n \\ &= \sum_{n=2}^\infty n(n-1) c_n x^{n-2} - \sum_{n=2}^\infty n(n-1) c_n x^n -2 \sum_{n=1}^\infty n c_n x^n + 2\sum_{n=0}^\infty c_n x^n \\ &= \sum_{n=0}^\infty (n+2)(n+1) c_{n+2} x^n - \sum_{n=2}^\infty n(n-1) c_n x^n -2 \sum_{n=1}^\infty n c_n x^n + 2\sum_{n=0}^\infty c_n x^n \\ &= 2 (c_0+c_2) + 6 c_3 x + \sum_{n=2}^\infty \big( (n+2)(n+1) c_{n+2} - (n+2)(n-1) c_n \big) x^n \ . \end{aligned} c0c_0 dan c1c_1 bebas, c2=c0c_2 = -c_0, c3=0c_3 =0, dan kita memperoleh hubungan rekurensi

cn+2=(n+2)(n1)(n+2)(n+1)cn=n1n+1cnc_{n+2} = \frac{(n+2)(n-1)}{(n+2)(n+1)} c_n = \frac{n-1}{n+1} c_n untuk n2n\geq 2. Oleh karena itu, c5=c32=0\displaystyle c_5 = \frac{c_3}{2} = 0, c7=2c53=0\displaystyle c_7 = \frac{2c_5}{3} = 0, dan seterusnya. Dengan demikian, cn=0c_n = 0 untuk setiap n>2n>2 yang ganjil. Selain itu, c4=13c2=13c0\displaystyle c_4 = \frac{1}{3} c_2 = - \frac{1}{3} c_0, c6=35c4=15c0\displaystyle c_6 = \frac{3}{5} c_4 = - \frac{1}{5} c_0, c8=57c6=17c0\displaystyle c_8 = \frac{5}{7} c_6 = - \frac{1}{7} c_0, dan seterusnya. Dengan induksi, kita memperoleh bahwa cn=c0n1\displaystyle c_n = - \frac{c_0}{n-1} untuk setiap n>1n>1 yang genap. Solusi pertama diberikan oleh c1=1c_1=1 dan c0=0c_0=0. Kita memperoleh y1(x)=xy_1(x) = x. Solusi kedua yang bebas secara linear dari solusi pertama diberikan oleh c1=0c_1=0 dan c0=1c_0=1. Kita memperoleh y2=n=0x2n2n1\displaystyle y_2 = \sum_{n=0}^\infty \frac{-x^{2n}}{2n-1} untuk |x|<1|x|<1. Solusi umumnya adalah y=Cx+Dn=0x2n2n1\displaystyle y = C\, x + D\,\sum_{n=0}^\infty \frac{-x^{2n}}{2n-1} untuk |x|<1|x|<1 dan konstanta sembarang CC dan DD.

Contoh. Selesaikan persamaan diferensial biasa

xy+y+xy=0.xy''+ y' + x y = 0 \ .

Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=1x\displaystyle P(x) = \frac{1}{x} dan Q(x)=1\displaystyle Q(x) = 1. Dengan demikian, xP(x)=1\displaystyle x P(x) = 1 dan x2Q(x)=x2\displaystyle x^2 Q(x) = x^2. Diperoleh α0=1\alpha_0 = 1 dan β0=0\beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ2=0\displaystyle\lambda^2 + (\alpha_0 - 1)\lambda + \beta_0 = \lambda^2 = 0. Hanya terdapat satu akar dengan multiplisitas aljabar 22, yaitu λ=0\lambda=0. Salah satu solusi berbentuk y1=n=0anxn\displaystyle y_1 = \sum_{n=0}^\infty a_n x^n dan solusi lain yang bebas secara linear terhadap solusi pertama berbentuk y2=ln(x)y1(x)+n=1bnxn\displaystyle y_2 = \ln(x) y_1(x) + \sum_{n=1}^\infty b_n x^n.

Jika y=n=0anxn\displaystyle y = \sum_{n=0}^\infty a_n x^n disubstitusikan ke dalam xy+y+xy=0\displaystyle xy''+ y' + xy = 0, diperoleh

0=xy+y+xy=n=2n(n1)anxn1+n=1nanxn1+n=0anxn+1=n=0(n+2)(n+1)an+2xn+1+n=1(n+2)an+2xn+1+n=0anxn+1=a1+n=0((n+2)2an+2+an)xn+1.\begin{aligned} 0 & = xy''+ y' + xy = \sum_{n=2}^\infty n(n-1) a_n x^{n-1} + \sum_{n=1}^\infty n a_n x^{n-1} + \sum_{n=0}^\infty a_n x^{n+1} \\ &= \sum_{n=0}^\infty (n+2)(n+1) a_{n+2} x^{n+1} + \sum_{n=-1}^\infty (n+2)a_{n+2} x^{n+1} + \sum_{n=0}^\infty a_n x^{n+1} \\ &= a_1 + \sum_{n=0}^\infty \left( (n+2)^2 a_{n+2} +a_n\right) x^{n+1} \ . \end{aligned} Koefisien a0a_0 bebas dan a1=0a_1=0. Relasi rekurensinya menjadi an+2=1(n+2)2an\displaystyle a_{n+2} = -\frac{1}{(n+2)^2} a_n untuk n0n \geq 0. Diperoleh a3=132a1=0\displaystyle a_3 = \frac{-1}{3^2} a_1 = 0, a5=152a3=0\displaystyle a_5 = \frac{-1}{5^2} a_3 = 0, dan seterusnya. Dengan demikian, berdasarkan induksi, a2n1=0a_{2n-1} = 0 untuk n>0n>0. Selain itu, a2=a022\displaystyle a_2 = \frac{-a_0}{2^2}, a4=142a2=a022×42\displaystyle a_4 = -\frac{1}{4^2}a_2 = \frac{a_0}{2^2\times 4^2}, a6=162a4=a022×42×62\displaystyle a_6 = -\frac{1}{6^2}a_4 = \frac{-a_0}{2^2\times 4^2\times 6^2}, dan seterusnya. Berdasarkan induksi, diperoleh a2n=(1)na0(2nn!)2\displaystyle a_{2n} = (-1)^n \frac{a_0}{(2^n n!)^2} untuk n>0n>0. Dengan memilih a0=1a_0=1, diperoleh solusi pertama

y1=n=0(1)n1(2nn!)2x2n.y_1 = \sum_{n=0}^\infty (-1)^n \frac{1}{ (2^n n!)^2}x^{2n} \ .

Jika y=ln(x)y1(x)+n=1bnxn\displaystyle y = \ln(x)y_1(x) + \sum_{n=1}^\infty b_n x^n disubstitusikan ke dalam xy+y+xy=0\displaystyle xy''+ y' + xy = 0, diperoleh

0=xy+y+xy=x(ln(x)y1(x)+2y1(x)xy1(x)x2)+n=2n(n1)bnxn1+(ln(x)y1(x)+y1(x)x)+n=1nbnxn1+xln(x)y1(x)+n=1bnxn+1=ln(x)(xy1(x)+y1(x)+xy1(x))=0+2y1(x)+n=0(n+2)(n+1)bn+2xn+1+n=1(n+2)bn+2xn+1+n=1bnxn+1=4n=1(1)nn(2nn!)2x2n1+b1+4b2x+n=1((n+2)2bn+2+bn)xn+1.\begin{aligned} 0 & = xy''+ y' + xy = x\left( \ln(x) y_1''(x) + \frac{2y_1'(x)}{x}- \frac{y_1(x)}{x^2}\right) + \sum_{n=2}^\infty n(n-1) b_n x^{n-1}\\ &\qquad + \left( \ln(x) y_1'(x) + \frac{y_1(x)}{x}\right) + \sum_{n=1}^\infty n b_n x^{n-1} + x \ln(x) y_1(x) + \sum_{n=1}^\infty b_n x^{n+1} \\ &= \ln(x) \underbrace{\left( x y_1''(x) + y_1'(x) + x y_1(x)\right)}_{=0} +2 y_1'(x) + \sum_{n=0}^\infty (n+2)(n+1) b_{n+2} x^{n+1} \\ &\qquad + \sum_{n=-1}^\infty (n+2) b_{n+2} x^{n+1} + \sum_{n=1}^\infty b_n x^{n+1} \\ &=4\sum_{n=1}^\infty (-1)^n \frac{n}{ (2^n n!)^2}x^{2n-1} + b_1 + 4 b_2x + \sum_{n=1}^\infty \big( (n+2)^2 b_{n+2} + b_n\big) x^{n+1} \ . \end{aligned} Jika suku-suku dalam xn\displaystyle x^n dengan nn genap dan suku-suku dalam xn\displaystyle x^n dengan nn ganjil dikelompokkan, diperoleh

0=b1+n=1((2n+1)2b2n+1+b2n1)x2n\begin{aligned} 0 &= b_1 + \sum_{n=1}^\infty \left((2n+1)^2 b_{2n+1}+b_{2n-1}\right) x^{2n} \end{aligned} dan

0=4n=1(1)nn(2nn!)2x2n1+4b2x+n=2((2n)2b2n+b2n2)x2n1=(4b21)x+n=2(4(1)nn(2nn!)2+(2n)2b2n+b2n2)x2n1.\begin{aligned} 0 &= 4\sum_{n=1}^\infty (-1)^n \frac{n}{ (2^n n!)^2}x^{2n-1} + 4b_2 x + \sum_{n=2}^\infty \big( (2n)^2 b_{2n} + b_{2n-2}\big) x^{2n-1} \\ &= (4b_2-1)x + \sum_{n=2}^\infty\left( 4(-1)^n \frac{n}{(2^n n!)^2} + (2n)^2 b_{2n} + b_{2n-2}\right) x^{2n-1} \ . \end{aligned} Dengan demikian, b1=0b_1=0 dan b2n+1=1(2n+1)2b2n1\displaystyle b_{2n+1} = - \frac{1}{(2n+1)^2}b_{2n-1} untuk n>0n >0. Berdasarkan induksi, diperoleh b2n+1=0b_{2n+1} = 0 untuk n0n \geq 0. Selain itu, b2=1/4b_2 = 1/4 dan b2n=1(2n)2(b2n2+(1)n4n(2nn!)2)\displaystyle b_{2n} = -\frac{1}{(2n)^2} \left( b_{2n-2} + (-1)^n \frac{4n}{(2^n n!)^2}\right) untuk n>1n>1. Oleh karena itu,

b4=142(14+8(222!)2)=124(2!)2(1+12)b6=162(b412(233!)2)=126(3!)2(26(3!)262b413)=126(3!)2(26(3!)224(2!)262(1+12)13)=126(3!)2(1+12+13),\begin{aligned} b_4 &= -\frac{1}{4^2} \left( \frac{1}{4} + \frac{8}{(2^2 2!)^2}\right) = -\frac{1}{2^4 (2!)^2}\left( 1 + \frac{1}{2} \right) \\ b_6 &= -\frac{1}{6^2} \left( b_4 - \frac{12}{(2^3 3!)^2}\right) = -\frac{1}{2^6 (3!)^2} \left( \frac{2^6(3!)^2}{6^2}b_4 -\frac{1}{3}\right) \\ &= -\frac{1}{2^6 (3!)^2} \left( -\frac{2^6(3!)^2}{2^4 (2!)^2 6^2}\left( 1 + \frac{1}{2} \right) -\frac{1}{3}\right) = \frac{1}{2^6 (3!)^2} \left( 1 + \frac{1}{2} + \frac{1}{3}\right) \ , \end{aligned} dan seterusnya. Berdasarkan induksi, diperoleh b2n=(1)n122n(n!)2(1+12+13++1n)\displaystyle b_{2n} = \frac{(-1)^{n-1}}{2^{2n}(n!)^2}\left(1+\frac{1}{2} +\frac{1}{3} + \ldots + \frac{1}{n}\right) untuk n>0n > 0.

Solusi kedua yang bebas secara linear terhadap y1y_1 diberikan oleh

y2=ln(x)y1(x)+n=1(1)n122n(n!)2(1+12+13++1n)x2n.y_2 = \ln(x) y_1(x) + \sum_{n=1}^\infty \frac{(-1)^{n-1}}{2^{2n}(n!)^2}\left(1+\frac{1}{2} +\frac{1}{3} + \ldots + \frac{1}{n}\right) x^{2n} \ .

Solusi umumnya adalah y=Cy1(x)+Dy2(x)\displaystyle y = C\, y_1(x) + D\, y_2(x) untuk konstanta sebarang CC dan DD.

Sebenarnya, y1(x)=J0(x)y_1(x) = J_0(x), yaitu “fungsi Bessel jenis pertama berorde nol”. Dengan normalisasi baku fungsi Bessel jenis kedua Y0Y_0, solusi yang dibangun di sini memenuhi

y2(x)=π2Y0(x)+(ln2γ)J0(x),y_2(x)=\frac{\pi}{2}Y_0(x)+(\ln 2-\gamma)J_0(x), dengan γ\gamma konstanta Euler. Jadi, y2y_2 bukan persis Y0Y_0, tetapi berbeda karena normalisasi dan suatu kelipatan J0J_0. Kita akan mengkaji fungsi-fungsi ini pada bagian mendatang.

Contoh. Selesaikan persamaan diferensial biasa

xy+y+2y=0.xy''+ y' + 2 y = 0 \ . Kita mempunyai persamaan berbentuk y+P(x)y+Q(x)y=0\displaystyle y'' + P(x)\, y' + Q(x) y = 0, dengan P(x)=1x\displaystyle P(x) = \frac{1}{x} dan Q(x)=2x\displaystyle Q(x) = \frac{2}{x}. Dengan demikian, xP(x)=1\displaystyle x P(x) = 1 dan x2Q(x)=2x\displaystyle x^2 Q(x) = 2x. Diperoleh α0=1\alpha_0 = 1 dan β0=0\beta_0= 0. Persamaan indisialnya adalah λ2+(α01)λ+β0=λ2=0\displaystyle\lambda^2 + (\alpha_0 - 1)\lambda + \beta_0 = \lambda^2 = 0. Hanya terdapat satu akar dengan multiplisitas aljabar 22, yaitu λ=0\lambda=0.

Kali ini, kita tidak akan menggunakan rumus kedua dalam (5.1.3) untuk memperoleh dua solusi bebas secara linear. Sebagai gantinya, kita akan menggunakan metode yang melibatkan turunan terhadap λ\lambda yang disebutkan pada awal bagian ini.

Kita mencari solusi berbentuk y=n=0cnxn+λ\displaystyle y = \sum_{n=0}^\infty c_n x^{n+\lambda}. Jika ungkapan ini disubstitusikan ke dalam xy+y+2y=0\displaystyle xy''+ y' + 2y = 0, diperoleh

0=xy+y+2y=n=0(λ+n)(λ+n1)cnxλ+n1+n=0(λ+n)cnxλ+n1+2n=0cnxλ+n=n=1(λ+n+1)(λ+n)cn+1xλ+n+n=1(λ+n+1)cn+1xλ+n+2n=0cnxλ+n=λ2c0xλ1+n=0((λ+n+1)2cn+1+2cn)xλ+n.\begin{aligned} 0 & = xy''+ y' + 2y = \sum_{n=0}^\infty (\lambda +n)(\lambda + n-1) c_n x^{\lambda + n-1} + \sum_{n=0}^\infty (\lambda +n) c_n x^{\lambda + n-1} + 2\sum_{n=0}^\infty c_n x^{\lambda + n} \\ &= \sum_{n=-1}^\infty (\lambda +n+1)(\lambda +n) c_{n+1} x^{\lambda +n} + \sum_{n=-1}^\infty (\lambda +n+1)c_{n+1} x^{\lambda +n} + 2\sum_{n=0}^\infty c_n x^{\lambda +n} \\ &= \lambda^2 c_0 x^{\lambda -1} + \sum_{n=0}^\infty \left( (\lambda +n+1)^2 c_{n+1} + 2c_n\right) x^{\lambda +n} \ . \end{aligned} Dari relasi rekurensi cn+1=2cn(λ+n+1)2\displaystyle c_{n+1} = - \frac{2c_n}{(\lambda +n+1)^2} untuk n0n \geq 0, berdasarkan induksi diperoleh

cn=(2)nc0(λ+n)2(λ+n1)2(λ+1)2c_n = \frac{(-2)^n c_0} {(\lambda +n)^2(\lambda +n-1)^2\cdots (\lambda + 1)^2} untuk n1n \geq 1. Dengan nilai-nilai cnc_n ini, fungsi

y(x)=c0xλ+n=1(2)nc0(λ+n)2(λ+n1)2(λ+1)2xn+λy(x) = c_0 x^\lambda + \sum_{n=1}^\infty \frac{(-2)^n c_0} {(\lambda +n)^2(\lambda +n-1)^2\cdots (\lambda + 1)^2} x^{n+\lambda} memenuhi xy+y+2y=λ2c0xλ1\displaystyle xy''+ y' + 2y = \lambda^2 c_0 x^{\lambda -1}. Jika λ=0\lambda = 0, diperoleh solusi

y1(x)=c0+n=1(2)nc0n2(n1)212xn=c0n=0(2)n(n!)2xn.y_1(x) = c_0 + \sum_{n=1}^\infty \frac{(-2)^n c_0} {n^2 (n-1)^2 \cdots 1^2}x^n = c_0 \sum_{n=0}^\infty \frac{(-2)^n} {(n!)^2} x^n \ . Perhatikan bahwa λ2c0=0\displaystyle\lambda^2 c_0 = 0 dapat berlaku dengan c00c_0 \neq 0 karena λ=0\lambda =0; yaitu, c0c_0 bebas.

Untuk memperoleh solusi kedua yang bebas secara linear terhadap solusi pertama, kita meninjau

y(x)=c0xλ+n=1(2)nc0(λ+n)2(λ+n1)2(λ+1)2xn+λ=xλ(c0+n=1(2)nc0(λ+n)2(λ+n1)2(λ+1)2xn).\begin{aligned} y(x) &= c_0 x^\lambda + \sum_{n=1}^\infty \frac{(-2)^n c_0} {(\lambda +n)^2(\lambda +n-1)^2\cdots (\lambda + 1)^2} x^{n+\lambda} \\ &= x^\lambda \left(c_0 + \sum_{n=1}^\infty \frac{(-2)^n c_0} {(\lambda +n)^2(\lambda +n-1)^2\cdots (\lambda + 1)^2} x^n\right) \ . \end{aligned} Solusi kedua ini diberikan oleh

y2(x)=yλ(x,λ)|λ=0=(xλln(x)(c0+n=1(2)nc0(λ+n)2(λ+n1)2(λ+1)2xn))|λ=0+(xλn=1(2)nc0λ((λ+n)2(λ+n1)2(λ+1)2)xn)|λ=0=ln(x)y1(x)+n=1(2)nc0((2(λ+n)12(λ+n1)12(λ+1)1)(λ+n)2(λ+n1)2(λ+1)2)|λ=0xn=ln(x)y1(x)+c0n=1(2)n+1(n!)2(1+12+13++1n)xn\begin{aligned} y_2(x) &= \frac{\partial{y}}{\partial \lambda}(x,\lambda)\big|_{\lambda=0} = \left(x^\lambda\ln(x) \left(c_0 + \sum_{n=1}^\infty \frac{(-2)^n c_0} {(\lambda +n)^2(\lambda +n-1)^2\cdots (\lambda + 1)^2} x^n\right)\right)\bigg|_{\lambda=0} \\ &\qquad + \left(x^\lambda \sum_{n=1}^\infty (-2)^n c_0 \frac{\partial}{\partial \lambda}{\Big((\lambda +n)^{-2}(\lambda +n-1)^{-2}\cdots (\lambda + 1)^{-2}\Big)} x^n\right)\bigg|_{\lambda=0} \\ &= \ln(x) y_1(x) + \sum_{n=1}^\infty (-2)^n c_0 \Big(\left( -2(\lambda +n)^{-1}-2(\lambda +n-1)^{-1} - \ldots -2 (\lambda + 1)^{-1}\right) \\ &\qquad (\lambda +n)^{-2}(\lambda +n-1)^{-2}\cdots (\lambda + 1)^{-2}\Big)\bigg|_{\lambda=0} x^n \\ &= \ln(x) y_1(x) + c_0 \sum_{n=1}^\infty \frac{(-2)^{n+1}}{(n!)^2} \left(1 + \frac{1}{2} + \frac{1}{3} + \ldots + \frac{1}{n}\right) x^n \end{aligned}

Solusi umumnya adalah y=Cy1(x)+Dy2(x)\displaystyle y = C\, y_1(x) + D\, y_2(x) untuk konstanta sebarang CC dan DD.

Atribusi dan hak. Karya sumber oleh Benoit Dionne dilisensikan berdasarkan Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International. Terjemahan dan perubahan yang dicatat menggunakan lisensi komponen yang sama. Produksi terjemahan dan edisi dibantu OpenAI Codex gpt-5.6-sol, Ultra; semua kredit penulis dan kontributor manusia tetap dipertahankan. Ini bukan terbitan atau dukungan resmi Benoit Dionne maupun University of Ottawa.

← Kembali ke Program Matematika