Let us analyze Airy's equation y''= xy. This equation arises in many different contexts, most commonly as a limiting case of Schrodinger's equation.
For x near zero, the approximate equation is y''=0, and this has the general solution y=
![[Graphics:Images/AiryEGs_gr_2.gif]](Images/AiryEGs_gr_2.gif)
![[Graphics:Images/AiryEGs_gr_4.gif]](Images/AiryEGs_gr_4.gif)
![[Graphics:Images/AiryEGs_gr_5.gif]](Images/AiryEGs_gr_5.gif)
The solution agrees with the simple approximation on a decently large interval about zero. Eyeballing it, we can't tell the difference between them on [-.2, .4].
![[Graphics:Images/AiryEGs_gr_6.gif]](Images/AiryEGs_gr_6.gif)
![[Graphics:Images/AiryEGs_gr_7.gif]](Images/AiryEGs_gr_7.gif)
![[Graphics:Images/AiryEGs_gr_9.gif]](Images/AiryEGs_gr_9.gif)
![[Graphics:Images/AiryEGs_gr_10.gif]](Images/AiryEGs_gr_10.gif)
Thus is "good enough" means 5%, then this approximation is good enough on the interval (-.8,.8)
Now let's turn to x very large negative. If x ~ -
, then we can try to approximate with y —
. We'll try ω=5. For starters, set c1=1 (no scaling), c2=0 (no phaseshift)
![[Graphics:Images/AiryEGs_gr_13.gif]](Images/AiryEGs_gr_13.gif)
![[Graphics:Images/AiryEGs_gr_14.gif]](Images/AiryEGs_gr_14.gif)
We see that the Airy function does have approximately the same period, and that we can scale (c1=.55) and phase shift (c2=-2.7)
![[Graphics:Images/AiryEGs_gr_15.gif]](Images/AiryEGs_gr_15.gif)
![[Graphics:Images/AiryEGs_gr_16.gif]](Images/AiryEGs_gr_16.gif)
Pretty neat huh? What is the error?
![[Graphics:Images/AiryEGs_gr_17.gif]](Images/AiryEGs_gr_17.gif)
![[Graphics:Images/AiryEGs_gr_18.gif]](Images/AiryEGs_gr_18.gif)
Let's see if this solution is good on bigger intervals:
![[Graphics:Images/AiryEGs_gr_19.gif]](Images/AiryEGs_gr_19.gif)
![[Graphics:Images/AiryEGs_gr_20.gif]](Images/AiryEGs_gr_20.gif)
Now it is starting to get away from us, mainly because the period is changing. What else changes towards the end points of the intervals?
Do you think that for higher frequencies you will get a better or worse approximation (meaning: valid on a longer or shorter interval)? Guess, and experiment.
Now let's go to the other extremem, x ~
, then we can try to approximate with y —
. Let's use again k=5.
![[Graphics:Images/AiryEGs_gr_23.gif]](Images/AiryEGs_gr_23.gif)
![[Graphics:Images/AiryEGs_gr_24.gif]](Images/AiryEGs_gr_24.gif)
Well, our approximation is off by many orders of magnitude, so let's start fiddling with the knobs......
![[Graphics:Images/AiryEGs_gr_25.gif]](Images/AiryEGs_gr_25.gif)
![[Graphics:Images/AiryEGs_gr_26.gif]](Images/AiryEGs_gr_26.gif)
There the shapes look to be about the same. Now let's twiddle the other dial....
![[Graphics:Images/AiryEGs_gr_27.gif]](Images/AiryEGs_gr_27.gif)
![[Graphics:Images/AiryEGs_gr_28.gif]](Images/AiryEGs_gr_28.gif)
Not bad.
Let's solve the equation for a variety of initial conditions.
![[Graphics:Images/AiryEGs_gr_29.gif]](Images/AiryEGs_gr_29.gif)
![[Graphics:Images/AiryEGs_gr_30.gif]](Images/AiryEGs_gr_30.gif)
![[Graphics:Images/AiryEGs_gr_32.gif]](Images/AiryEGs_gr_32.gif)
![[Graphics:Images/AiryEGs_gr_33.gif]](Images/AiryEGs_gr_33.gif)
![[Graphics:Images/AiryEGs_gr_34.gif]](Images/AiryEGs_gr_34.gif)
![[Graphics:Images/AiryEGs_gr_35.gif]](Images/AiryEGs_gr_35.gif)
![[Graphics:Images/AiryEGs_gr_36.gif]](Images/AiryEGs_gr_36.gif)