The Student Room Group

Multiplying power series

I have a differential equation and I am solving it using power series. I am now at this point where I have got:

n=0(ax)nn!λ=0aλxλ+k1 \displaystyle \sum_{n = 0}^{\infty} \frac{{(-ax)}^n}{n!} \cdot \sum_{\lambda = 0}^{\infty} a_{\lambda}x^{\lambda + k - 1}

The reason I have got it like that was because my original equation read:

eaxxy \displaystyle \frac{e^{-ax}}{x} y and so the minus 1 represents the x on the denominator, the y is obviosuly represented by the power series on the right and I have represented e^(-ax) as a power series as well.

How would I go about simplifying this? So far, I haven't said what n or lambda are so could I simply swap the n's in the exponential power series for lambda and be left with

λ=0((ax)λλ!aλxλ+k1) \displaystyle \sum_{\lambda = 0}^{\infty} \left( \frac{{(-ax)}^{\lambda}}{\lambda!} \cdot a_{\lambda}x^{\lambda + k - 1} \right) and then try and simplify like this?

Although, whilst I'm writing this, I'm thinking if theres a λ! \lambda ! at the bottom, how would I shift the index to help me find the reccurance relation (which I think is the indicial equation?).
Reply 1
It's not true in general that anxnbnxn=anbnxn\sum a_nx^n \sum b_nx^n = \sum a_nb_nx^n, as you seem to be implying. What's the full differential equation you're trying to solve?
Reply 2
As Nuodai remarked - you are confused about how to multiply sums. In technical terms - you are using the fallacy of (a variation on) the Freshman's Dream. Never mind power series, consider the finite sums

n=01anxn=a0+a1x\displaystyle \sum_{n=0}^1 a_n x^n = a_0 + a_1x and m=01bmxn=b0+b1x\displaystyle \sum_{m=0}^1 b_m x^n = b_0 + b_1x

Now

n=01[(anxn)(bnxn)]=a0b0+a1b1x2\displaystyle \sum_{n=0}^1 \left[ (a_n x^n) (b_n x^n)\right] = a_0b_0 + a_1b_1x^2

which clearly isn't equal to

(n=01anxn)(m=01bmxn)=(a0+a1x)(b0+b1x)=a0b0+a0b1x+a1b0x+a1b1x2\displaystyle \left(\sum_{n=0}^1 a_n x^n\right) \cdot\left(\sum_{m=0}^1 b_m x^n\right) = (a_0 + a_1x)(b_0 + b_1x) = a_0b_0 + a_0b_1x + a_1b_0x + a_1b_1x^2

...Anyway once you have dealt with that; as to (at least how I interpret what you are asking in your) final question:

Firstly reindex your power series expansion for yy by writing
m=0amxm+k1=m=0bmxm\displaystyle \sum_{m=0}^\infty a_m x^{m+k-1}=\sum_{m=0}^\infty b_m x^{m} where bm=am(k1)b_m = a_{m-(k-1)}

Then (n=0cnxn)(m=0bmxm)=n=0(i=0ncibni)xn\displaystyle \left(\sum_{n=0}^\infty c_n x^n \right)\cdot \left(\sum_{m=0}^\infty b_m x^m \right) = \sum_{n=0}^\infty \left(\sum_{i=0}^n c_i b_{n-i}\right) x^n

If you prefer to express this in terms of the original power series coefficients for yy then write:

n=0(i=0ncibni)xn=n=0(i=0nciani(k1))xn \displaystyle \sum_{n=0}^\infty \left(\sum_{i=0}^n c_i b_{n-i}\right) x^n = \sum_{n=0}^\infty \left(\sum_{i=0}^n c_i a_{n-i-(k-1)}\right) x^n
(edited 11 years ago)
Original post by nuodai
It's not true in general that anxnbnxn=anbnxn\sum a_nx^n \sum b_nx^n = \sum a_nb_nx^n, as you seem to be implying. What's the full differential equation you're trying to solve?


The equation is:

h22md2ydx2+(EAeaxx)y=0 \frac{h^2}{2m} \frac{d^2y}{dx^2} + (E - \frac{Ae^{-ax}}{x}) y = 0

Where h, m, A < 0 and a and E are constants.

Because its a Frobenius method question I used
Unparseable latex formula:

y = \sum^{\nifty}_{\lambda} = a_k x^{k + \lambda}

