Home

Project Description

Authors

Differential Equations
  with Mathematica

• Preface
• Contents
• Samples
• Ordering

Differential Equations
  with Maple

• Preface
• Contents
• Samples
• Ordering

Differential Equations
  with Matlab

• Preface
• Contents
• Samples
• Ordering
• Scripts
• Glossary
• Series

Differential Equations with MATLAB

Sample Problem

Problem Set D. Second Order Equations

Problem 1.

(a) For t close to zero, the equation resembles y" = 0, which has general solution y = c₁t + c₂. We refer to this as a "facsimile" solution. Graph a numerical solution to Airy's equation with initial conditions y(0)=0, y'(0)=1, and the facsimile solution (with the same initial data) on the interval (-2,2). How well do they match?

(b) For t close to -K^2 << 0, the equation resembles y "= -K^2y, and the corresponding facsimile solution is given by y = c₁ sin(Kt + c₂). Again using the initial conditions y(0)=0, y'(0)=1, plot a numerical solution of Airy's equation over the interval (-18,-14). Using the value K=4, try to find values of c₁ and c₂ so that the facsimile solution matches well with the actual solution. Why shouldn't we expect the initial conditions for Airy's equation to be the appropriate initial conditions for the facsimile solution?

(c) For t close to K^2 >> 0, Airy's equation resembles y" = K^2y, which has solution y = c₁ sinh(Kt + c₂). Plot a numerical solution of Airy's equation together with a facsimile solution (with K=4) on the interval (14,18). In analogy with part (b), you have to choose values for c₁ and c₂ in the facsimile solution. (Try c₁=1, c₂=0.)

(d) Plot the numerical solution of Airy's equation on the interval (-20,2). What does the graph suggest about the frequency and amplitude of oscillations as t --> infinity? Could any of that information have been predicted from the facsimile analysis?