Re: Series solution to DE

Another attempt using Frobenius method on the original DE.

I reach .

Solving the inidicial equation leads me to , or or . Though one would expect at most 2 different values for c since it's a second order DE. Something's definitely wrong with my attempt.

Any help is appreciated.

Re: Series solution to DE

Quote:

Originally Posted by

**arbolis** Consider the DE

. Give 2 solutions; one regular and worth 1 at the origin and the other of the form

where

and

are regular at the origin. Give the first 3 terms of the series of

and

.

My attempt: Divide the DE by x:

. In order to solve this DE, I had in mind to propose a solution of the form

where

would be the solution to this DE but when x tends to infinity. It turns out that this didn't simplify things as I'd hoped.

When

, the DE becomes

. I used Frobenius's method to solve this DE:

I assumed that

. I derivated this once and twice and plugged into the DE.

I eventually reached

.

The indicial equation leads to

or

. At first glance it looks like both solutions are acceptable.

So now I get a recurrence relation with

in terms of

and

which isn't what I hoped for. Maybe I shouldn't have proposed a solution of the form

? How would you tackle this problem?

A little 'excamotage' that can make the task more comfortable is to set...

(1)

... and then try to find . First You compute ...

(2)

(3)

... and then insert (1),(2) and (3) in the original DE obtaining...

(4)

... which seems 'easier to be attacked'. Now You search a series solution of (4) imposing proper constrain...

http://www.sv-luka.org/ikone/ikone180a.jpg

Marry Christmas from Serbia

Re: Series solution to DE

Thanks for the tips.

I propose a solution of the form . After derivating once and twice and plugging into the DE the indicial equation gives me or . And also . But the denominator blows up if and , thus only. Thus .

is arbitrary but different from 0. I'm guessing I have to choose it so that as requested by the statement of the problem. When I choose , I get , and in general, . So that . It seems like I got lucky this time when I chose the arbitrary constant since as required.

Is what I've done correct so far?

P.S.:This indeed satisfies the original DE. Wow! Thank you very much!

Post Scriptum 2: I have no idea how to find A(x) and B(x) for a second solution. I tried to write them as infinite series and plug into the DE but I get an enormous expression that doesn't look easily solvable. I'm missing a trick here.

Re: Series solution to DE

Quote:

Originally Posted by

**arbolis** Thanks for the tips.

I propose a solution of the form

. After derivating once and twice and plugging into the DE the indicial equation gives me

or

. And also

. But the denominator blows up if

and

, thus

only. Thus

.

is arbitrary but different from 0. I'm guessing I have to choose it so that

as requested by the statement of the problem. When I choose

, I get

,

and in general,

. So that

. It seems like I got lucky this time when I chose the arbitrary constant

since

as required.

Is what I've done correct so far?

P.S.:This indeed satisfies the original DE. Wow! Thank you very much!

Post Scriptum 2: I have no idea how to find A(x) and B(x) for a second solution. I tried to write them as infinite series and plug into the DE but I get an enormous expression that doesn't look easily solvable. I'm missing a trick here.

If a solution of the second order DE...

(1)

... is known, a second solution independent from it can be found with the following procedure. If u(x) and v(x) are solutions of (1), then is...

(2)

(3)

If You multiply the (2) by v, the (3) by u and do the difference You obtain...

(4)

The (4) is a DE the solution of which is...

(5)

... and from (5), deviding by You obtain...

(6)

Now the solution of (6) is...

(7)

... and that permits You to derive u(x) from v(x)...

(8)

You found as possible solution of (1) and from (8) You derive [finally!]...

(9)

Now all that You have to do is solving the integral in (9)...

http://www.sv-luka.org/ikone/ikone180a.jpg

Marry Christmas from Serbia

Re: Series solution to DE

In the previous post we are arrived to write...

(1)

The next step is of course the computation of the integral in (1). Proceeding with standard integration by parts we obtain...

(2)

... where Ei(*) is the so called Exponential Integral Function, the expansion of which is given by...

(3)

... where is the Euler's constant. It is requested to write y(x) in the form...

(4)

... and then to find the series expansion of a(x) and b(x). From (1) (2) and (3) we derive [neglecting the sign '-']...

(5)

At this point the effective computation of the first coefficients of a(x) and b(x) is not too complex and the task is left to the Reader...

http://www.sv-luka.org/ikone/ikone180a.jpg

Marry Christmas from Serbia