. I use 'x' rather than 'x - x_0' as x_0 = 0 because thats the singular point (I'm just writing this out to make sure I understand the reason why I've done it).

So when I worked it out, I out in that power series and got:

Unparseable latex formula:

E \cdot y - Ae^{-ax} \sum^{\nifty}_{\lambda} a_k x^{k + \lambda - 1}



Now I'm trying to get rid of the e^x and so I thought about using the maclaurin series of e^x and then simplifying.
Reply 4
Original post by claret_n_blue

Now I'm trying to get rid of the e^x and so I thought about using the maclaurin series of e^x and then simplifying.


Do you now understand how to reindex and multiply power series?
Reply 5
Original post by claret_n_blue
I don't get what you've done here. How did you get that on the rhs?


As in my first post - you are basically forgetting to cross multiply and add indices. If you multiply two sums you have to multiply each term of the first sum with each term of the second sum. In your formula - you just multiply the first term with the first term, the second with the second etc. and you also forget that anxnbnxn=anbnx2nanbnxna_nx^n \cdot b_nx^n = a_nb_n x^{2n} \neq a_nb_n x^{n}


Look at the example of a power series where all of the coefficients except for the first two terms are zero:

n=0anxn,m=0bmxm\displaystyle \sum_{n=0}^\infty a_nx^n, \sum_{m=0}^\infty b_mx^m

where ai=0=bia_i = 0 = b_i for all i2i \geq 2.

Thus both series are just finite sums of two terms.

Lets just expand:
m=0n=0manbmnxm\displaystyle \sum_{m=0}^\infty \sum_{n=0}^m a_n b_{m-n} x^m
out term by term. We get

a0b0+(a0b1+a1b0)x1+a1b1x2+a_0b_0 + (a_0b_1 +a_1b_0)x^1 + a_1b_1 x^2 + \cdots

where the \cdots is all zero since the terms a2,b2,a3,b3,a_2,b_2,a_3,b_3,\ldots are all zero by definition.

For eachxmx^m, we sum over all combinations of powers of xx from each term that multiply to give xmx^m. So for example, for x1x^1, we can multiply the constant term from the first sum and the xx term from the second and also the xx term from the first and the cconstant term from the second. We do things this way to collect like powers - nothing more and nothing less.
(edited 11 years ago)
Original post by Jake22
Do you now understand how to reindex and multiply power series?


Ok, so I understand it all now, so that leaves my differential equation to read:

y=n=0anxn+k,C=h22m\displaystyle y = \sum_{n = 0}^{\infty} a_n x^{n + k}, C = \frac{h^2}{2m}

so

Cn(n+k)(n+k1)xn+k1+En=0anxn+kn=0(i=0n(ax)ii!ani(k1)) \displaystyle C \cdot \sum_{n}^{\infty} (n + k)(n + k - 1) x ^{n + k - 1} + E \cdot \sum_{n = 0}^{\infty} a_n x^{n + k} - \sum_{n = 0}^{\infty} \left( \sum_{i = 0}^{n} \frac{(-ax)^i}{i!} \cdot a_{n - i - (k - 1)} \right)

Is that correct?
Reply 7
Forgetting about the differential equation for a moment* - it is not true that

(n=0(ax)nn!)(m=0amxm+k1)=n=0((ax)ii!ani(k1))\displaystyle \left(\sum_{n=0}^\infty \frac{(-ax)^n}{n!}\right) \cdot \left(\sum_{m=0}^\infty a_m x^{m+k-1}\right) = \sum_{n=0}^\infty \left(\frac{(-ax)^i}{i!}\cdot a_{n-i-(k-1)}\right)

Try to understand the summation notation for power series. Until you take the time to understand it you are just plugging things randomly into formulae you don't understand and are making mistakes. In particular - seperate your coefficients from your powers of xx.

So, rewrite n=0(ax)nn!\displaystyle \sum_{n=0}^\infty \frac{(-ax)^n}{n!} as n=0(a)nn!xn\displaystyle \sum_{n=0}^\infty \frac{(-a)^n}{n!}x^n and try again...


