Problem 39 - Advanced Engineering Mathematics

Question

(a) Find two power series solutions for \(y^{\prime \prime}+x y^{\prime}+y=0\) and express the solutions \(y_{1}(x)\) and \(y_{2}(x)\) in terms of summation notation. (b) Use a CAS to graph the partial sums \(S_{N}(x)\) for \(y_{1}(x)\). Use \(N=2,3,5,6,8,10 .\) Repeat using the partial sums \(S_{N}(x)\) for \(y_{2}(x)\). (c) Compare the graphs obtained in part (b) with the curve obtained using a numerical solver. Use the initial conditions \(y_{1}(0)=1, y_{1}^{\prime}(0)=0\), and \(y_{2}(0)=0, y_{2}^{\prime}(0)=1\). (d) Rexamine the solution \(y_{1}(x)\) in part (a). Express this series as an elementary function. Then use (5) of Section \(3.2\) to find a second solution of the equation. Verify that this second solution is the same as the power series solution \(y_{2}(x)\).

Step-by-Step Solution

Answer

Find power series solutions, graph partial sums, compare with numerical solver, express as elementary function.

Step 1: Setup the power series solution

Assume a power series solution of the form:\[ y(x) = \sum_{n=0}^{\infty} a_n x^n \]Compute its derivatives:\[ y'(x) = \sum_{n=1}^{\infty} n a_n x^{n-1} \]\[ y''(x) = \sum_{n=2}^{\infty} n(n-1) a_n x^{n-2} \]

Step 2: Substitute into the ODE

Substitute the series and its derivatives into the given differential equation:\[ \sum_{n=2}^{\infty} n(n-1) a_n x^{n-2} + x \sum_{n=1}^{\infty} n a_n x^{n-1} + \sum_{n=0}^{\infty} a_n x^n = 0 \].Re-index and combine terms to match powers of \(x\).

Step 3: Equate coefficients

Re-index and shift the summation indices to express terms with the same powers of \(x\). Equate the coefficients of like powers of \(x\) to zero to find recurrence relations for \(a_n\).

Step 4: Solve the recurrence relations

From the recurrence relations, you will find two linearly independent sets of coefficients \(a_n\), corresponding to \(y_1(x)\) and \(y_2(x)\).

Step 5: Express solutions as summations

Express the solutions \(y_1(x)\) and \(y_2(x)\) as power series in summation notation:\[ y_1(x) = \sum_{n=0}^{\infty} a_n x^n \]\[ y_2(x) = \sum_{n=0}^{\infty} b_n x^n \], where \(a_n\) and \(b_n\) are derived from Step 4.

Step 6: Graph partial sums with CAS

Using a Computer Algebra System (CAS), compute and graph the partial sums \(S_N(x)\) for different values of \(N\) for both \(y_1(x)\) and \(y_2(x)\).

Step 7: Compare with numerical solver

Use a numerical solver to compute a solution with the initial conditions specified, and compare it graphically with the graphs of the partial sums generated in Step 6.

Step 8: Examine elementary function form

Identify if \(y_1(x)\) can be expressed as a familiar elementary function.

Step 9: Find a second solution using reduction of order

Use the method of reduction of order on the elementary function form found in Step 8 to find a second solution, and verify it matches the power series solution for \(y_2(x)\).

Key Concepts

Ordinary Differential Equations

Ordinary Differential Equations (ODEs) are equations involving a function of one independent variable and its derivatives. In mathematics, they help us understand how certain quantities change over time or space. For example, the equation \(y'' + xy' + y = 0\) is a second-order linear ODE with variable coefficients. Solving such equations can be challenging, but it provides insights into the dynamics of systems, such as mechanical motion or electrical circuits.

In this particular exercise, we use the method of power series to find solutions to the ODE. This involves assuming solutions in the form of an infinite series and finding coefficients that satisfy the equation. Although ODEs can sometimes be solved analytically, many real-world problems require numerical approaches or specialized techniques like power series to yield meaningful solutions.

CAS (Computer Algebra Systems)

Computer Algebra Systems (CAS) are powerful tools that assist in solving complex mathematical problems. They are especially valuable for working with ODEs, where symbolic manipulation and simplification can become overwhelmingly complicated. In the context of this exercise, CAS software is used to graph the partial sums of power series solutions. This graphical representation helps visualize the behavior of solutions over a range of values.

Using a CAS allows us to explore solutions more deeply than manual calculation would permit. With CAS, we can quickly compute large sums, perform algebraic operations, and even verify the accuracy of our solutions. The flexibility and computational power of a CAS make it a preferred choice for students and professionals dealing with intricate differential equations.

Numerical Solver

A numerical solver is a computational tool designed to find approximate solutions to mathematical problems that may not be easily solvable analytically. For ODEs, numerical solvers employ algorithms to approximate the behavior of unknown functions over a specified domain.

In practice, numerical solvers produce results that offer a way to validate or contrast with solutions obtained through other means, such as power series. In this exercise, we compare graphs of partial sums, constructed with a CAS, against the output of a numerical solver, ensuring that the initial conditions are met. This comparison can highlight the strengths and limitations of each method, providing a comprehensive understanding of solution behavior.

Reduction of Order

Reduction of Order is a technique used in solving linear differential equations, particularly when a known solution exists and a second, independent solution is required. It serves as a method to construct another linearly independent solution from one that is already known.

By assuming the second solution has a particular structure dependent on the first, reduction of order simplifies the process of solving ODEs. In the given exercise, after expressing one solution as an elementary function, this method aids in directly deriving the second solution without starting from scratch. This is especially helpful for verifying results obtained via power series, ensuring that solutions satisfy the original differential equation.