Series solution for ode by undetermined coefficients

John 123
Messages
32
Reaction score
0

Homework Statement


Obtain the Taylor series solution up to and including order 3 of the following non linear ode
<br /> y&#039;=x^2+\sin y,y(0)=\frac{\pi}{2}<br />

Homework Equations


After substituting the power series form of sin(y) I get:
<br /> y&#039;=x^2+(y-\frac{y^3}{3!}+\frac{y^5}{5!}-\frac{y^7}{7!}...)<br />


The Attempt at a Solution


We require a series solution of the form:
<br /> y=a_0+a_1x+a_2x^2+a_3x^3+a_4x^4+...<br />
Then
<br /> y&#039;=a_1+2a_2x+3a_3x^2+4a_4x^3+...<br />
Substituting y and y' in the ode and equating coefficients of like powers of x gives:
<br /> a_1=a_0-\frac{a_0^3}{6}<br />
<br /> 2a_2=a_1-\frac{a_0^2a_1}{2}<br />
<br /> 3a_3=1+a_2-\frac{a_0a_1^2}{2}-\frac{a_2a_0^2}{2}<br />
Then expressing all the a's in terms of a_0 and using :
<br /> a_0=y(0)=\frac{\pi}{2}<br />
I get a result completely different to the book answer and cannot see my error.
BOOK ANSWER IS:
<br /> y=\frac{\pi}{2}+x+\frac{1}{6}x^2+<br />





 
Physics news on Phys.org
Could it be that a_1 = sin a_0 instead and you just forgot infinite minus two terms from the expression? I don't get the book answer though, quickly glancing it would seem to me that a_2 = 0.
 
Perhaps it should be \pi / 2 + x + x^3/6 ?
 
One more thing: it might be pretty hard to solve the way you are trying to. Instead you can be a bit tricky. You are given the initial condition for y(0). This allows you to determine a_0=y_0 instantly. Now it's easy to solve y'(0) from the DE, giving you again instantly a_1. Then you can simply differentiate the DE again, solving y''(0) and y'''(0) as well.
 
Hi Clamtrox
I believe my error is in expanding siny about y=0 instead of about y=pi/2.
Regards
John
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top