Differential Equations with Mathematica, Second Edition
Sample Problem
Problem Set D. Second Order Equations
Problem 1.
Airy's equation is the linear second order homogeneous equation y"=xy. Although it arises in a number of applications, including quantum mechanics, optics and waves, it cannot be solved exactly by the standard symbolic methods. In order to analyze the solution curves, let us reason as follows.(a) For x close to zero, the equation resembles y"=0, which has general solution y=c_1x+c_2. 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 x close to -K^2 << 0, the equation resembles y"=-K^2y, and the corresponding facsimile solution is given by y=c_1 sin(Kx+ c_2). 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_1 and c_2 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 x close to K^2 >> 0, Airy's equation resembles y"=K^2y, which has solution y=c_1 sinh(Kx+c_2). 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_1 and c_2 in the facsimile solution. (Try c_1=1, c_2=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 x->-infinity? Could any of that information have been predicted from the facsimile analysis?
Here is a Postscript image of a Mathematica Notebook solution. If your browser is configured to view Mathematica Notebooks, you can view the actual Mathematica Notebook solution.