Problem 3 - Advanced Engineering Mathematics

Question

Construct a table comparing the indicated values of \(y(x)\) using Euler's method, the improved Euler's method, and the RK4 method. Compute to four rounded decimal places. Use \(h=0.1\) and then use \(h=0.05\). $$ y^{\prime}=\sqrt{x+y}, \quad y(0.5)=0.5 $$ \(y(0.6), y(0.7), y(0.8), y(0.9), y(1.0)\)

Step-by-Step Solution

Answer

Compute values using Euler, Improved Euler, and RK4 methods for each \( y(x) \) at specified points with both step sizes, then tabulate the results.

Step 1: Euler's Method Overview

Euler's method updates the current value by using the derivative information. The formula is given by:\[ y_{n+1} = y_n + h \, f(x_n, y_n) \]where \( f(x, y) = \sqrt{x+y} \).We use \( h = 0.1 \) and \( h = 0.05 \).

Step 2: Initial Values

Start with the initial condition: \( y(0.5) = 0.5 \).

Step 3: Euler's Method with h = 0.1

For each step, apply Euler's method:- **At \( x = 0.6 \):** \[ y_{1} = 0.5 + 0.1 \sqrt{0.5 + 0.5} = 0.5 + 0.1 \, \cdot 1 = 0.6 \]- **At \( x = 0.7 \):** \[ y_{2} = 0.6 + 0.1 \sqrt{0.6 + 0.6} = 0.6 + 0.1 \, \cdot 1.0954 = 0.7095 \]- Continue similarly for \( x = 0.8, 0.9, 1.0 \).

Step 4: Improved Euler's Method Overview

Improved Euler (or Heun's method) corrects the overshoot of Euler's method by averaging slopes. Use:\[ y_{n+1} = y_n + \frac{h}{2} ( f(x_n, y_n) + f(x_{n+1}, y_{temp}) ) \]where \( y_{temp} = y_n + h \, f(x_n, y_n) \).

Step 5: Improved Euler's with h = 0.1

Following similar calculations as Euler Method, with correction:- **At \( x = 0.6 \):** - Compute temp step: \( y_{temp} = 0.5 + 0.1 \times 1 = 0.6 \) - Apply correction: \[ y_{1} = 0.5 + \frac{0.1}{2} (1 + \sqrt{0.6 + 0.6}) \approx 0.6048 \]- Continue for \( x = 0.7, 0.8, 0.9, 1.0 \).

Step 6: RK4 Method Overview

RK4 involves four evaluations to calculate \( y_{n+1} \). Calculate using these four stages:- \( k_1 = h f(x_n, y_n) \)- \( k_2 = h f(x_n + \frac{h}{2}, y_n + \frac{k_1}{2}) \)- \( k_3 = h f(x_n + \frac{h}{2}, y_n + \frac{k_2}{2}) \)- \( k_4 = h f(x_n + h, y_n + k_3) \)Then update \( y_{n+1} = y_n + \frac{1}{6}(k_1 + 2k_2 + 2k_3 + k_4) \).

Step 7: RK4 with h = 0.1

Apply RK4 step-by-step:- **At \( x = 0.6 \):** Compute \( k_1, k_2, k_3, k_4 \) and then the update: \[ k_1 = 0.1 \sqrt{0.5 + 0.5}, \; k_2 = 0.1 \sqrt{0.55 + y + k_1/2}, \; k_3 = 0.1 \sqrt{0.55 + y + k_2/2}, \; k_4 = 0.1 \sqrt{0.6 + y + k_3} \]Use these to compute \( y_{1} \), then proceed to the next values similarly.

Step 8: Iterations for h = 0.05

Repeat Steps 3, 5, and 7 but use a step size \( h = 0.05 \). This improves the accuracy of the methods.

Step 9: Constructing the Table

Collect the result for each \( y(x) \) value from each method and list them in a table, showing results for both \( h = 0.1 \) and \( h = 0.05 \) to four decimal places.

Key Concepts

Improved Euler's Method

The Improved Euler's Method, also known as Heun’s method, is a refined version of the basic Euler's method. The main goal of this technique is to reduce the error inherent in the standard Euler's method.

In the Improved Euler's Method, each step is taken by averaging the slopes at the beginning and the end of the interval. This process helps correct the trajectory based on the curvature of the solution.

• First, compute a provisional value of the solution at the next step, known as the ' predictor', using Euler's method.
• Calculate the slope at this predicted point, and then average it with the slope at the beginning of the interval.
• Finally, update your solution using this average slope.

The formula can be summarized as:\[y_{n+1} = y_n + \frac{h}{2} \left( f(x_n, y_n) + f(x_{n+1}, y_{temp}) \right)\]Where:

• \( y_{temp} = y_n + h \, f(x_n, y_n) \) is the temporary solution estimate.

This method offers a marked improvement over basic Euler's technique as it accounts for the changing rate of the solution more accurately.

Runge-Kutta Method (RK4)

The Runge-Kutta Method, particularly the fourth-order version (RK4), is one of the most popular and widely used techniques for solving ordinary differential equations due to its excellent accuracy. RK4 is more precise because it uses several intermediate computations to estimate the change in the solution more accurately over each step.

The process contains four key calculations:

• Compute the initial slope \( k_1 = h \, f(x_n, y_n) \).
• Estimate the slope at the midpoint using \( k_2 = h \, f(x_n + \frac{h}{2}, y_n + \frac{k_1}{2}) \).
• Again, estimate the slope at the midpoint but with a revised function value using \( k_3 = h \, f(x_n + \frac{h}{2}, y_n + \frac{k_2}{2}) \).
• Compute the slope at the end of the interval with \( k_4 = h \, f(x_n + h, y_n + k_3) \).

The value of the function at the next point, \( y_{n+1} \), is then calculated as:\[y_{n+1} = y_n + \frac{1}{6} (k_1 + 2k_2 + 2k_3 + k_4)\]This method leverages multiple estimates (slopes) in various scenarios to improve accuracies, essentially balancing Euler method's speed with better accuracy, making it a robust choice for complex systems.

Numerical Approximation Techniques

Numerical approximation techniques provide a way to solve differential equations when a known analytical solution is not available, or when dealing with complex equations that cannot be easily resolved analytically.

These methods help us approximate the functions by numerically calculating the successive points of a solution curve. Here are some popular numerical methods:

• Euler's Method: A straightforward approach using the derivative, although it can be less accurate.
• Improved Euler's Method: Bridges the gap between simplicity and better accuracy by correcting the slope at each step.
• Runge-Kutta Methods: Offers a balance between computational intensity and accuracy with higher-order methods like RK4 being widely used.

Numerical methods are invaluable in many fields—such as engineering, physics, and finance—where they help predict future behaviors based on current data, aiding in simulations and complex dynamic systems.

Differential Equations

Differential equations form the foundation of many complex phenomena in natural sciences, engineering, and economics. These equations describe the relationship between functions and their derivatives, representing how a certain quantity changes over time or space.

Typical forms include:

• Ordinary Differential Equations (ODEs): Equations involving functions of one variable and their derivatives.
• Partial Differential Equations (PDEs): Involve multiple independent variables and partial derivatives, often used in multivariable systems.

Understanding differential equations gives insight into the processes modeled by them. For instance:

• In physics, they describe dynamics of systems like oscillations in a pendulum or electric circuits.
• In biology, they help model population dynamics, describing changes in species numbers over time.
• In finance, differential equations model the evolution of financial options prices.

Numerical solutions provide practical means to obtain answers to differential equations where analytical solutions are not feasible, enriching our ability to model real-world systems.