Frobenius Series Solution
Frobenius Series Solution
Assignment
Assignment Topic:
Series solution_ Frobenius method: Elementary discussion of series solution and solution of
Differential Equation using Frobenius method.
Submitted To:
Mst. Rupale Khatun
Assistant Professor
Dept. of Mathematics,
Rajshahi University of Engineering & Technology.
Submitted By:
Roll: 1703001-1703014, 1703016-1703030
Section : A
Dept. of Computer Science and Engineering.
∑∞ 𝑛 2
𝑛=0 𝐶𝑛 𝑥 = 𝐶0 + 𝐶1 𝑥 + 𝐶2 (𝑥) +…………..
there is a number ρ ∈ [0, ∞] (meaning: ρ > 0 and can be infinity) such that
• the power series converges for all x such that |x – x0| < ρ;
• the power series diverges for all x such that |x – x0| > ρ.
Series Solution:
The power series method is used to seek a power series solution to certain differential equations.
In general, such a solution assumes a power series with unknown coefficients, then substitutes that
solution into the differential equation to find a recurrence relation for the coefficients.
The basic idea to finding a series solution to a differential equation is to assume that we can
write the solution as a power series in the form,
∞
y(x)= ∑𝑛=0 𝐶𝑛 (𝑥 − 𝑥0 )𝑛 ………………(1)
Page |2
and then try to determine what the Cn’s need to be. We will only be able to do this if the
point x=x0, is an ordinary point. We will usually say that (1) is a series solution around x=x0.
Why should we try to solve the equation around singular points? If we want to know what happens
around a singular point, why couldn’t we just pick a regular point nearby, find the solution, and
then study the solution near the singular point? Afterall, this solution around the regular point
should converge for points near the singular point. Unfortunately, the above strategy will not give
us what we want. Because as x approaches the singular point, the convergence of the power series
(obtained by solving the equation at a nearby regular point) becomes worse and worse. And letting
x approach the singular point would not give us useful information.
Problem-Solving Strategy:
1)
Definition 1: A function f(z) is said to be analytic in a region R of the complex plane if f(z) has
a derivative at each point of R and if f(z) is single valued.
Corollary: If f(z) is analytic at a point z, then f(z) has continuous derivatives of all order at
the point z.
Page |3
be a complex function. Since x = (z + z̅)/2 and y = (z − z̅)/2i, substituting for x and y gives
𝑑𝑢 𝑑𝑣
= − 𝑑𝑥 (3)
𝑑𝑦
Equations (2, 3) are known as the Cauchy-Riemann equations. They are a necessary condition for
f = u + iv to be analytic.
Theorem:
𝜕2 𝑢 𝜕2 𝑢
If f(z) is analytic, then + =0 (4)
𝜕2 𝑥 𝜕𝑦2
𝜕2 𝑣 𝜕2 𝑣
+ =0 (5)
𝜕2 𝑥 𝜕𝑦2
The real and imaginary parts of an analytic function are harmonic conjugate functions, i.e.,
solutions to Laplace equation and satisfy the Cauchy Riemann equations (2,3).
3.Example:
𝑥−𝑖𝑦
𝑓(𝑥)=𝑥2 +𝑦2 𝑓(𝑧) =𝑧 3 , z = (x +i y)
Page |4
Example:
1
f(x) = where x ∈ ℝ
1−𝑥
In this function all values except x = 1 are ordinary points cause the function is defined by all
those values except 1.
𝑑2𝑦 𝑑𝑦
+ P(x) + Q(x)y = 0 ……………(i)
𝑑𝑥2 𝑑𝑥
A point x = a is said to be an ordinary point of a differential equation (i) if the functions P(x) and
Q(x) are analytic at the point x = a,
i.e. they have power series expansions in non-negative powers of (x-a) with a positive radius of
convergence. Note that in the above definition the coefficient of the highest order derivative term
is taken as 1.
If x = a is not an ordinary point of the equation, then it is termed as a singular point of the equation.
Example:
𝑑2𝑦 1 𝑑𝑦 1
+ + y = 0
𝑑𝑥2 𝑥−1 𝑑𝑥 𝑥
P(x) is undefined for x = 1 and Q(x) is undefined for x = 0. So the differential equation given in
example is defined for all values except 0 and 1.
𝑑2𝑦 𝑑𝑦
𝑥3 + 𝑥 + y = 0
𝑑𝑥2 𝑑𝑥
𝑑2𝑦 1 𝑑𝑦 1
⇒ + + y = 0
𝑑𝑥2 𝑥2 𝑑𝑥 𝑥3
Here we divided both sides of the equation with 𝑥 3 just to form the shape of the differential
equation of the form (i).
Both P(x) and Q(x) are undefined at x = 0. So that differential equation is defined for all values
except 0.
𝑑2 𝑦 𝑑𝑦
+ 𝑃(𝑥) + 𝑄(𝑥)𝑦 = 0 ……….. (i)
𝑑𝑥2 𝑑𝑥
Singular Point: If P(x) or Q(x) or both P(x) and Q(x) are non-analytic at the point x = a, then
x = a is called singular point of equation ( i )
Explanation: The function f (z) = ez/z is analytic throughout the complex plane—for all values
of z—except at the point z = 0, where the series expansion is not defined because it contains the
term 1/z. The series is:
where the factorial symbol (k!) indicates the product of the integers from k down to 1. When the
function is bounded in a neighborhood around a singularity, the function can be redefined at the
point to remove it; hence it is known as a removable singularity. In contrast, the above function
tends to infinity as z approaches 0; thus, it is not bounded and the singularity is not removable (in
this case, it is known as a simple pole).
Example Math:
𝑑 2𝑦 𝑑𝑦
2𝑥 2 2 + 𝑥 + (𝑥 2 + 1)𝑦 = 0
𝑑𝑥 𝑑𝑥
1 (𝑥2 +1)
𝑃(𝑥) = 𝑄(𝑥) =
2𝑥 2𝑥2
Here, as for x = 0 both P(x) and Q(x) is infinite/undefined, it is called a Singular Point.
Page |7
Definition:
A point x = x0 is called an ordinary point of the equation,
If the point x = x0 is not an ordinary point of the differential equation (1), then it is called a singular
point of the differential equation (1). And a singular point x = x0 of the equation (1) is called a
regular singular point of the differential equation (1) if both (x = x0) P(x) and (x = x0)2 Q(x) are
analytic at x = x0.
Example-1:
Let the equation, x2y" + x y' + (x2 – n2)y = 0
Comparing (1) with y" + P(x)y' + Q(x)y = 0, we get, P(x) = 1/x and Q(x) = (x2 – n2) / x2.
Here (x – 0)P(x) = 1 and (x – 0)2 Q(x) = x2 – n2, showing that both (x – 0)P(x) and (x – 0)2 Q(x)
are analytic at x = 0.
Hence, x = 0 is a regular singular point of the equation, x 2y" + x y' + (x2 – n2)y = 0.
Example-2:
Let the equation, 9x(1 – x)y" - 12y' + 4y = 0 ………(1)
4 4
𝑦 ′′ − 𝑦′ + 𝑦 =0
3𝑥(1 − 𝑥 ) 9𝑥(1 − 𝑥 )
4 4
Comparing it with y" + P(x)y' + Q(x)y = 0, we have P(x) = − and Q(x) = .
3𝑥(1−𝑥) 9𝑥 (1−𝑥)
Since P(x) and Q(x) are not both analytic at x = 0, so x = 0 is not ordinary point of equation (1).
−4 4𝑥
Again x P(x) = 3(1−𝑥) and x2Q(x) = 9(1−𝑥) , showing that both P(x) and Q(x) are analytic at x = 0.
Now, from the above equation, there are two functions P(x) and Q(x) by which regular singular
point and irregular singular point can be determined.
If case 1 is finite and case 2 is infinite then it is called irregular singular point.
Definition 2: According to the second definition, t0 is called an irregular singular point of (I) if
there is no number σ>0 such that every solution x(t) grows not faster than |t−t0|-σ as t→t0 . The
case t0=∞ can be reduced to the case t0=0 by the transformation t→t-1.
At regular singular point singularity can be removed. But at irregular singular point singularity can
not be removed.
Example:
x3 y" + x2 y' +3y=0 where, x=0 is a singular point. Prove that, x=0 is irregular singular point.
1
Where, p(x)=
𝑥
3
Q(x)=
𝑥3
For, x0 = 0,
Case 1 : P(x)(x-x0)
1
= (x-0)
𝑥
=1
Case 1 is finite.
Case 2 : Q(x)(x-x0)2
= 3
(x-0)2
𝑥3
3
=x
3
=0
Case 2 is infinite.
[Proved]
P a g e | 10
and the functions F(x) and G(x) are analytic at x=0. Then the following method for solving (1) is
called Frobenius Method
i.e.: y = x k ∑∞
𝑚=0 cm x
m
-----------------------------------(2)
Thus we take
: y = ∑∞
𝑚=0 cm x
m+k
where, c0 ≠0 -------------------------(3)
y’ = ∑∞
𝑚=0 cm (m +k) xm+k-1
And y” = ∑∞
𝑚=0 cm (m +k)(m+k-1)x
m+k-2
-----------------------(4)
Step 3:
Equating to zero the coefficient of the smallest power of x in the identity obtained in step 2
above, we obtain a quadric equation is k. The quadric equation so obtained is called the indicial
equation
Step 4:
We equate to zero the coefficient of general power in the identity obtained in step 2. The
equation so obtained will be called the recurrence relation
Step 6:
If the recurrence relation connects c m and cm-2 then in general determine c 1 by equation to zero
the coefficient of the next higher power. On the other hand, if the recurrence relation connects c m
and cm-1 this step may be omitted.
Step 7:
After getting various coefficients with the help of steps 5 and 6 above, solution is obtained by
substituting these in (2) and (3) above. Other necessary working will be shown in details.
P a g e | 12
Let, k1 and k2 be the roots of the indicial equation. If k 1 and k2 do not differ by an integer. Then
in general, two independent solutions u and v are obtained by putting k=k 1 and k2 in the series
for y. Then the general solution is y= au + b v , where a and b are arbitrary constants.
Where c’0,c’1,….are the new values of the coefficients corresponding to k=k 2. Since k1-k2 is not
an integer , u/v is never constant. Thus u and v are two independent solutions. The general
solution will be,
If the indicial equation has two equal roots k 1=k2=k, we obtain two independent solutions by
substituting this value of k in y and d y/dk.
The indicial equation has double root k if, (a 0 - 1)2 - 4b0 = 0 giving, r = (1-a0)/2
To determine another solution, we use method of variation of parameters, that is, we replace the
constant c in the solution cu by a function w(x) to be determined such that,
v1 = u1w + uw1
𝑢1 𝑥 𝑘−1 [𝑘𝑐0 + (𝑘 + 1) 𝑐1 𝑥 + ⋯ ] 𝑘
= = +⋯
𝑢 𝑥 𝑘[𝑐0 + 𝑐1 𝑥 + ⋯ ] 𝑥
Then,
2𝑘 + 𝑎0
𝑤2 + ( + ⋯ ) 𝑤1 = 0
𝑥
1
𝑤2 + ( + ⋯ ) 𝑤1 = 0
𝑥
Or,
𝑤2 1
= −𝑥 + ⋯
𝑤1
1
Integrating, log 𝑤1 = − log 𝑥 + ⋯ or, 𝑤1 = 𝑥 𝑒 (𝑥)
Expanding the exponential function in powers of x and integrating once more, we obtain the
expression of w always in the following form,
𝑤 = log 𝑥 + 𝑘1 𝑥 + 𝑘2 𝑥 2 + ⋯
Putting this value of w in v(x), we obtained desired another independent solution v(x). Then
general solution will be, y = au + b v, where a and b are arbitrary constants.
P a g e | 14
Case 3: Roots of indicial equation unequal, differing by an integer and making a coefficient
of y infinite.
If indicial equation has two unequal roots k 1 and k2 differing by an integer, and if some of the
coefficients of y become infinite when k=k 2, we modify the form of y by replacing c 0 by do (k-
k2). We then obtain two independent solutions by putting k=k2 in the modified form of y and
dy/dk. The result of putting k=k1 in y gives a numerical multiple of that obtained by putting k=k 2
and hence we reject the solution obtained by putting k=k 1 in y.
Let the roots k1 and k2 of the indicial equation differ by an integer, say k 1=k and k2=k-p, where p
is a positive integer, then as before one solution is,
For k2 it’s may not be possible to determine another independent solution v(x) . In such case, we
determine v(x) by using method outlined in case-2. Thus as before we first obtained,
2𝑘 + 𝑎0
𝑤2 + ( + ⋯ ) 𝑤1 = 0
𝑥
𝑤2 𝑝+ 1
= −( + ⋯)
𝑤1 𝑥
1 𝑟1 𝑟𝑝
𝑤1 = + + ⋯ + + 𝑟𝑝+1 + 𝑟𝑝+2 𝑥 + ⋯
𝑥 𝑝+1 𝑥 𝑝 𝑥
Integrating again,
1
𝑤=− − ⋯ + 𝑟𝑝 𝑙𝑜𝑔 𝑥 + 𝑟𝑝+1 𝑥 + ⋯
𝑝𝑥 𝑝
putting this value of w in, v(x) = u(x) w(x), we obtain another solution v(x) and as usual the
general solution y = au + b v, where a and b are arbitrary constants.
P a g e | 15
Case 4: Roots of indicial equation unequal, differing by an integer and making a coefficient
of y indeterminate.
If the indicial equation has two roots k1 and k2 (say k1<k2) and if one of the coefficients of y
becomes indeterminate when k=k2, the complete solution is given by putting k=k 2 in y, which
then contains two arbitrary constants. The result of putting k=k1 in y merely gives a numerical
multiple of one of the series contained in the first solution. Hence, we reject the solution obtained
by putting k=k1.
Case 1: [1703009]
Rule:
Let k1 and k2 be the roots of the indicial equation. It k1 and k2 do not differ by an integer. Then,
in general, two independent solutions u and v are obtained by putting k=k1 and k2 in the series
for y. Then the general solution is v = au + b v .Where a and b are arbitrary constants.
Example 1: [1703010]
Solve the series: 9x(1-x)y′′ - 12y′ + 4y = 0
Solution:
Given equation is 9x(1-x)y′′ - 12y′ + 4y = 0 ……………. (1)
Dividing by 9x(1-x), equation 1 can be written in standard form
𝑑2 𝑦 4 𝑑𝑦 4𝑦
- 3𝑥(1−𝑥) + 9𝑥(1−𝑥) = 0
𝑑𝑥2 𝑑𝑥
And, 𝑦′′ = ∑∞
𝑚 =0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 ……….(3)
9x(1-x)∑∞
𝑚 =0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 - 12∑∞
𝑚 =0 𝑐𝑚 (𝑘 + 𝑚)𝑥 𝑘+𝑚−1 +
𝛼 𝑘+𝑚
4∑𝑚 = 0 𝑐𝑚 𝑥 =0
Or, 9x ∑∞
𝑚 =0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 - 9𝑥 2 ∑∞
𝑚=0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 −
𝑘+𝑚−2 ∞ 𝑘+𝑚−1 ∞ 𝑘+𝑚
1)𝑥 -12∑𝑚 = 0 𝑐𝑚 (𝑘 + 𝑚)𝑥 + 4∑𝑚 = 0 𝑐𝑚 𝑥 =0
Or, ∑∞
𝑚 =0 𝑐𝑚 {9(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) − 12(𝑘 + 𝑚)}𝑥 𝑘+𝑚−1 +
∑∞
𝑚 =0 𝑐𝑚 {4 − 9(𝑘 + 𝑚)(𝑘 + 𝑚 − 1)}𝑥 𝑘+𝑚 = 0 ………… (4)
Which is an identity in x. Equating to zero the coefficient of the smallest power of x, namely xk-
1
,(7) gives the indicial equation,
Which are unequal and not differing by an integer. To obtain the recurrence relation, we equate
to zero the coefficient of xk+m-1 and obtain,
3cm (k +m)(3k+3m-7)-cm-1[3k+3(m-1)-4]*[3k+3(m-1)+1]=0
3𝑘+3𝑚−2
Or cm = cm-1 ........................................(9)
3(𝑘+𝑚)
𝑐 3𝑘+1
Taking m=1 in (9) gives c1= 30 ∗ ……………….(10)
𝑘+1
P a g e | 17
3𝑘 +4 3𝑘 +4 𝑐0 3𝑘+1
Next, taking m=2 in (9) gives c2= c1=3(𝑘+2) ∗ ∗ …………….(11)
3(𝑘+2) 3 𝑘+1
And so on, Putting these values in (2) i.e., y=x k(c0+c1x+c2x2 +……..) gives
3𝑘+1 (3𝑘+1)(3𝑘+4)
y=c0xk[1+3(𝑘+1)x+ 𝑥 2 +………………..] ……(12)
32
1 1∗4
Putting k=0 and replacing 𝑐0 by b in (12),y=a(1+3x+3∗6 𝑥 2 +………..)=au,say
7
7 8 8∗11
Next, putting k=3 and replacing c0 by b in (12), y=b𝑥 3 (1+10x+10∗13 𝑥 2 +………..)=bv,say.
Example- 2: [1703012]
Solve the following differential equation in series:
2𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + (1 − 𝑥 2 )𝑦 = 0
Solution:
Given 2𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + (1 − 𝑥 2 )𝑦 = 0 (1)
1 1−𝑥2 1
Dividing by 𝑥 2 , (1) takes standard form 𝑦 ′′ − (2𝑥 ) 𝑦 ′ + { 2𝑥2 } y so that both 𝑥𝑃 (𝑥 ) = −(2)
and 𝑥 2 𝑄(𝑥 ) = (1 − 𝑥 2 )/2. Thus both P(x) and Q(x) are analytic at x=0 and so x=0 is a
regular singular point. Let the series solution of (1) be of the form
∞
(𝑦) 𝑛 = ∑𝑚=0 𝐶𝑚 𝑥 𝑘+𝑚 where 𝐶0 ≠ 0 (2)
∞
Therefore, (𝑦) ′ = ∑𝑚=0 𝐶𝑚 (𝑘 + 𝑚)𝑥 𝑘+𝑚−1 and
∞
(𝑦) ′′ = ∑ 𝐶 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 (3)
𝑚=0 𝑚
P a g e | 18
𝑚=0 𝑚=0
or
∞
∞ ∞
− ∑ 𝐶𝑚 𝑥 𝑘+𝑚+2 = 0
𝑚=0
or
∞ ∞
𝑘+𝑚
∑ 𝐶𝑚 {2(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) − (𝑘 + 𝑚) + 1}𝑥 − ∑ 𝐶𝑚 𝑥 𝑘+𝑚+2 = 0
𝑚=0 𝑚=0
or
∞ ∞
∑𝑚=0 𝐶𝑚 {2(𝑘 + 𝑚)2 − 3(𝑘 + 𝑚) + 1}𝑥 𝑘+𝑚 − ∑𝑚=0 𝐶𝑚 𝑥 𝑘+𝑚+2 = 0
∞ ∞
or ∑𝑚=0 𝐶𝑚 (2𝑘 + 2𝑚 − 1)(𝑘 + 𝑚 − 1)𝑥 2 − ∑𝑚=0 𝐶𝑚 𝑥 𝑘+𝑚+2 = 0 (4)
(4) is an identity. Equating to zero the coefficient of the smallest power of x namely 𝑥 𝑘 , (4)
gives the indical equation 𝐶0 (2𝑘 − 1)(𝑘 − 1) = 0 or (2k-1)(k-1)=0 [ ∴ 𝐶0 ≠ 0 ]
[1703013]
Which are unequal and not differing by an integer. To obtain the recurrence
𝑐𝑚−2
Giving 𝑐𝑚 = ------------------------------------(6)
(2𝑘 +2𝑚−1)(𝑘+𝑚−1)
To obtain 𝑐1, we now equate to zero the co-efficient of 𝑥 𝑘+1 and get
1
𝑐1(2k+1)k = 0, so that 𝑐1 = 0 for both roots k = 1 and k = of the indicial
2
We have, 𝑐1 = 𝑐3 = 𝑐5 = ⋯ … … … . = 0 ---------------------------(7)
gives
𝑥2 𝑥4
x = 𝑐0 𝑥 𝑘 [1 + ( + + ⋯ ] --------(10)
𝑘+1)(2𝑘+3 ) (𝑘+1)(𝑘+3)(2𝑘+3)(2𝑘+5)
1
Next putting k = 2 and replacing 𝑐0 by b in (10) gives
1
𝑥2 𝑥4
y = b𝑥 2 [1 + + + ⋯ … ] = bv, say (Ans)
2.3 2.3.4.7
Example-3:
Solve the following differential equation in series.
3𝑥𝑦 ′′ + 2𝑦 ′ + 𝑦 = 0
Solution:
From we substitute equation (1) by equation (2) (3) (4) and we found
3x∑∞
𝑘=0 𝑎𝑘 (𝑚 + 𝑘 )(𝑚 + 𝑘 − 1)𝑥
𝑚+𝑘−2 + 2 ∑∞ 𝑎 (𝑚 + 𝑘 )𝑥 𝑚+𝑘−1 + ∑∞ 𝑎 𝑥 𝑚+𝑘 = 0
𝑘=0 𝑘 𝑘=0 𝑘
➔ 3∑∞
𝑘=0 𝑎𝑘 (𝑚 + 𝑘 )(𝑚 + 𝑘 − 1)𝑥
𝑚+𝑘−1
+ 2∑∞
𝑘=0 𝑎𝑘 (𝑚 + 𝑘 )𝑥
𝑚+𝑘−1
+ ∑∞
𝑘=0 𝑎𝑘 𝑥
𝑚+𝑘
=0
➔ ∑∞
𝑘=0 𝑎𝑘 (𝑚 + 𝑘 )𝑥
𝑚+𝑘−1 [ (
3 𝑚 + 𝑘 − 1) + 2] + ∑∞
𝑘=0 𝑎𝑘 𝑥
𝑚+𝑘
=0
➔ ∑∞
𝑘=0 𝑎𝑘 (𝑚 + 𝑘 )𝑥
𝑚+𝑘−1 [3𝑚 + 3𝑘 − 1] + ∑∞ 𝑎 𝑥 𝑚+𝑘 = 0 -----(5)
𝑘=0 𝑘
i.e. 𝑥 𝑚−1
𝑎0 𝑚(3𝑚 − 1) = 0[𝑘 = 0]
-> m = 0,1/3
y = ∑∞
𝑘=0 𝑎𝑘 𝑥
𝑚+𝑘
= 𝑥 𝑚 (𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + ⋯ )
P a g e | 21
[1703014]
Equating the coefficients of 𝑥 𝑚 ,
𝑎1 (𝑚 + 1)(3𝑚 + 2) + 𝑎0 = 0 [𝑘 = 1]
𝑎0
∴ 𝑎1 = −
(𝑚 + 1)(3𝑚 + 2)
𝑎𝑘+1 (𝑚 + 𝑘 + 1)(3𝑚 + 3𝑘 + 2) + 𝑎𝑘 = 0
𝑎𝑘
∴ 𝑎𝑘+1 =
(𝑚 + 𝑘 + 1)(3𝑚 + 3𝑘 + 2)
and so on.
𝑎0 𝑎 𝑎
For 𝑚 = 0, 𝑎1 = − 2
; 𝑎2 = 200 ; 𝑎3 = − 480
0
1 𝑎0 𝑎 𝑎
For 𝑚 = 3, 𝑎1 = − ; 𝑎2 = 560 ; 𝑎3 = − 1680
0
4
Here, 𝑦 = 𝑥 𝑚 (𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 + ⋯ )
𝑎0 𝑎0 2 𝑎0 3
∴ (𝑦) 𝑚=0 = 𝑥 0 (𝑎0 − 𝑥+ 𝑥 − 𝑥 +⋯)
2 20 480
1 1 2 1 3
∴ (𝑦)𝑚=0 = 𝑎0 (1 − 𝑥 + 𝑥 − 𝑥 + ⋯)
2 20 480
1 𝑎0 𝑎0 2 𝑎0 3
∴ (𝑦)𝑚= 1 = 𝑥 3 (𝑎0 − 𝑥+ 𝑥 − 𝑥 + ⋯)
3 4 56 1680
1 1 4 1 7 1 10
∴ (𝑦) 1 = 𝑎0 (𝑥 3 − 𝑥 3 + 𝑥3 − 𝑥 3 +⋯)
𝑚=
3 4 56 1680
Case-2: [1703016]
Rule:
If the indicial equation has two equal roots (i.e., 𝑘 = 𝑘1 ), we obtain two independent solutions
𝛿𝑦
by substituting this value of 𝑘 in 𝑦 and 𝛿𝑘 .
𝛿𝑦
𝑦 = 𝑐1 (𝑦)𝑚1 + 𝑐2 ( )𝑚1
𝛿𝑚
Example 1:
1
Comparing it with 𝑦" + 𝑃(𝑥)𝑦′ + 𝑄(𝑥)𝑦 = 0, here we find: 𝑃(𝑥 ) = , 𝑄(𝑥 ) = 1; so that
𝑥
𝑥𝑃 (𝑥 ) = 1 and 𝑥 2 𝑄(𝑥 ) = 𝑥 2 ; which are analytic at 𝑥 = 0. So, 𝑥 = 0 is a regular singular point.
Let the series solution of the given equation be of the form,
𝑦 = ∑∞
𝑚=0 𝐶𝑚 𝑥
𝑘+𝑚
… (2)
∴ 𝑦′ = ∑∞
𝑚=0 𝐶𝑚 (𝑘 + 𝑚)𝑥
𝑘+𝑚−1 and 𝑦" = ∑∞
𝑚=0 𝐶𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥
𝑘+𝑚−2 … (3)
Where, 𝐶0 ≠ 0. Putting the above values of 𝑦, 𝑦 ′ , 𝑦" into equation (1) gives,
∞ ∞ ∞
∞ ∞
𝑘+𝑚−1
∴ ∑ 𝐶𝑚 {(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) + (𝑘 + 𝑚)}𝑥 + ∑ 𝐶𝑚 𝑥 𝑘+𝑚+1 = 0
𝑚=0 𝑚=0
∴ ∑∞ 2 𝑘+𝑚−1 + ∑∞
𝑚=0 𝐶𝑚 (𝑘 + 𝑚) 𝑥 𝑚=0 𝐶𝑚 𝑥
𝑘+𝑚+1 = 0 … (4)
Which is an identity. Equating to zero, the coefficient of the smallest power of x, namely 𝑥 𝑘−1
gives the indicial equation,
𝐶0 𝑘 2 = 0
Since, 𝐶0 ≠ 0, 𝑘2 = 0 ∴ 𝑘 = 0, 0
P a g e | 23
[1703017]
For recurrence relation, equating to zero the coefficient of 𝑥 𝑘+𝑚+ 1 in (4) gives
−𝐶𝑚−2
𝐶𝑚 (𝑘 + 𝑚)2 +𝐶𝑚−2 = 0 or, 𝐶𝑚 = … (5)
(𝑘+𝑚) 2
𝐶2 𝐶0
𝐶4 = − =
(𝑘+4) 2 (𝑘+2)2 (𝑘+4) 2
𝐶4 𝐶0
𝐶6 = − =
(𝑘+6) 2 (𝑘+2)2 (𝑘+2) 2 (𝑘+6) 2
To get another independent solution, substituting, (7) into the L.H.S of (1) and simplifying, we
find
𝜕 2 𝑑 2 𝑦 𝑑𝑦 𝐶0 𝜕 2 𝑘
[𝑥 + + 𝑥𝑦] = (𝑘 𝑥 )
𝜕𝑘 𝑑𝑥 2 𝑑𝑥 𝑥 𝜕𝑘
𝑑2 𝑑 𝜕𝑦 𝐶0
∴ [𝑥 2 + + 𝑥] = [2𝑘𝑥 𝑘 + 𝑘 2 𝑥 𝑘 log 𝑥] … (10)
𝑑𝑥2 𝑑𝑥 𝜕𝑘 𝑥
The presence of the factor k in each term on R.H.S of (10) shows that a second solution is
𝜕𝑦
( ) , where y is given by (7) . Differentiating partially w.r.t ‘k’ (7) gives
𝜕𝑘 𝑘=0
P a g e | 24
𝜕𝑦 𝑥2 𝑥4 𝑥6 𝑥2 −2
= 𝐶0 𝑥 𝑘 log 𝑥 {1 − (𝑘+2) 2 + (𝑘+2) 2 (𝑘+4) 2 − (𝑘+2) 2(𝑘+4) 2 (𝑘+6) 2 + ⋯ } + 𝐶0 𝑥 𝑘 {− (𝑘+2) 2 × 𝑘+2 +
𝜕𝑘
𝑥4 2 2 𝑥6 2 2 2
(𝑘+2) 2 (𝑘+4 )2
[− − ]−( [− − − ] + ⋯} … (11)
𝑘+2 𝑘+4 𝑘+2) 2(𝑘+4) 2 (𝑘+6) 2 𝑘+2 𝑘+4 𝑘+6
𝑑 1
[To find 𝑑𝑥 [ (𝑘+2) 2(𝑘+4) 2(𝑘+6) 2 ], we have proceeded as follows,
1
𝑧=[ ]
(𝑘 + 2)2 (𝑘 + 4)2 (𝑘 + 6)2
1 𝑑𝑧 2 2 2
=− − −
𝑧 𝑑𝑘 𝑘+2 𝑘+4 𝑘+6
𝑑𝑧 2 2 2
∴ = 𝑧 (− − − )
𝑑𝑘 𝑘+2 𝑘+4 𝑘+6
1
Where, 𝑧 = [ (𝑘+2) 2(𝑘+4) 2(𝑘+6) 2 ]; etc. Other terms can be similarly differentiated easily]
𝑥2 𝑥4 1
= 𝑏 [𝑢 log 𝑥 + {22 − 22 42 (1 + 2) + ⋯ }]
Example 2: [1703018]
Solve in series (𝑥 − 𝑥 2 )𝑦 ′′ + (1 − 5𝑥 )𝑦 ′ − 4𝑦 = 0.
(1 − 5𝑥 ) ′ 4
𝑦 ′′ + { } 𝑦 − { } 𝑦 = 0.
(𝑥 − 𝑥 2 ) (𝑥 − 𝑥 2 )
(1−5𝑥)
Comparing it with 𝑦 ′′ + 𝑃(𝑥 )𝑦 ′ + 𝑄(𝑥 )𝑦 = 0, here 𝑃 (𝑥 ) = { (𝑥−𝑥2 )} and
4 (1−5𝑥) 4𝑥
𝑄(𝑥 ) = − {( } so that 𝑥𝑃(𝑥 ) = and 𝑥 2 𝑄(𝑥 ) = − ( . Since 𝑥𝑃(𝑥 ) and 𝑥 2 𝑄(𝑥 )
𝑥−𝑥2 ) (1−𝑥) 1−𝑥)
are analytic at 𝑥 = 0, hence 𝑥 = 0 is a regular singular point of (1).
∞ ∞
or
∞
or
∞ ∞
or
∞ ∞
which is an identity. Equating to zero the coefficient of the smallest power of 𝑥, namely 𝑥 𝑘−1 ,
(4) gives the identical equation
And so on. Putting these in (2),i.e. y = xk(C0 + C1x + C2x2 + C3x3 + … ), gives
(𝑘+2) 2 (𝑘+3) 2 (𝑘+4) 2
y = xkC0[1 + (𝑘+1) 2 x + (𝑘+1) 2x2 + (𝑘+1) 2x3 + … ] …(6)
(𝑘+2) 2 2 2 (𝑘+3) 2 2 2 2
+ C0xk [(𝑘+1) 2 x {𝑘+2 –𝑘+1} + x {𝑘+3 – 𝑘+1} + …]
(𝑘+1) 2
=bv, say
Example 3: [1703020]
Solve the equation (𝑥 − 𝑥 2 )𝑦 ′′ + (1 − 𝑥 )𝑦 ′ − 𝑦 = 0 in series near x=0.
Solution:
Given, (𝑥 − 𝑥 2 )𝑦 ′′ + (1 − 𝑥 )𝑦 ′ − 𝑦 = 0 … (1)
Since x=0 is a regular singular point then let the series solution of equation (1) be of form,
∞
y=𝑥 𝑘 (𝐶0 + 𝐶1 𝑥 + 𝐶2 𝑥 2 + ⋯ ) = ∑𝑚=0 𝐶𝑚 𝑥 𝑘+𝑚 , 𝐶0 ≠ 0 … (2)
𝑦 = ∑∞
𝑚=0(𝑘 + 𝑚)𝐶𝑚 𝑥
𝑘+𝑚−1
and 𝑦" = ∑∞
𝑚=0(𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝐶𝑚 𝑥
𝑘+𝑚−2
… (3)
(𝑥 − 𝑥 2 ) ∑∞
𝑚=0(𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝐶𝑚 𝑥
𝑘+𝑚−2 +(1-x) ∑∞ (𝑘 + 𝑚) 𝐶 𝑥 𝑘+𝑚−1 -∑∞
𝑚=0 𝑚 𝑚=0 𝐶𝑚 𝑥
𝑘+𝑚 =0
or, ∑∞
𝑚=0[(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) + (𝑘 + 𝑚)]𝐶𝑚 𝑥
𝑘+𝑚−1 -∑∞ [(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) +
𝑚=0
+(𝑘 + 𝑚)1]𝐶𝑚 𝑥 𝑘+𝑚 =0
or, ∑∞ 2
𝑚=0(𝑘 + 𝑚) 𝐶𝑚 𝑥
𝑘+𝑚−1 ∑∞
- 𝑚=0[(𝑘 + 𝑚)2 + 1]𝐶𝑚 𝑥 𝑘+𝑚 =0…………….(4)
Equation (4) is an identity. Equating to zero the coefficient of the smallest power of x, namely
𝑥 𝑘−1 , we obtain the indicial equation –
𝐶0 𝑘 2 = 0
So that, k= 0, 0 as 𝐶0 ≠ 0
[1703021]
Putting m= 1, 2, 3… in (5) we have
𝑘2 +1
C1 = (k+1)2 C0 … (6)
𝜕𝑦 (𝑘2 +1) ( 𝑘2 +1){(𝑘+1) 2 +1} 2 (𝑘2 +1) {(𝑘+1) 2+1} { (𝑘+2) 2+1}
= C0 xk log x { 1+ x+ x + x3 + ⋯ } +
𝜕𝑥 (𝑘+1) 2 (𝑘+1) 2 (𝑘+2) 2 (𝑘+1) 2(𝑘+2) 2(𝑘+3) 2
2 2.5
= b{ log x ((1+x+ 4 x 2 + 4.9 x 3 + ⋯)+ (-2x-x2-…)}
= b v, say
Case 3: [1703022]
Rule:
If the indicial equation has two unequal roots 𝑘1 and 𝑘2 (say 𝑘1 > 𝑘2) differing by an integer, and
if some of the coefficients of y become infinite when 𝑘 = 𝑘1 .We modify the form y by replacing
𝐶0 by 𝐷0 (𝑘 − 𝑘2 ) . We then obtain two independent solutions by putting 𝑘 = 𝑘2 in the modified
𝜕𝑦
form of y and . The result of putting 𝑘 = 𝑘1gives a numerical multiple of that obtained by
𝜕𝑘
putting 𝑘 = 𝑘2 and hence we reject the solution obtained by putting 𝑘 = 𝑘1in y.
Example-1: [1703023]
Solve in series 𝑥(1 − 𝑥 )𝑦 ′′ − 3𝑥𝑦 ′ − 𝑦 = 0 near x = 0.
Solution:
3𝑥 𝑥
So that 𝑥𝑃 (𝑥 ) = − 1−𝑥 𝑎𝑛𝑑 𝑥 2 𝑄(𝑥 ) = − 1−𝑥 .
Since 𝑥𝑃 (𝑥 ) and 𝑥 2 𝑄(𝑥 ) are both analytic at x = 0, so x = 0 is a regular singular point of (1).
𝑦′ = ∑ 𝑐𝑚 (𝑘 + 𝑚)𝑘+𝑚+1
𝑚=0
P a g e | 30
∝
′′
𝑦 = ∑ 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 … … … … … … … … … . (3)
𝑚=0
Or, ∑∝𝑚=0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−1 − ∑∝𝑚=0 𝑐𝑚 {(𝑘 + 𝑚) 2 + 2(𝑘 + 𝑚) + 1}𝑥 𝑘+𝑚 = 0
𝑐0 𝑘(𝑘 − 1) = 0
Or, 𝑘 (𝑘 − 1) = 0 since 𝑐0 ≠ 0
[1703024]
These are unequal and differ by an integer. Next, to find the recurrence relation we equate to
zero the coefficient of 𝑥 𝑘+𝑚−1 and obtain
( )
Putting m = 1 in (5) gives 𝐶1 = { 𝑘 + 1 ⁄𝑘 } 𝐶0 ………………………. (6)
𝑘+2 𝑘+2
Putting m = 2 in (5) and using (6) gives 𝐶2 = 𝐶1 = 𝐶0……… (7)
𝑘+1 𝑘
𝑘+3 𝑘+3
Putting m = 3 in (5) and using (7) gives 𝐶3 = 𝐶2 = 𝐶0 ……….. (8)
𝑘+2 𝑘
If we put k = 0 in (9), we find that due to presence of the factor k in their denominators, the
coefficients becomes infinite. To remove this difficulty, we write 𝐶0 = 𝑘 do in (9). Then (9)
becomes
𝑦 = 𝐶0 (𝑥 + 2𝑥 + 3𝑥 2 + ⋯ ) ……………………………………………….… (12)
Which is not distinct (i.e. not linearly because ratio of the two series in (11) and (12) is a
constant from (11). Hence (12) will not serve the purpose of a second solution. In such a case the
𝜕𝑦
second independent solution is given by (𝜕𝑥 )𝑘=0 Differentiating (10) partially w.r.t ‘k’
𝜕𝑦
= 𝑑0 𝑥 𝑘𝑙𝑜𝑔𝑥 [𝑘 + (𝑘 + 1)𝑥 + (𝑘 + 2) 𝑥 2 + ⋯ ] + 𝑑0 𝑥 𝑘[1 + 𝑥 + 𝑥 2 + ⋯ ]……….. (13)
𝜕𝑘
P a g e | 32
𝜕𝑦
( ) = 𝑏𝑙𝑜𝑔𝑥(𝑥 + 2𝑥 2 + 3𝑥 3 + ⋯ ) + 𝑏(1 + 𝑥 + 𝑥 2 + ⋯ )
𝜕𝑘 𝑘=0
𝜕𝑦
(𝜕𝑘) 𝑘=0 = 𝑏[𝑢𝑙𝑜𝑔𝑥 + (1 + 𝑥 + 𝑥 2 + ⋯ )] = 𝑏𝑣, 𝑏𝑦(11) …………. (14)
[1703025]
Example: 02
Solve in series the Bessel’s equation of order 2, near x=0, 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + (𝑥 2 − 4)𝑦 = 0.
Solution:
Given 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + (𝑥 2 − 4)𝑦 = 0.
1 (𝑥2 −4)
Dividing by 𝑥 2 , 𝑦 ′′ + (𝑥 ) 𝑦 ′ + [( )]y=0.
𝑥2
1
Comparing it with 𝑦 ′′ + 𝑃(𝑥 )𝑦 ′ + 𝑄(𝑥 )𝑦 = 0. ℎ𝑒𝑟𝑒 𝑃(𝑥 ) = and Q(x) = (𝑥 2 − 4)/𝑥 2 so
𝑥
that xP(x) =1 and 𝑥 2 𝑄(𝑥 ) = 𝑥 2 − 4. 𝑆𝑖𝑛𝑐𝑒 𝑥𝑃(𝑥 )𝑎𝑛𝑑 𝑥 2 𝑄(𝑥 )𝑎𝑟𝑒 𝑏𝑜𝑡ℎ 𝑎𝑛𝑎𝑙𝑦𝑡𝑖𝑐 𝑎𝑡 𝑥 = 0.
Hence x=0 is a regular singular point of (1). Let the series solution of (1) be
∞
y = ∑𝑚=0 𝑐𝑚 𝑥 𝑘+𝑚 , where 𝑐0 ≠ 0.
∞ ∞
y′ = ∑𝑚=0 𝑐𝑚 (𝑘 + 𝑚)𝑥 𝑘+𝑚−1 and y′′ = ∑𝑚=0 𝑐𝑚 (𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝑥 𝑘+𝑚−2 …
(3)
∞
∞
Or ∑ 𝑐𝑚 {(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) + (𝑘 + 𝑚) − 4}𝑥 𝑘+𝑚 + ∑𝑚=0 𝑐𝑚 𝑥 𝑘+𝑚+2 =
𝑚=0
∞ ∞
0 or∑𝑚=0 𝑐𝑚 {(𝑘 + 𝑚)2 − 4}𝑥 𝑘+𝑚 + ∑𝑚=0 𝑐𝑚 𝑥 𝑘+𝑚+2 = 0.
∞ ∞
Or ∑𝑚=0 𝑐𝑚 (𝑘 + 𝑚 + 2)(𝑘 + 𝑚 − 2)𝑥 𝑘+𝑚 + ∑𝑚=0 𝑐𝑚 𝑥 𝑘+𝑚+2 = 0 …. (4)
Which is an identify Equation to zero the coefficient of the smallest power of x. namely 𝑥 𝑘 ,
given the indicial equation 𝑐𝑜 (k+2) (k-2) =0 or (k+2) (k-2) =0 [𝑠𝑖𝑛𝑐𝑒 , 𝑐𝑜 ≠ 0]
[1703026]
This gives k = 2 and k = -2. This are unequal and differ by an integer. For the recurrence relation
we equate to zero the coefficient of 𝑥 𝑘+𝑚 and get
1
𝐶𝑚 (k+m+2) (k+m-2) + 𝐶𝑚−2 = 0 so 𝐶𝑚 = - 𝐶𝑚−2 …………………………….
(𝑘+𝑚+2)(𝑘+𝑚−2)
(5)
To determine𝐶1, we equate to zero the coefficient of xk+1 and get 𝐶1 (k+3) (k-1) = 0 giving C1 = 0
for both the roots k = 2 and k = -2 of the identical equation. Now using C1 = 0 and (5), we get
𝐶1 = 𝐶3 = 𝐶5 = 𝐶7 = ⋯ = 0................................................. (6)
𝑥2 𝑥4 𝑥6
Putting k = 2 in (7) yields, y = 𝐶0 𝑥 2 {1- 2.6 + 2.4.8.6 − 2.4.62 8.10 +……} = 𝐶0 𝜔 , say …………….
(8)
P a g e | 34
Next if we putt k = -2 in (7), the coefficient of 𝑥 4 , 𝑥 6 , becomes infinite. To get rids of this
difficulty, we put 𝐶0 = 𝐷0 (𝑘 + 2) in (7) and obtained the modified solution as,
(𝑘+2)𝑥2 𝑥4
y = 𝐷0 𝑥 𝑘{(𝑘 + 2) − + 𝑘(𝑘+4)(𝑘+6) −
𝑘(𝑘+4)
𝑥6
2 +………}………………………. (9)
𝑘 (𝑘+4) (𝑘+6)(𝑘+8)
Showing that 𝜔 and u are dependent solutions. Hence we must find one more independent solution
in order to obtain the required general solution.
To get another independent solution, substituting (9) into the L.H.S. of (1) and simplifying. We
find 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + (𝑥 2 − 4)𝑦 = 𝐷0 (𝑘 − 2)(𝑘 + 2)2 𝑥 𝑘 ………….. (11)
𝜕 𝑑 2𝑦 𝑑𝑦 𝜕
[ 𝑥2 2 + 𝑥 + (𝑥 2 − 4)𝑦] = 𝐷0 [(𝑘 + 2)2 . 𝑥 𝑘 . (𝑘 − 2)]
𝜕𝑘 𝑑𝑥 𝑑𝑥 𝜕𝑘
𝑑2 𝑑 𝜕𝑦
Or [ 𝑥 2 𝑑𝑥2 + 𝑥 𝑑𝑥 + (𝑥 2 − 4)] 𝜕𝑘 = 2𝐷0 (𝑘 + 2). 𝑥 𝑘 . (𝑘 − 2) + 𝐷0 (𝑘 + 2)2 [𝑥 𝑘 log 𝑥 . (𝑘 −
2) + 𝑥 𝑘 − 1 ]………………………………….. (12)
The presence of the factor (𝑘 + 2) in each term on R.H.S. of (12) shows that a second solution is
𝜕𝑦
(𝜕𝑘) 𝑘=2 where y is given by (9). Differentiating (9) partially w.r.t. ‘k’ we get
𝜕𝑦 (𝑘 + 2)𝑥 2 𝑥4 𝑥6
= 𝐷0 𝑥 𝑘 log 𝑥. {(𝑘 + 2) − + − + ⋯}
𝜕𝑘 𝑘 (𝑘 + 4) 𝑘 (𝑘 + 4)(𝑘 + 6) 𝑘(𝑘 + 4)2 (𝑘 + 6)(𝑘 + 8)
+
(𝑘 + 2) 𝑥 2 1 1 1 𝑥4 1 1 1
𝐷0 𝑥𝑘 {1 − ( − − )+ (− − − )
𝑘 (𝑘 + 4) 𝑘 + 2 𝑘 𝑘 + 4 𝑘 (𝑘 + 4)(𝑘 + 6) 𝑘 𝑘+4 𝑘+6
𝑥6 1 2 1 1
− 2 × (− − − − ) + ⋯… ……… .. }
𝑘 (𝑘 + 4) (𝑘 + 6)(𝑘 + 8) 𝑘 𝑘+4 𝑘+6 𝑘 +8
P a g e | 35
𝑑 𝑘 +2 𝑘 +2 𝑘 +2
[To find𝑑𝑘 {𝑘(𝑘+4)}, we proceed as follows: Take z = 𝑘(𝑘+4) so that log 𝑧 = log[𝑘(𝑘+4)]
𝜕𝑦 𝑥4 𝑥6 𝑥2
= 𝑏𝑥 −2 log{0 − 0𝑥 2 + − + ⋯ } + 𝑏𝑥 −2 {1 +
𝜕𝑘 𝑘=−2 (−2)(2)(4) ( −2)(2) 2 (4)(6) 22
𝑥4 1 𝑥6 1 1 1
+ 2 ( )+ 3 ( −1+ − )+⋯
2 .4 4 2 . 4.6 2 4 6
𝑥 2 𝑥4
= 𝑏 [𝑢 log 𝑥 + 𝑥 −2 {1 + 2 + 2 2 − ⋯ }] = 𝑏𝑣
2 2 .4
by (10)
The required solution is 𝑦 = 𝑎𝑢 + 𝑏𝑣, where ‘a’ and ‘b’ are arbitrary constants.
Case 4: [1703027]
Rule:
If the indicial equation has two roots k1 and k2 (say k1<k2) and if the one of the coefficients of y
becomes indeterminate when k=k2, the complete solution is given by putting k=k 2 in y, which
then contains two arbitrary constants. The result of putting k=k 1 in y merely gives a numerial
multiple of one of the series contained in the first solution. Hence, we reject the solution obtained
by putting k=k1.
P a g e | 36
Example-1: [1703029]
Find solution near x = 0 of 𝑥 2 𝑦 ′′ +(𝑥 + 𝑥 2 )𝑦 ′ +(𝑥 − 9)𝑦 = 0
Solution:
𝑦 ′′ + {(1 + 𝑥 ) ∕ 𝑥 }𝑦 ′ + {(𝑥 − 9) ∕ 𝑥 2 }𝑦 = 0
Since both 𝑥𝑃 (𝑥 ) and 𝑥 2 𝑄(𝑥 ) are analytic at x = 0, hence x = 0 is a regular singular point of (1).
Let the series solution of (1) be
∞
𝑦 = 𝑥 𝑘(𝐶0 + 𝐶1 𝑥 + 𝐶2 𝑥 2 + . . 𝑎𝑑. 𝑖𝑛𝑓. ) = ∑𝑚=0 𝐶𝑚 𝑥 𝑘+𝑚 , where 𝐶0 ≠ 0 …(2)
∞
So 𝑦 ′ = ∑𝑚=0(𝑘 + 𝑚)𝐶𝑚 𝑥 𝑘+𝑚−1
∞
and 𝑦 ′′ = ∑𝑚=0(𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝐶𝑚 𝑥 𝑘+𝑚−2 …(3)
+ (𝑥 − 9) ∑ 𝐶𝑚 𝑥 𝑚+𝑘 = 0
𝑚=0
∞ ∞
Or ∑𝑚=0(𝑘 + 𝑚)(𝑘 + 𝑚 − 1)𝐶𝑚 𝑥 𝑘+𝑚 + ∑𝑚=0(𝑘 + 𝑚) 𝐶𝑚 𝑥 𝑘+𝑚 +
∞ ∞ ∞
∞ ∞
Or ∑𝑚=0{(𝑘 + 𝑚)(𝑘 + 𝑚 − 1) + (𝑘 + 𝑚) − 9}𝐶𝑚 𝑥 𝑘+𝑚 + ∑𝑚=0(𝑘 + 𝑚 + 1)𝐶𝑚 𝑥 𝑘+𝑚+1 = 0
P a g e | 37
∞ ∞
Or ∑𝑚=0 {(𝑘 + 𝑚)2 − 32 } 𝐶𝑚 𝑥 𝑘+𝑚 + ∑𝑚=0(𝑘 + 𝑚 + 1) 𝐶𝑚 𝑥 𝑚+𝑘+1 = 0
∞
Or ∑∞
𝑚=0(𝑘 + 𝑚 + 3)(𝑘 + 𝑚 − 3) 𝐶𝑚 𝑥
𝑘+𝑚
+ ∑𝑚=0 (𝑘 + 𝑚 + 1) 𝐶𝑚 𝑥 𝑚+𝑘+1 = 0…….(4)
(𝑘+2) (𝑘+1)(𝑘+2)
𝐶2 = - (𝑘+5)(𝑘−1) 𝐶1=(𝑘+4)(𝑘−2)(𝑘+5)(𝑘−1) 𝐶0,
−(𝑘+3) (𝑘+1)(𝑘+2)(𝑘+3)
𝐶3 = (𝑘+6)k 𝐶2 = - 𝑘(𝑘+4)(𝑘−2)(𝑘+5)(𝑘−1)(𝑘+6) 𝐶0 ,And so on.
Example-2. [1703030]
Solve in series (1-x2) 𝑦 ′′ -x𝑦 ′ +4y=0
Solution:
Given,
Clearly x=0 is an ordinary point of (1). Let the series solution of (1) be of the form
y′ =∑𝛼𝑚=0 𝐶𝑚 xk+m-1
which is an identity. Equating to zero the coefficient of the smallest power of x, namely xk-2
To get the recurrence relation, we equate to zero the coefficient of x k +m- 2 .Thus
Next, we equate to zero the coefficient of xk-1 and get C1(k+1) k=0 ...(6)
If we take k=0, (6) shows that C1 is indeterminate. With k=0 and using (5), We can express C2,
C4, C6…… in terms of C1 if we assume C1 to be finite. Thus,
Cm ={(m-4)/(m-1)} Cm-2
And C3=(-C1)/ (2) =-C1/2 C5=C3/ (4) =-C1/8, C7=(3C5)/ (6) =-C1/16.
which is the required series solution C0 and C1 being two arbitrary constants.