University of Florida/Egm4313/s12 Report 3, Problem 3.5

From testwiki
Jump to navigation Jump to search

Problem 3.5

Solved by: Andrea Vargas

Problem Statement

Given y3y+2y=4x26x5

1. Obtain the coefficients of x,x2,x3,x5
2.Verify all equations by long-hand expansion. Use the series before adjusting the indexes.

3. Put the system of equations in an upper triangular matrix.

4. Solve for c0,c1,c2,c3,c4,c5 by using back substitution.

5. Using the initial conditions y(0)=1,y(0)=0 find y(x) and plot it

Solution

1. To obtain the equations for the coefficients we use the following equation from (1) p7-11:

j=03[cj+2(j+2)(j+1)3cj+1(j+1)+2cj]xj3c5(5)x4+2[c4x4+c5x5]=4x26x5

Finding the coefficients of x where j=1:

[c3(3)(2)3(c2)(2)+2c1]x1=[6c36c2+2c1]x

Finding the coefficients of x2 where j=2:

[c4(4)(3)3(c3)(3)+2c2]x2=[12c49c3+2c2]x2

Finding the coefficients of x3 where j=3:

[c5(5)(4)3(c4)(4)+2c3]x3=[20c412c3+2c2]x3

Finding the coefficients of x5 where j=5:

c5x5

2. To verify all the above equations by long hand expansion we use the following equation from (4) p7-12:
j=25cj×j×(j1)×xj23j=15cj×j×xj1+2j=05cj×xj=4x26x5

Finding the coefficients when j=0:

2c0

Finding the coefficients when j=1:

3c1+2c1x1

Finding the coefficients when j=2:

2c26c2x+2c2x2

Finding the coefficients when j=3:

6c3x9c3x2+2c3x3

Finding the coefficients when j=5:

20c5x315c5x4+2c5x5

By collecting these terms we can compare them to the equations of part 1.

Coefficients of x:

[c3(3)(2)3(c2)(2)+2c1]x1=[6c36c2+2c1]x

Coefficients of x2:

[c4(4)(3)3(c3)(3)+2c2]x2=[12c49c3+2c2]x2

Coefficients of x3:

[c5(5)(4)3(c4)(4)+2c3]x3=[20c412c3+2c2]x3

Coefficients of x5:

c5x5

We can see that we obtain the same system of equations to solve for the coefficients with both methods.

3.Constructing the coefficients matrix:
[2320000266000029120000212200000215000002]

Then, the system becomes:

[2320000266000029120000212200000215000002][c0c1c2c3c4c5]=[004006]

4. Solving for the coefficients:
2c5=6c5=3
2c4(15)(3)=0c4=452
2c312(452)+20(3)=0c3=105
2c29(105)+12(452)=4c2=6712
2c16(6712)+6(105)=0c1=13832
2c03(13832)+2(6712)=0c0=28074

Particular solution

This yields the particular solution:

                                                            yp(x)=3x5452x4105x36712x213832x28074

Homogeneous Solution

y'hy'h+2yh=0
λ23λ+2=0
Then, we can find the characteristic equation:
(λ2)(λ1)
λ=2,1
Then the solution for the homogeneous equation becomes:

                                                                                      yh(x)=C1e2x+C2ex

General Solution

Using the given initial conditions y(0)=1,y(0)=0 we find the overall solution:

y(x)=yh+yp
y(x)=C1e2x+C2ex3x5452x4105x36712x213832x28074
y(x)=2C1e2x+C2ex15x490x3315x2671x13832
Using the initial conditions to solve for C1 and C2

1=C1+C228074
0=C1+C213832
C1=454C2=714

The general solution becomes

                                                     y(x)=454e2x+714ex3x5452x4105x36712x213832x28074

Plot

Below is a plot of the solution:

File:R3 5.jpg

--Egm4313.s12.team11.vargas.aa 05:56, 21 February 2012 (UTC)

Template:CourseCat