next up previous contents
Next: Example Up: First Order ODEs Symbolically Previous: Example 2   Contents

Using hw3plot.m

Tying graphical, symbolic, and numerical representations together is an important goal of this course. hw3plot.m is designed to bring these together on one page.
  1. Solve the indicated ODE by hand using whichever method is appropriate to the textbook section containing the problem. This is to be done on scratch paper.
  2. Make the appropriate entries in init.m for ftx and ftxsln. Make an educated guess for initax-values. Be sure some value is entered for each of t0, x0, C. Save then type init in the command window to initialize the variables (or copy-and-paste in the Emporium).
  3. Type slopef(ftx,N) to check that the plot window size is appropriate. Adjust initax in init.m if necessary. Save then type init in the command window (or copy-and-paste in the Emporium).
  4. Type sol1 and choose initial conditions (i.e., (t0,x0)) that produce a solution illustrating the ``interesting'' behaviour of the ODE. You may need to experiment here.
  5. Note the t0 and x0-values returned by sol1 at the bottom of the plot and enter these values in init.m.
  6. Solve for C by substituting t0 and x0 into the general symbolic ODE solution you obtained by hand, then enter the result into init.m.
  7. Save init.m then type init in the command window to initialize the new variables (or copy-and-paste in the Emporium).
  8. Type hw3plot. A ``y'' answer to ``Do you already have a separate figure window for hw3plot? (y/n)'' creates a new MATLAB Figure window specifically for hw3plot. If you have already used hw3plot in this session, an appropriate window should already exist, and you can answer ``n ``.
    In your Figure window, you should have three plots.
    1. Upper right contains Figure 1, the slopefield for ftx.
    2. Lower left contains Figure 2, your symbolic solution for the C obtained with the initial conditions (t0,x0). This is plotted over the slope field.
    3. Lower right contains Figure 3, the numerical solution of $x'=f(t,x)$ determined via sol1 plotted over the slope field.
  9. If the solutions to Figures 2 and 3 do not match, there are generally two possibilities: (a) the existence and uniqueness theorems governing ODEs have been violated; or (b) you did not solve the ODE properly or made an error solving for C. In both cases, the numerical solution from MATLAB is usually correct. In case (a), your symbolic plot should look like the numerical plot except that yours will contain more of the graph. Your graph has ``jumped'' over a discontinuity. There is nothing to change at this time, but you should include an explanation to this effect. If your graph looks differnt in other ways, you are probably in case (b). Check your work and any of Items 1-8 that are appropriate for errors.

    If the solutions to Figures 2 and 3 are the same, print a copy of the display.
  10. In the upper left of your printout, neatly and orderly rewrite the work you did to solve the ODE in Item 1 and to find C.



Subsections
next up previous contents
Next: Example Up: First Order ODEs Symbolically Previous: Example 2   Contents
Michael Renardy
2000-05-12