# First term for particular solution of series solution of inhomogenous DEWatch

Announcements
#1
I am solving a differential equation:
x*y''(x)+(x+1)*y'(x)+x*y(x) = x^2+1
Firstly,I tried to find the complementary solution.
Using Frobenius' method and expanding around x=0,
the first solution is
1-(1/4)*x^2+(1/18)*x^3+(1/192)*x^4-(11/3600)*x^5+(29/103680)*x^6+O(x^7)
the second solution is
ln(x)*(1-(1/4)*x^2+(1/18)*x^3+(1/192)*x^4-(11/3600)*x^5+(29/103680)*x^6+O(x^7))-x+(1/2)*x^2-(1/108)*x^3-(41/1152)*x^4+(1529/216000)*x^5-(1/345600)*x^6+O(x^7)
Since the indicial equation has repeated solutions,both equals to zero.
After that,I tried to find the particular solution.
I assumed the particular solution has the form of a series,
y=c0+c1(x)+c2(x^2)+O(x^3)
and tried to compare the coefficients.
Although I obtained the particular solution
x*(1-(1/4)*x+(1/18)*x^2+(1/192)*x^3-(11/3600)*x^4+(29/103680)*x^5+O(x^6))
,which is correct after checking using Maple,but I assumed c0=0 to obtain this solution,but why?
And also,when finding P(x) in y2=y1(ln x)+P(x)(after substituting the solution of indicial equation)
,I also had to assume that the term without the unknown is 0,is it because the indicial equation has repeated roots?
Sorry for my bad English and my ignorance since I am self-studying in this topic.
By the way,I am self-studying based on this book:
(Undergraduate Lecture Notes in Physics) Lev Kantorovich - Mathematics for Natural Scientists_ Fundamentals and Basics-Springer (2015)
Thanks.
0
4 weeks ago
#2
I doubt many people are going to read this because it looks absolutely horrendous typed out! Can you perhaps upload a picture of it instead?
0
#3
(Original post by Y12_FurtherMaths)
I doubt many people are going to read this because it looks absolutely horrendous typed out! Can you perhaps upload a picture of it instead?
OK...
The complete solution ,y=c0(y1)+c1(y2)+yp
0
#4
(Original post by Y12_FurtherMaths)
I doubt many people are going to read this because it looks absolutely horrendous typed out! Can you perhaps upload a picture of it instead?
Is it OK? Or you need additional information?
0
4 weeks ago
#5
Is it OK? Or you need additional information?
I'm afraid I'm unable to help you with this. Is this university maths?
0
#6
(Original post by Y12_FurtherMaths)
I'm afraid I'm unable to help you with this. Is this university maths?
I don't know,I am self studying for this and I still don't have a uni degree.But I think series solution is not covered in A Level syllabus.
0
4 weeks ago
#7
I don't know,I am self studying for this and I still don't have a uni degree.But I think series solution is not covered in A Level syllabus.
It doesn't seem like people will be able to help you im afraid
0
#8
(Original post by Y12_FurtherMaths)
It doesn't seem like people will be able to help you im afraid
Never mind,nobody is obligated to help me here.Thanks.By the way,in my previous reply,I mean that I am not entering uni yet.
1
4 weeks ago
#9
I am solving a differential equation:
x*y''(x)+(x+1)*y'(x)+x*y(x) = x^2+1
Firstly,I tried to find the complementary solution.
Using Frobenius' method and expanding around x=0,
the first solution is
1-(1/4)*x^2+(1/18)*x^3+(1/192)*x^4-(11/3600)*x^5+(29/103680)*x^6+O(x^7)
the second solution is
ln(x)*(1-(1/4)*x^2+(1/18)*x^3+(1/192)*x^4-(11/3600)*x^5+(29/103680)*x^6+O(x^7))-x+(1/2)*x^2-(1/108)*x^3-(41/1152)*x^4+(1529/216000)*x^5-(1/345600)*x^6+O(x^7)
Since the indicial equation has repeated solutions,both equals to zero.
After that,I tried to find the particular solution.
I assumed the particular solution has the form of a series,
y=c0+c1(x)+c2(x^2)+O(x^3)
and tried to compare the coefficients.
Although I obtained the particular solution
x*(1-(1/4)*x+(1/18)*x^2+(1/192)*x^3-(11/3600)*x^4+(29/103680)*x^5+O(x^6))
,which is correct after checking using Maple,but I assumed c0=0 to obtain this solution,but why?
And also,when finding P(x) in y2=y1(ln x)+P(x)(after substituting the solution of indicial equation)
,I also had to assume that the term without the unknown is 0,is it because the indicial equation has repeated roots?
Sorry for my bad English and my ignorance since I am self-studying in this topic.
By the way,I am self-studying based on this book:
(Undergraduate Lecture Notes in Physics) Lev Kantorovich - Mathematics for Natural Scientists_ Fundamentals and Basics-Springer (2015)
Thanks.
This is university maths. It goes deep into the theory of Linear Differential Equations.

So to answer your question, I'll run you down the problem quickly.

Firstly, you put the ODE into canonical form

.

It is clear that is a regular singular point, therefore Method of Frobenius is the way to go.

Solving the homogeneous case, you found the solution

.

The second solution, because the roots repeat, is

.

If you substitute this into the homogeneous case and compare coefficients, the first two comparisons would yield you

.... (from comparing coefficients of )

and

(*) .... (from comparing coefficients of )

