Series solution for ode by undetermined coefficients

Click For Summary

Homework Help Overview

The discussion revolves around finding a Taylor series solution up to order 3 for a nonlinear ordinary differential equation (ODE) given by y' = x^2 + sin(y) with the initial condition y(0) = π/2.

Discussion Character

  • Exploratory, Assumption checking, Problem interpretation

Approaches and Questions Raised

  • Participants discuss the form of the series solution and the implications of substituting the power series for sin(y). There are attempts to equate coefficients of like powers of x to derive relationships between coefficients a_0, a_1, a_2, and a_3. Some participants question whether the initial condition was applied correctly and suggest alternative approaches to determine the coefficients.

Discussion Status

The discussion is ongoing, with participants exploring different interpretations of the problem and questioning the assumptions made in the original poster's approach. Some guidance has been offered regarding the use of the initial condition to simplify the solution process.

Contextual Notes

There is a noted concern about the expansion of sin(y) around y=0 instead of y=π/2, which may affect the results. Additionally, participants are considering the implications of missing terms in the series expansion.

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
 

Similar threads

  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 16 ·
Replies
16
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 5 ·
Replies
5
Views
1K
Replies
8
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
12
Views
2K