contents.m

这是《Numerical Methods with MATLAB: Implementation and Application》一书的配书程序（Matlab）

%
% binSearch         Binary search to find index i such that x(i)<= xhat <= x(i+1)
% binSearch.c       C source code for binSearch.mex
% compInterp        Compare flops for interpolation with different polynomial bases
% compNotKnot       Compare implementations of splines with not-a-knot end conditions
%                   NMM toolbox routine splint, is compared with built in spline
%                   function in the interpolation of y=x*exp(-x) on [0,5]
% compSplintFlops   Show flops savings due to sparse storage in spline interp
% compSplinePlot    Compare cubic splines w/ different end conditions
%                   Approximations to y = x*exp(-x) are constructed and plotted
% demoGasLag        Interpolate gasoline price data with Lagrange polynomials
% demoGasNewt       Interpolate gasoline price data with Newton polynomials
% demoGasVand       Interpolate gas price data using monomial basis functions
% demoGasVandShift  Interpolate gas price data using monomials and shifted dates
% demoHermite       Piecewise cubic Hermite interpolation of y = x*exp(-x) on [0,5]
% demoInterp1       Use built in interp1 function on data sampled from 'humps'
% demoLinterp       Script demonstrating linterp function with thermocouple data
% demoSplineFE      Spline approx to y = x*exp(-x) with fixed slope end conditions
% demoWiggle        Plot wiggle from increasing order of polynomial interpolant
% divdiffTable      Construct a table of divided difference coefficients
% hermint           Piecewise cubic Hermite interpolation
% lagrint           Interpolation with Lagrange polynomials of arbitrary degree
% linterp           Piecewise linear interpolation in a table of (x,y) data
% newtint           Interpolation with a Newton polynomias of arbitrary degree
% splint            Cubic spline interpolation with various end conditions
% splintFE          Cubic spline interpolation with fixed slop end conditions
% splintFull        Cubic-spline interpolation with various end conditions; full matrix storage
%
% ----------------------------------------------------------------
%
% Contents of "linalg" directory in the NMM toolbox
%
% Cholesky        Cholesky factorization of a symmetric, positive definite matrix
% demoNewtonSys   Solve a 2-by-2 nonlinear system by Newton's method
% demoSparse      Script demonstrating some sparse matrix commands
% demoSSub        Solve a 2-by-2 nonlinear system by successive substitution
% GEPivShow       Show steps in Gauss elimination with partial pivoting and
%                 back substitution
% GEshow          Show steps in Gauss elimination and back substitution
%                 No pivoting is used.
% luNopiv         LU factorization without pivoting
% luNopivVec      LU factorization without pivoting - vectorized implementation
% luPiv           LU factorization with partial pivoting
% newtonSys       Newton's method for systems of nonlinear equations.
% pumpCurve       Coefficients of quadratic pump curve given head and flow rate
% rotvec          Rotates a three dimensional vector
% rotvectTest     Use the rotvec function to rotate some vectors
% solveSpeed      Measure elapsed time and flop rate for solving Ax=b
% tridiag         Create tridiagonal matrix from two or three scalars or vectors.
% tridiags        Create sparse tridiagonal matrix from two or three scalars or vectors.
% vectorSequence  Behavior of a vector sequence x.^k in different p-norms
%
% ----------------------------------------------------------------
%
% Contents of "ode" directory in the NMM toolbox
%
% compEM       Compare Euler and Midpoint for solution of dy/dx = -y;  y(0) = 1
% compEMRK4    Compare flops and accuracy of Euler, Midpoint and RK4 methods
%              for the solution of  dy/dx = -y;  y(0) = 1
% demoEuler    Integrate dy/dx = x^2 - y;  y(0) = 1 with Euler's method
% demoODE45    Integrate dy/dx = -y;  y(0) = 1 with ode45
% demoRK4      Integrate  dy/dx = -y;  y(0) = 1 with RK4 method
% demoSystem1  Driver for integration of system defined in rhsSys1
% odeEuler     Euler's method for integration of a single, 1st order ODE
% odeMidpt     Integrate a single ODE with the Midpoint Method
% odeRK4       Integrate a single ODE with the fourth order Runge-Kutta method
% odeRK4sys    Fourth order Runge-Kutta method for systems of first order ODEs
%              Non vectorized version
% odeRK4sysv   Use 4th order RK4 method to integrate a system of first order ODEs
%              Vectorized version.
% odeRK4v      Integrate a single ODE with the fourth order Runge-Kutta method
%              Allows pass through of parameters to ODE-defining routine
% predprey     Integrate coupled ODEs for a two-species predator-prey simulation
% rhs1         Evaluate right hand side of dy/dx = x^2 - y 
% rhs2         Evaluate right hand side of dy/dx = -y 
% rhsDecay     Define rhs of dydx = -alpha*y with a variable alpha
% rhspop2      Evaluate rhs of the ODEs for 2 species predator-prey system
% rhssmd       Evaluate right hand sides of ODEs for a spring-mass-damper system
% rhsSteelHeat Evaluate right hand side of ODE for heat treating simulation
% rhsSys1      Evaluate rhs of the system
%              dy1/dx = -y1*exp(1-x) + 0.8*y2,  dy2/dx = y1 - y2^3
% smdsys       Solve 2nd order system of ODEs for a spring-mass-damper system
% steelHeat    Set up and solve ODE describing heat treating of a steel bar.
% steelTest    Script to verify that solutions obtained by ode45 are
%              independent of the interpolative refinement of the solution.
%
% ----------------------------------------------------------------
%
% Contents of "orgbug" directory in the NMM toolbox
%
% dailyAve      Compute average daily flow rate from hourly flow data
%               Used for analysis of Glen Canyon Dam data
% indent        Script demonstrating advantages of indentation
% multiply      Compute product of two arguments
% noneg         Returns x if x>0, otherwise prints an error message and stops
% noneg2        Returns x if x>0, otherwise prints an error message and stops.
%               The error message contains the value of x.
% nonegp        Returns x if x>0, otherwise prints a message and pauses
% quadroot      Roots of quadratic equation and demo of keyboard command
% riverReport   Compute and plot summary of river flow data
% riverReport2  Compute and plot summary of river flow data.  Compatible with
%               Student version of MATLAB 5.
% weeklyAve     Compute average weekly flow rate from hourly flow data
%               Used for analysis of Glen Canyon Dam data
%
% ----------------------------------------------------------------
%
% Contents of "program" directory in the NMM toolbox
%
% addmult      Compute sum and product of two matrices
% binSearch.c  C source code for binSearch.mex --- see also binSearch.m
%              in Chapter 9, "Interpolation"
% cvcon        Use switch construct to assign constants in curve fit for cv of fluids
% demoArgs     Variable numbers of input and output parameters
% demoBreak    Show how the break command causes exit from a while loop
% demoCopy     Demonstrate vectorized copy operations
%              Messages describing various copy operations are given
%              along with the results of those operations.  See text
%              in Appendix A for additional description.
% demoLogic    Function m-file file to demonstrate logical operators
% demoLoop     Script file to demonstrate for and while loops
% demoReturn   Show how the return command causes exit from a function
% demoXcosx    Demonstrate use of an inline function object with the fsum function
% easyplot     Script to plot data in file xy.dat
% elvis        Script file to demonstrate output from fprintf
% fsum         Computes sum of f(x) values at n points in  a <= x <= b
% H2Odensity   Density of saturated liquid water
% H20Props     Thermophysical properties of saturated liquid water
% inputAbuse   Use annoying input messages to compute sum of three variables 
% myCon        Script to define useful constants in the workspace
% plotData     Plot (x,y) data from columns of an external file
% plotfun      Plot sin(x), cos(x), and sin(x)*cos(x) for a prescribed
% polyGeom     Compute area and perimeter of a regular polygon
% sincos       Evaluates sin(x)*cos(x) for any input x
% takeout      Script to display restaurant telephone numbers. 
% threesum     Adds three variables and returns the result
% trigplot     Script to plot sin(x), cos(x), and sin(x)*cos(x)
% twosum       Adds two variables and prints the result
%
% ----------------------------------------------------------------
%
% Contents of "rootfind" directory in the NMM toolbox
%
% bisect      Use bisection to find a root of the scalar equation f(x) = 0
% brackPlot   Find and plot brackets for roots of a function.
% demoBisect  Use bisection to find the root of x - x^(1/3) - 2
% demoNewton  Use Newton's method to find the root of f(x) = x - x^(1/3) - 2
% fx3         Evaluates f(x) = x - x^(1/3) - 2
% fx3n        Evaluate f(x) = x - x^(1/3) - 2 and dfdx for Newton algorithm
% legs        Evaluate f(theta) for picnic leg geometry.  Used with root-finders.
% legsn       Evaluate f(theta) and fprime(theta) for the picnic leg problem.
%             Used with the Newton's method.
% legsnNG     Evaluate f(theta) and fprime(theta) for the picnic leg problem.
%             without using global variables.  Used with newtonNG
% legz        Evaluate f(theta) for picnic leg geometry.  Use pass-through
%             parameters in fzero to send w, h, and b values to this function.
%             Used with fzero
% newton      Newton's method to find a root of the scalar equation f(x) = 0
% newtonNG    Newton's method to find a root of the scalar equation f(x) = 0
%             This version can pass optional input arguments to the f(x) function
% secant      Shell of the secant function.  Implementation left to the reader
% tablen      Use Newton's method to find dimensions of picnic table legs
% tablenNG    Use Newton's method to find dimensions of picnic table legs
%             Variable input arguments are passed through to newtonNG
% tablez      Use fzero to find dimensions of picnic table legs
%
% ----------------------------------------------------------------
%
% Contents of "utils" directory in the NMM toolbox
%
% addpwd           Add the current directory to the path
% chop10           Round a floating point number to n base-10 digits.
% drawPlane        Draws a plane in 3D
% loadColData      Import a file containing header text, column titles and data
% loadColDateData  Import a file containing header text, column titles and data.
%                  First column of data is a date in one of MATLAB's datenum formats
% makeHTMLindex    Creates HTML file listing contents of NMM toolbox
% myArrow          Draw 2D arrows with filled tip(s).
% myArrow3         Draw 3D arrows with filled head. Size and color of
%                  arrowhead can be specified
% nmmCheck         Verify installation of NMM toolbox
% nmmVersion       Returns version number of NMM toolbox currently installed
% yesNoAnswer      Prompt user with a question, and return 1 for 'yes' and 0 for 'no'