350 Handouts and M-files

• Numerical Computing with MATLAB (OUR TEXTBOOK by Cleve Moler which can be downloaded for FREE)
• Getting Started in MATLAB (pdf) (some very basic info to get you off the ground)
• A Very Elementary MATLAB Tutorial from The MathWorks
• Experiments with MATLAB (by Cleve Moler), the textbook for MATH 100 which can be downloaded for FREE.
• For more introductory material and tutorials see the main course page

• Slides for lectures.
  The printer-friendly versions collapse all overlays.
  Advantages: fewer pages, smaller files;
  Disadvantage: a couple of overlays are missing because they print on top of each other in the presentation.
  Chapter 1 Chapter 1 (print) (minor updates, Jan.18)
  Chapter 2 Chapter 2 (print) (some links updated, Jan.26)
  Chapter 3 Chapter 3 (print) (minor updates, Feb.20; Maple maplet and Mathematica notebooks for Mairhuber-Curtis theorem can be found below)
  Chapter 4 Chapter 4 (print) (minor updates, March 2)
  Chapter 5 Chapter 5 (print)
  Chapter 6 Chapter 6 (print) (updated March 27, minor update March 31)
  Chapter 7 Chapter 7 (print) The PDF file contains some embedded VBS scripts that run MATLAB code directly from the PDF. In order to be able to do this you will have to be on an MS Windows machine and download the VBS script file CallMATLAB.vbs into the same folder as the PDF file. You will also need to have MATLAB installed on that machine. Alternatively, you can just cut and paste the MATLAB code fragments into the MATLAB command window.

• MATLAB scripts for the first skydive model:
  • SkydiveDemo.m SkydiveDemo.html
  • Skydive.m
  • Euler.m
• MATLAB scripts for the second skydive model (also uses Euler.m):
  • Skydive2Demo.m Skydive2Demo.html
  • Skydive2.m
• MATLAB scripts for the predator-prey model:
  • LynxHareDemo.m LynxHareDemo.html
  • LynxHare.m
• Maple worksheet on Taylor approximation: Taylor.mw
• Minimal modification of Moler's floatgui.m: Myfloatgui.m (updated Jan.27)
• MATLAB script to demonstrate roundoff error: RoundoffDemo.m
• MATLAB script to illustrate loss of significant digits: LossOfSignificanceDemo.m LossOfSignificanceDemo.html
• Introductory Matlab scripts
  • matrix_vector.m matrix_vector.html (definition of matrices and vectors)
  • arithmetic.m arithmetic.html (simple arithmetic with matrices and vectors)
  • lin_sys.m lin_sys.html (solving linear systems)
  • submatrices.m submatrices.html (definition of submatrices)
  • statistics.m statistics.html (timing in Matlab)
  • loops1.m loops1.html (for-loops, and why to avoid them)
  • loops2.m loops2.html (more on for-loops)
  • plotting.m plotting.html (simple plotting)
  • efficiency.m efficiency.html (some tricks for efficient coding)
  • user_defined_functions.m user_defined_functions.html (how to use user-defined functions)
  • circle.m (a user-defined function)
  • anonymous_functions.m anonymous_functions.html (how to use anonymous functions)
  • function_functions.m function_functions.html (how to use anonymous functions)
• MATLAB script to illustrate use of inverse for division: division.m
• Maple worksheet on partial pivoting: PartialPivoting.mw
• MATLAB scripts from Numerical Computing with MATLAB (NCM) to illustrate the solution of linear systems:
  bslashtx.m
  lutx.m
• MATLAB scripts for the solution of tridiagonal linear systems:
  tridisolve.m (from NCM)
  TridisolveDemo.m TridisolveDemo.html (a little script that calls tridisolve)
• MATLAB scripts for pagerank calculations:
  pagerank.m (from NCM, with a small modification by me)
  TinyWeb.m (a little script that does a simple pagerank calculation and also calls pagerank)
• MATLAB scripts for polynomial interpolation:
  interpgui.m (from NCM)
  InterpolationDemo.m (a little script that calls interpgui)
  polyinterp.m (from NCM)
  PolyinterpDemo.m PolyinterpDemo.html (a little script that calls polyinterp)
• MATLAB script for piecewise linear interpolation:
  piecelin.m (from NCM)
  PiecelinDemo.m (a little script that calls piecelin)
• MATLAB script for piecewise cubic Hermite interpolation:
  PchipDemo.m PchipDemo.html (a little script that calls the built-in pchip function)
