Problem 22 - Advanced Engineering Mathematics

Question

(a) Experiment with acalculator to find an interval \(0 \leq \theta<\theta_{1}\), where \(\theta\) is measured in radians, for which you think \(\sin \theta \approx \theta\) is a fairly good estimate. Then use a graphing utility to plot the graphs of \(y=x\) and \(y=\sin x\) on the same coordinate axes for \(0 \leq x \leq \pi / 2\). Do the graphs confirm your observations with the calculator? (b) Use a numerical solver to plot the solutions curves of the initial-value problems $$ \begin{aligned} \quad \frac{d^{2} \theta}{d t^{2}}+\sin \theta=0, & \theta(0)=\theta_{0}, \theta^{\prime}(0)=0 \\ \text { and } \quad \frac{d^{2} \theta}{d t^{2}}+\theta=0, \quad \theta(0)=\theta_{0}, \theta^{\prime}(0)=0 \end{aligned} $$ for several values of \(\theta_{0}\) in the interval \(0 \leq \theta<\theta_{1}\) found in part (a). Then plot solution curves of the initialvalue problems for several values of \(\theta_{0}\) for which \(\theta_{0}>\theta_{1}\)

Step-by-Step Solution

Answer

For \( 0 \leq \theta < 0.2 \), \( \sin \theta \approx \theta \) is a good approximation. Graphs confirm this observation.

Step 1: Experiment with the Calculator

Use a calculator to evaluate \( \sin \theta \) and compare it to \( \theta \) for various values of \( \theta \) within the range \( 0 \leq \theta < \pi/2 \). Start with small increments, like \( 0.1 \) or \( 0.01 \). Notice at which point \( \sin \theta \) deviates significantly from \( \theta \). You will find that \( \sin \theta \approx \theta \) holds well for \( \theta < 0.2 \) radians.

Step 2: Graphing Functions

Use a graphing utility to plot the graphs of \( y = x \) and \( y = \sin x \) in the interval \( 0 \leq x \leq \pi/2 \). Observe the two graphs; they should be close to each other for small values of \( x \). The graphs confirm that \( \sin \theta \approx \theta \) is a good approximation when \( 0 \leq \theta < 0.2 \) radians.

Step 3: Solve and Plot ODEs for \( \theta_{0} < \theta_{1} \)

Use a numerical solver to solve the differential equations \( \frac{d^{2} \theta}{d t^{2}}+\sin \theta=0 \) and \( \frac{d^{2} \theta}{d t^{2}}+\theta=0 \) with initial conditions \( \theta(0)=\theta_{0}, \theta'(0)=0 \). Choose several values of \( \theta_{0} \) within \( 0 \leq \theta_{0} < 0.2 \). Plot the solutions and observe their behavior.

Step 4: Solve and Plot ODEs for \( \theta_{0} > \theta_{1} \)

Similarly, solve the same differential equations but choose initial values \( \theta_{0} > 0.2 \). Plot these solutions and compare how they differ from the ones where initial values were \( \theta_{0} < 0.2 \). The solutions for \( \theta_{0} > 0.2 \) will likely demonstrate larger oscillations or deviations from linear behavior.

Key Concepts

Approximation Techniques

Approximation techniques are incredibly useful when dealing with functions that can be complex to evaluate directly. In the context of differential equations, finding where \( \sin \theta \approx \theta \) simplifies our calculations significantly when \( \theta \) is small. We'll often find that simple linear approximations, such as this, work sufficiently well in the real world to provide a sensible answer for small ranges of \( \theta \). This is because the Taylor series expansion of the sine function around zero starts with \( \sin \theta = \theta - \frac{\theta^3}{6} + \cdots \), where the first term is a linear function of \( \theta \).

• Use for small angles: Linear approximations for sine are useful for angles less than approximately 0.2 radians, as the higher-order terms become negligible, simplifying the calculation.


• Saves computational resources: By using approximations, we can solve equations faster with less computational overhead.


• Effective domain: For the sine function's linear approximation, the domain is typically for small angles, providing a precise estimate in this limited scope.

Graphical Analysis

Graphical analysis involves plotting functions to visually analyze their behavior, particularly when theoretical analysis becomes cumbersome. In our exercise, plotting \( y=x \) alongside \( y=\sin x \) reveals how close these functions are to each other in their behavior when \( x \) is small.

• Visual verification: By graphing \( y=x \) and \( y=\sin x \), one can clearly see how accurately \( \sin \theta \approx \theta \) holds for \( \theta < 0.2 \).


• Identifying deviation: Graphs allow us to quickly identify where the approximation begins to break down and deviations occur, aiding in setting boundaries for feasible approximations.


• Understanding dynamics: Graphical inspection provides insight into how different functions represent phenomena dynamically over a range of interest.

Thus, plotting graphically confirms hypothesis and provides a clear validation method for approximations.

Numerical Solvers

Numerical solvers are critical in solving differential equations that are hard to solve analytically. These solvers use algorithms to approximate solutions through iterative processes. For the exercise, numerical solvers help visualize the solutions of the second-order differential equations by generating solution curves for initial-value problems.

• Iterative Solution: Numerical solvers break down complex differential equations into smaller, manageable steps, iterating towards an approximate solution.


• Versatility: They handle a broad range of problems, including initial-value problems, where initial conditions are given.


• Exploration: Allowing exploration of how solution behaviors change with different parameters, giving insight where theoretical solutions aren't feasible.


• Tools: Software like MATLAB, Mathematica, or Python’s scipy library are commonly used tools to implement numerical solvers.

Applying these solvers provides valuable insights into dynamics that would otherwise be hidden in unsolvable equations.

Initial Value Problems

Initial value problems (IVPs) are a crucial part of studying differential equations. These are problems where solutions are calculated based on specified initial conditions. In this exercise, the focus is on solving second-order linear differential equations using given initial values.

• Specificity: IVPs are tailored for particular scenarios by setting initial conditions that align with what's being modeled.


• Practicality: They mirror real-world problems where systems begin at known states and progress over time.


• Integration with solvers: Numerical solvers require initial conditions to iterate and find approximate solutions.


• Sensitivity: IVPs help analyze the sensitivity of solutions to initial conditions, revealing how small changes can affect overall behavior.

By understanding initial value problems, one can effectively model and predict outcomes in dynamic systems, making this approach invaluable in applied mathematics and physics.