*You have to forget about the differential equation for the time being. You are making the classical undergraduate mistake of trying to solve a problem x by using y without understanding what y is or how to manipulate it. In other words get comfortable with basic manipulation of power series before you use them to solve problems. In particular, be able to add, multiply, divide and reindex them.
Original post by Jake22
You have to forget about the differential equation for the time being. You are making the classical undergraduate mistake of trying to solve a problem x by using y without understanding what y is or how to manipulate it. In other words get comfortable with basic manipulation of power series before you use them to solve problems. In particular, be able to add, multiply, divide and reindex them.


The problem I'm having is in the letter. After what you wrote in that post, what I have done is this:

Unparseable latex formula:

\displaystyle y = \sum_{\m = 0}^{\infty} a_m x^{m + k - 1} \, , \, e^{-ax} = \sum_{n = 0}^{\infty} \frac{(-a)^n}{n!} \cdot x^n



So eaxy= e^{-ax} \cdot y =

Unparseable latex formula:

\displaystyle \sum_{n = 0}^{\infty} \frac{(-a)^n}{n!} \cdot x^n \cdoy\sum_{m = 0}^{\infty} a_m x^{m + k - 1}



Now the bit I am stuck on is the actual multiplication. I understand how every term in the first sum must be multiplied by every term in the second term but I don't get how to write that?

I'm trying to use your very first point to help me write it like this but I don't get what 'i' represents and how you got it.
Reply 9
Original post by claret_n_blue


Now the bit I am stuck on is the actual multiplication. I understand how every term in the first sum must be multiplied by every term in the second term but I don't get how to write that?

I'm trying to use your very first point to help me write it like this but I don't get what 'i' represents and how you got it.


Reread my posts #3 and #6 - you need to try and see how it works with finite sums before you see how it works in general with power series. I wrote out an example there. The point about the formula I mentioned is that we collect like terms to present the final answer in the standard form for a power series i.e. we collect together all the terms involving x, collect all the terms involving x^2,... etc.

i is just an index - it doesn't matter which letter you take.

Again, let n=0anxn\displaystyle \sum_{n=0}^\infty a_n x^n and m=0bmxm\displaystyle \sum_{m=0}^\infty b_m x^m be power series. Then if we multiply both together, we do the following:

multiply the first term of the first series by everything in the second:

a0x0m=0bmxm=m=0a0bmx0xm=m=0a0bmxm\displaystyle a_0x^0 \sum_{m=0}^\infty b_m x^m = \sum_{m=0}^\infty a_0b_m x^0x^m=\sum_{m=0}^\infty a_0b_m x^m

then multiply the second term of the first series by everything in the second:

a1x1m=0bmxm=m=0a1bmx1xm=m=0a1bmxm+1\displaystyle a_1x^1 \sum_{m=0}^\infty b_m x^m = \sum_{m=0}^\infty a_1b_m x^1x^m = \sum_{m=0}^\infty a_1b_m x^{m+1}

and so on and on. We then have to add all of these expressions together.

More generally, when we multiply the nnth term of the first series by everything in the second - we get

anxnm=0bmxm=m=0anbmxnxm=m=0anbmxm+n\displaystyle a_nx^n \sum_{m=0}^\infty b_m x^m = \sum_{m=0}^\infty a_nb_m x^nx^m = \sum_{m=0}^\infty a_nb_m x^{m+n}

so the whole sum is therefore

m=0a0bmxm+m=0a1bmxm+1++m=0anbmxm+n+\displaystyle \sum_{m=0}^\infty a_0b_m x^{m} + \sum_{m=0}^\infty a_1b_m x^{m+1} + \cdots + \sum_{m=0}^\infty a_nb_m x^{m+n} + \cdots

in summation notation, we could write this as:

n=0(m=0anbmxm+n)\displaystyle \sum_{n=0}^\infty \left( \sum_{m=0}^\infty a_nb_m x^{m+n} \right)

but this isn't that useful to us since we want to present the answer in the standard form of a power series. Instead of adding all those series like that - we collect like terms:

Question: Which terms in the previous sum are all just some coefficient times, say, x3x^3?

Answer: The general term is anbmxm+na_nb_m x^{m+n}. This is 'coefficient times x3x^3' precisely when m+n=3m+n = 3. Thus the relevant terms are a0b3x0+3,a1b2x1+2,a2b1x2+1,a3b0x3+0a_0b_3 x^{0+3}, a_1b_2 x^{1+2}, a_2b_1 x^{2+1},a_3b_0 x^{3+0}. If we write down the sum of all those terms, we just have:

