
Differential Equations with MATLABSample ProblemProblem Set D. Second Order EquationsProblem 1.Airy's equation is the linear second order homogeneous equation y" = ty. 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 t close to zero, the equation resembles y" = 0, which has general solution y = c_{1}t + 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 t close to K^2 << 0, the equation resembles y "= K^2y, and the corresponding facsimile solution is given by y = c_{1} sin(Kt + 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 t close to K^2 >> 0, Airy's equation resembles y" = K^2y, which has solution y = c_{1} sinh(Kt + 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 t > infinity? Could any of that information have been predicted from the facsimile analysis? Here is a published MATLAB Mfile containing the solution of this problem. 