Any further comparisons will not isolate for you, in fact none of them will add on any extra restrictions on it. So in fact, you can choose it arbitrarily as long as equation (*) holds. Of course, whenever we say arbitrarily, we usually mean try to choose something that makes the problem simple. So of course, we can just choose . This implies , and so all of our are now fully defined.

Same approach applies to when you try a particular solution

.

You run into the issue of having no restriction on but you want your to be defined, so just choose it to be zero.

Hence you obtain the solution that Maple is throwing at you

.
Last edited by RDKGames; 4 weeks ago
1
#10
(Original post by RDKGames)
This is university maths. It goes deep into the theory of Linear Differential Equations.

So to answer your question, I'll run you down the problem quickly.

Firstly, you put the ODE into canonical form

.

It is clear that is a regular singular point, therefore Method of Frobenius is the way to go.

Solving the homogeneous case, you found the solution

.

The second solution, because the roots repeat, is

.

If you substitute this into the homogeneous case and compare coefficients, the first two comparisons would yield you

.... (from comparing coefficients of )

and

(*) .... (from comparing coefficients of )

Any further comparisons will not isolate for you, in fact none of them will add on any extra restrictions on it. So in fact, you can choose it arbitrarily as long as equation (*) holds. Of course, whenever we say arbitrarily, we usually mean try to choose something that makes the problem simple. So of course, we can just choose . This implies , and so all of our are now fully defined.

Same approach applies to when you try a particular solution

.

You run into the issue of having no restriction on but you want your to be defined, so just choose it to be zero.

Hence you obtain the solution that Maple is throwing at you

.
Well,do you mean that both the terms can be arbitrarily chosen,so 0 is chosen to make the calculation more easier? If we choose another number other than zero,does the resulted series is just the series using 0 plus multiple of y1?
0
#11
(Original post by RDKGames)
.... (from comparing coefficients of )
And forgive me,is this instead of ?
0
4 weeks ago
#12
Well,do you mean that both the terms can be arbitrarily chosen,so 0 is chosen to make the calculation more easier? If we choose another number other than zero,does the resulted series is just the series using 0 plus multiple of y1?
In the series you find that but ultimately depends on what is. In fact, every for depends on what is. Once you set that, every other coefficient gets fixed in place.

This coefficient is an arbitrary choice for us, because we have no restriction for it. Any suitable value of it will yield every other value and hence give us a solution for . If that makes sense to you, then the next step for you to understand is that is the simplest case for us pick, so just do that.

If you pick a non-zero value then you're going to get the same sum but expressed differently.

You do the same thing when you determine the solution for . You had the choice of picking arbitrary but you picked since the obtained recurrence relation must satisfy

so (from comparing terms) and . Similarly, depends on , and so on. You fixed and obtained all the coefficients you see in .

And forgive me,is this instead of ?
No. Once you substitute in a series into this particular ODE, you should notice that we always begin with , then we have the (constant) terms, and then , and so on... but for my explanation I'm only considering and .

The coefficients turn out to be (when you substitute in ) and (when you substitute in ) which yields us and . But once again I stress the fact that we have no restrictions on hence we choose them arbitrarily so we just pick and .

I say they're arbitrary, for you shouldn't (for instance) pick because then the entire series for is just zero -- a trivial solution.
Last edited by RDKGames; 4 weeks ago
0
#13
(Original post by RDKGames)
In the series you find that but ultimately depends on what is. In fact, every for depends on what is. Once you set that, every other coefficient gets fixed in place.

This coefficient is an arbitrary choice for us, because we have no restriction for it. Any suitable value of it will yield every other value and hence give us a solution for . If that makes sense to you, then the next step for you to understand is that is the simplest case for us pick, so just do that.

If you pick a non-zero value then you're going to get the same sum but expressed differently.

You do the same thing when you determine the solution for . You had the choice of picking arbitrary but you picked since the obtained recurrence relation must satisfy

so (from comparing terms) and . Similarly, depends on , and so on. You fixed and obtained all the coefficients you see in .

No. Once you substitute in a series into this particular ODE, you should notice that we always begin with , then we have the (constant) terms, and then , and so on... but for my explanation I'm only considering and .

The coefficients turn out to be (when you substitute in ) and (when you substitute in ) which yields us and . But once again I stress the fact that we have no restrictions on hence we choose them arbitrarily so we just pick and .

I say they're arbitrary, for you shouldn't (for instance) pick because then the entire series for is just zero -- a trivial solution.
Yes,just like what you had said,during my calculations I found that b1 is determined and b2 depends on b0.OK,after your explainations I accepted the choice of choosing zero.

Well,I think you use x^(-1) and I use x^(0) because you substituted the series to the canonical form of the equation but I substituted it to the original equation in the first post.Am I right?
0
4 weeks ago
#14
Well,I think you use x^(-1) and I use x^(0) because you substituted the series to the canonical form of the equation but I substituted it to the original equation in the first post.Am I right?
Yep.
0
#15
(Original post by RDKGames)
Yep.
0
X

new posts
Back
to top
Latest
My Feed

### Oops, nobody has postedin the last few hours.

Why not re-start the conversation?

see more

### See more of what you like onThe Student Room

You can personalise what you see on TSR. Tell us a little about yourself to get started.

### University open days

• Bournemouth University
Midwifery Open Day at Portsmouth Campus Undergraduate
Wed, 18 Dec '19
• The University of Law
Open Day – GDL and LPC - Chester campus Postgraduate
Sat, 4 Jan '20
• University of East Anglia
Mon, 6 Jan '20

### Poll

Join the discussion

#### Did you vote in the 2019 general election?

Yes (387)
44.38%
No (90)
10.32%
I'm not old enough (395)
45.3%