m=03amb3mxm+(3m)=m=03amb3mx3\displaystyle \sum_{m=0}^3 a_mb_{3-m}x^{m + (3-m)} = \sum_{m=0}^3 a_mb_{3-m}x^{3}

Now, similarly, all the terms of the form 'coefficient times xnx^n' are just the sum:

m=0nambnmxn\displaystyle \sum_{m=0}^n a_mb_{n-m}x^{n}

Obviously, the total product is the sum of all the terms 'coefficient times x^0', 'coefficient times x^1', 'coefficient times x^2' etc. i.e.

m=00amb0mx0+m=01amb1mx1++m=0nambnmxn+\displaystyle \sum_{m=0}^0 a_mb_{0-m}x^{0} + \sum_{m=0}^1 a_mb_{1-m}x^{1} + \cdots + \sum_{m=0}^n a_mb_{n-m}x^{n} + \cdots

or in summation notation:

n=0(m=0nambnmxn)\displaystyle \sum_{n=0}^\infty \left(\sum_{m=0}^n a_mb_{n-m}x^{n}\right)

I cannot emphasise enough: your problem is not understanding the summation notation. You must do examples with finite sums where you write out the same thing in series notation and in conventional sum notation if you want to understand how to manipulate infinite sums. e.g.

n=03anxn=a0+a1x+a2x2+a3x3\displaystyle \sum_{n=0}^3 a_n x^n = a_0 + a_1x + a_2 x^2 + a_3 x^3

Now multiply this by the series m=03bmxm=b0+b1x+b2x2+b3x3\displaystyle \sum_{m=0}^3 b_m x^m = b_0 + b_1x + b_2 x^2 + b_3 x^3

presenting your answer as a sum of powers of xx where you collect together like terms writing it out both ways.
(edited 11 years ago)
Original post by ttoby
.
Original post by Jake22
.


Hello. So I did "everything" and after I multiplied these power series, I showed it to my tutor and he told me a simpler way would be to expand the e^x power series and multiplt term by term. So I did all that and after shifting the index I ended up with:

(h22m(k2k)b0)xk2+(h22m(k2+k)b0Ab1)xk1+λ=2[h22m(λ+k1)(λ+k)bλAbλ1+(EAa)bλ2]xλ+k2=0\left( \frac{h^2}{2m}(k^2 - k)b_0\right)x^{k - 2} +\left(\frac{h^2}{2m}(k^2 +k) b_0 - Ab_1 \right)x^{k -1} + \sum_{\lambda = 2}^{\infty} \left[ \frac{h^2}{2m}(\lambda + k - 1)(\lambda + k)b_{\lambda} -A b_{\lambda - 1} + (E - Aa)b_{\lambda - 2} \right] x^{\lambda +k - 2} = 0

From here, I said my indicial equation is (k)(k-1) = 0 and I need to use k = 1 to solve the equation because the question says "By direct substitution, use the larger root of the indicial equation for ...."

So I subbed in k = 1 into my reccurance relation and ended up with

bλ2bλ1=AC(λ2+λ)EAabλ[br]\frac{b_{\lambda - 2}}{b_{\lambda - 1}} = \frac{A - C (\lambda^2 + \lambda)}{E - Aa} b_{\lambda}[br]

Where A, C=h22mC = \frac{h^2}{2m} , E and a are constants.

Do you think this is correct? I'm not sure how to go from here and find my solution. I have a solution that I need to work to so I can show that if it helps.

EDIT: I just re-tried it and it's clearly wrong because I've rearranged it to get my reccurance relation wrong. What do I do from:

λ=2[h22m(λ+k1)(λ+k)bλAbλ1+(EAa)bλ2]xλ+k2=0 \displaystyle \sum_{\lambda = 2}^{\infty} \left[ \frac{h^2}{2m}(\lambda + k - 1)(\lambda + k)b_{\lambda} -A b_{\lambda - 1} + (E - Aa)b_{\lambda - 2} \right] x^{\lambda +k - 2} = 0

I know you do:

h22m(λ+k1)(λ+k)bλAbλ1+(EAa)bλ2=0 \frac{h^2}{2m}(\lambda + k - 1)(\lambda + k)b_{\lambda} -A b_{\lambda - 1} + (E - Aa)b_{\lambda - 2} = 0 but I don't get how to reaarange from here.
(edited 11 years ago)

Quick Reply

Latest