• MATLAB scripts for cubic spline interpolation:
  SplineDemo.m SplineDemo.html (a little script that calls the built-in spline function)
  splinetx_natural.m (a modification of NCM's splinetx.m for cubic natural splines)
• Maplet for the proof of the Mairhuber-Curtis theorem:
  MairhuberCurtis.maplet (to be run with the Maplet viewer)
  MairhuberCurtis.mws (Maple source code)
• Mathematica version of the proof of the Mairhuber-Curtis theorem:
  MairhuberCurtis.nb (Mathematica notebook)
• MATLAB scripts for bisection method:
  BisectDemo.m BisectDemo.html (a little script that calls bisect.m)
  bisect.m (bisection method as explained in NCM)
• MATLAB scripts for Newton's method:
  NewtonDemo.m NewtonDemo.html (a little script that calls newton.m)
  newton.m (Newton's method as explained in NCM)
• Maple worksheet for graphical illustration of Newton's method:
  NewtonDemo.mw
• MATLAB scripts for the secant method:
  SecantDemo.m SecantDemo.html (a little script that calls secant.m)
  secant.m (secant method as explained in NCM)
• Maple worksheet for graphical illustration of secant method:
  SecantDemo.mws
• MATLAB scripts for the inverse quadratic interpolation method:
  IQIDemo.m IQIDemo.html (a little script that calls iqi.m)
  iqi.m (inverse quadratic interpolation method as explained in NCM)
• MATLAB scripts for zero finding:
  FzeroDemo.m FzeroDemo.html (a little script that calls fzerotx.m)
  fzerotx.m (zero finding routine from NCM)
• MATLAB scripts for multivariate Newton's method:
  NewtonmvDemo.m NewtonmvDemo.html (a little script that calls newtonmv.m)
  newtonmv.m (multivariate Newton's method)
  missile_f.m (function for example 1)
  missile_j.m (Jacobian for example 1)
  circhyp_f.m (function for example 2)
  circhyp_j.m (Jacobian for example 2)
• MATLAB scripts for function minimization:
  FminDemo.m FminDemo.html (a little script that calls fmintx.m)
  fmintx.m (minimization routine from NCM)
• MATLAB scripts for least squares fitting:
  LSQquad.m LSQquad.html (script that shows it isn't good to solve the normal equations)
  LSQquadOrtho.m LSQquadOrtho.html (script that shows orthogonal basis does better)
  PolyfitDemo.m PolyfitDemo.html (script that calls built-in polyfit and polyval)
  censusgui2010.m (from NCM - updated for 2010 Census, illustrates different curve fitting methods)
  LSQquadQR.m LSQquadQR.html (script that shows how to use QR decomposition, minor update 3/18)
  LSQquadSVD.m LSQquadSVD.html (script that shows how to use SVD, minor update 3/24)
  SVD_movie.m (script that shows how to use SVD for image compression)
  imagesvd.m (similar function from NCM)
• Maple worksheet for various basic numerical integration methods:
  Integration.mw
• MATLAB scripts for numerical integration:
  QuadDemo.m QuadDemo.html (script that calls NCM file quadtx and also uses besselintegrand.m)
  besselintegrand.m (function with integrand called by QuadDemo)
  quadtx.m (from NCM, textbook version of built-in function quad)
  quadgui.m (from NCM, illustrates adaptive quadrature)
  QuadgkDemo.m (compares MATLAB functions quad and quadgk on a few integrals)
  • Maple worksheet for various basic numerical differentiation methods:
    Differentiation.mw
  • MATLAB scripts to illustrate Euler's method:
    • EulerDemo350.m EulerDemo350.html (script that calls Euler.m)
    • Euler.m
  • Maple worksheet for Euler's method:
    EulerDemo.mw
  • MATLAB scripts to illustrate second-order Runge-Kutta method:
    • RK2Demo.m RK2Demo.html (script that calls RK2.m)
    • RK2.m
  • MATLAB scripts to illustrate fourth-order Runge-Kutta method:
    • RK4Demo.m RK4Demo.html (script that calls RK4.m)
    • RK4.m
  • MATLAB script to illustrate the convergence behavior of Euler, second- and fourth-order Runge-Kutta methods:
    • EulerRKConvergenceDemo.m EulerRKConvergenceDemo.html (calls Euler.m, RK2.m and RK4.m)
  • MATLAB scripts for the twobody problem:
    • TwobodyDemo.m TwobodyDemo.html (script that calls twobody.m)
    • twobody.m
  • Maple worksheet for match lighting problem:
    MatchDemo.mw
  • MATLAB script for match lighting problem:
    • MatchDemo.m MatchDemo.html
    • flame.m animated solution from NCM
  • MATLAB scripts for the second skydive model with event handling:
    • Skydive3Demo.m Skydive3Demo.html (script that calls Skydive3.m and uses Skydive3Event for event handling in ode23)
    • Skydive3.m
    • Skydive3Event.m
  • MATLAB scripts for van der Pol problem:
    • VanderPolDemo.m VanderPolDemo.html
    • VanderPolPlots.m VanderPolPlots.html
  Review for Tests
  • Review for Midterm 1:
    Midterm exam from Spring 08 semester (note that only problems 1, 2, 3 and 7 are relevant for our first midterm)
    Solutions for S08 midterm
  • Review for Midterm 2:
    Look at the rest of the problems from the S08 midterm above
    Extra review problems for Midterm 2
  • Review for Final:
    Final exam from Spring 08 semester
    Solutions for S08 final
  Solutions of Tests
  • Solutions for Midterm 1
  • Solutions for Midterm 2

  ---

  Last updated: April 21, 2011
  Greg Fasshauer (fasshauer at iit.edu)