nmmfiles.html

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

  <a href="integrate/trapezoid.m">trapezoid</a>        Composite trapezoid rule
  <a href="integrate/trapzDat.m">trapzDat</a>         Composite trapezoid rule for arbitrarily spaced discrete data
  <a href="integrate/trapzDatTest.m">trapzDatTest</a>     Verify trapzDat function for different types of input
  <a href="integrate/xemx.m">xemx</a>             Evaluate x*exp(-x), where x is a scalar or vector
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="interact">Contents of interact directory</a></h2>
<pre>
  <a href="interact/Tfield.m">Tfield</a>           Load and plot temperature field stored in custom format text file
  <a href="interact/demoContour.m">demoContour</a>      Demonstration of a contour plot
  <a href="interact/demoSubplot.m">demoSubplot</a>      Demonstration of subplot with four sine functions
  <a href="interact/demoSurfTypes.m">demoSurfTypes</a>    Different types of surface plots for z = 2 - x^2 - y^2
  <a href="interact/sourceSinkPlot.m">sourceSinkPlot</a>   Surface and contour plots of source/sink pair
  <a href="interact/tempPlot.m">tempPlot</a>         Plot monthly temperature variations for Portland, OR.
  <a href="interact/tempStats.m">tempStats</a>        Load data from PDXtemp.dat and compute simple statistics
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="interpolate">Contents of interpolate directory</a></h2>
<pre>
  <a href="interpolate/binSearch.m">binSearch</a>        Binary search to find index i such that x(i)<= xhat <= x(i+1)
  <a href="interpolate/compInterp.m">compInterp</a>       Compare flops for interpolation with different polynomial bases
  <a href="interpolate/compNotKnot.m">compNotKnot</a>      Compare implementations of splines with not-a-knot end conditions
  <a href="interpolate/compSplinePlot.m">compSplinePlot</a>   Compare end conditions for cubic-spline interpolants
  <a href="interpolate/compSplintFlops.m">compSplintFlops</a>   Show flops savings due to sparse storage in spline interp
  <a href="interpolate/demoGasLag.m">demoGasLag</a>       Interpolate gasoline price data with Lagrange polynomials
  <a href="interpolate/demoGasNewt.m">demoGasNewt</a>      Interpolate gasoline price data with Newton polynomials
  <a href="interpolate/demoGasVand.m">demoGasVand</a>      Interpolate gas price data using monomial basis functions
  <a href="interpolate/demoGasVandShift.m">demoGasVandShift</a>   Interpolate gas price data using monomials and shifted dates
  <a href="interpolate/demoHermite.m">demoHermite</a>      Cubic Hermite interpolation of y = x*exp(-x) for 0 <= x <= 8
  <a href="interpolate/demoInterp1.m">demoInterp1</a>      Use built in interp1 function on data sampled from 'humps'
  <a href="interpolate/demoLinterp.m">demoLinterp</a>      Script demonstrating linterp function with thermocouple data
  <a href="interpolate/demoSplineFE.m">demoSplineFE</a>     Spline approx to y = x*exp(-x) with fixed slope end conditions
  <a href="interpolate/demoWiggle.m">demoWiggle</a>       Plot wiggle from increasing order of polynomial interpolant
  <a href="interpolate/divDiffTable.m">divDiffTable</a>     Construct a table of divided-difference coefficients
  <a href="interpolate/hermint.m">hermint</a>          Piecewise-cubic Hermite interpolation
  <a href="interpolate/lagrint.m">lagrint</a>          Interpolation with Lagrange polynomials of arbitrary degree
  <a href="interpolate/linterp.m">linterp</a>          Piecewise linear interpolation in a table of (x,y) data
  <a href="interpolate/newtint.m">newtint</a>          Interpolation with Newton polynomials of arbitrary degree
  <a href="interpolate/splint.m">splint</a>           Cubic-spline interpolation with various end conditions
  <a href="interpolate/splintFE.m">splintFE</a>         Cubic-spline interpolation with fixed slope end conditions
  <a href="interpolate/splintFull.m">splintFull</a>       Cubic-spline interpolation with various end conditions; full matrix storage
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="linalg">Contents of linalg directory</a></h2>
<pre>
  <a href="linalg/Cholesky.m">Cholesky</a>         Cholesky factorization of a symmetric, positive definite matrix
  <a href="linalg/GEPivShow.m">GEPivShow</a>        Show steps in Gauss elimination with partial pivoting and
  <a href="linalg/GEshow.m">GEshow</a>           Show steps in Gauss elimination and back substitution
  <a href="linalg/JfReservoir.m">JfReservoir</a>      Jacobian and f vector for three reservoir system
  <a href="linalg/demoNewtonSys.m">demoNewtonSys</a>    Solve a 2-by-2 nonlinear system by Newton's method
  <a href="linalg/demoSSub.m">demoSSub</a>         Solve a 2-by-2 nonlinear system by successive substitution
  <a href="linalg/demoSparse.m">demoSparse</a>       Script demonstrating some sparse matrix commands
  <a href="linalg/demothreeRes.m">demothreeRes</a>     Solve the three reservoir problem with Newton's method
  <a href="linalg/luNopiv.m">luNopiv</a>          LU factorization without pivoting
  <a href="linalg/luNopivVec.m">luNopivVec</a>       LU factorization without pivoting - vectorized implementation
  <a href="linalg/luPiv.m">luPiv</a>            LU factorization with partial pivoting
  <a href="linalg/newtonSys.m">newtonSys</a>        Newton's method for systems of nonlinear equations.
  <a href="linalg/pumpCurve.m">pumpCurve</a>        Coefficients of quadratic pump curve given head and flow rate
  <a href="linalg/rotvec.m">rotvec</a>           Rotates a three dimensional vector
  <a href="linalg/rotvecTest.m">rotvecTest</a>       Use the rotvec function to rotate some vectors
  <a href="linalg/solveSpeed.m">solveSpeed</a>       Measure elapsed time and flop rate for solving Ax=b
  <a href="linalg/tridiag.m">tridiag</a>          Create tridiagonal matrix from two or three scalars or vectors.
  <a href="linalg/tridiags.m">tridiags</a>         Create sparse tridiagonal matrix from two or three scalars or vectors
  <a href="linalg/vectorSequence.m">vectorSequence</a>   Behavior of a vector sequence x.^k in different p-norms
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="ode">Contents of ode directory</a></h2>
<pre>
  <a href="ode/compEM.m">compEM</a>           Compare Euler and Midpoint for solution of dy/dt = -y;  y(0) = 1
  <a href="ode/compEMRK4.m">compEMRK4</a>        Compare flops and accuracy of Euler, Midpoint and RK4 methods
  <a href="ode/demoEuler.m">demoEuler</a>        Integrate dy/dt = t - 2*y;  y(0) = 1 with Euler's method
  <a href="ode/demoODE45.m">demoODE45</a>        Integrate dy/dx = -y;  y(0) = 1 with ode45
  <a href="ode/demoODE45args.m">demoODE45args</a>    Integrate dy/dt = -alpha*y;  y(0) = 1 with variable alpha
  <a href="ode/demoODE45opts.m">demoODE45opts</a>    Integrate dy/dx = -y;  y(0) = 1 with ode45 and options
  <a href="ode/demoPredprey.m">demoPredprey</a>     Coupled ODEs for a two-species predator-prey simulation
  <a href="ode/demoRK4.m">demoRK4</a>          Integrate  dy/dt = t-2*y;  y(0) = 1 with RK4 method
  <a href="ode/demoSmd.m">demoSmd</a>          Second order system of ODEs for a spring-mass-damper system
  <a href="ode/demoSteel.m">demoSteel</a>        Solve ODE describing heat treating of a steel bar using ode45
  <a href="ode/demoSystem.m">demoSystem</a>       Solve system of two coupled first order ODEs
  <a href="ode/odeEuler.m">odeEuler</a>         Euler's method for integration of a single, first order ODE
  <a href="ode/odeMidpt.m">odeMidpt</a>         Midpoint method for integration of a single, first order ODE
  <a href="ode/odeRK4.m">odeRK4</a>           Fourth order Runge-Kutta method for a single, first order ODE
  <a href="ode/odeRK4sys.m">odeRK4sys</a>        Fourth order Runge-Kutta method for systems of first order ODEs
  <a href="ode/odeRK4sysv.m">odeRK4sysv</a>       Fourth order Runge-Kutta method for systems of first order ODEs
  <a href="ode/odeRK4v.m">odeRK4v</a>          Fourth order Runge-Kutta method for a single, first order ODE
  <a href="ode/rhs1.m">rhs1</a>             Evaluate right hand side of dy/dt = t - 2*y 
  <a href="ode/rhs2.m">rhs2</a>             Evaluate right hand side of dy/dt = -y 
  <a href="ode/rhsDecay.m">rhsDecay</a>         Evaluate rhs of dy/dt = -alpha*y with a variable alpha.
  <a href="ode/rhsPop2.m">rhsPop2</a>          Right hand sides of coupled ODEs for 2 species predator-prey system
  <a href="ode/rhsSmd.m">rhsSmd</a>           Right-hand sides of coupled ODEs for a spring-mass-damper system
  <a href="ode/rhsSteelHeat.m">rhsSteelHeat</a>     Right hand side of first order ODE for heat treating simulation
  <a href="ode/rhsSys.m">rhsSys</a>           Right hand side vector for two, coupled, first order ODEs
  <a href="ode/testSteel.m">testSteel</a>        Verify that solutions obtained by ode45 are independent of
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="orgbug">Contents of orgbug directory</a></h2>
<pre>
  <a href="orgbug/dailyAve.m">dailyAve</a>         Compute average daily flow rate from hourly flow data
  <a href="orgbug/indent.m">indent</a>           Script demonstrating advantages of indentation
  <a href="orgbug/multiply.m">multiply</a>         Compute product of two arguments
  <a href="orgbug/noneg.m">noneg</a>            Returns x if x>0, otherwise prints an error message and stops
  <a href="orgbug/noneg2.m">noneg2</a>           Returns x if x>0, otherwise prints an error message and stops.
  <a href="orgbug/nonegp.m">nonegp</a>           Returns x if x>0, otherwise prints a message and pauses
  <a href="orgbug/quadroot.m">quadroot</a>         Roots of quadratic equation and demo of keyboard command
  <a href="orgbug/riverReport.m">riverReport</a>      Compute and plot summary of river flow data
  <a href="orgbug/riverReport2.m">riverReport2</a>     Compute and plot summary of river flow data.  Compatible
  <a href="orgbug/weeklyAve.m">weeklyAve</a>        Compute average weekly flow rate from hourly flow data
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="program">Contents of program directory</a></h2>
<pre>
  <a href="program/H2Odensity.m">H2Odensity</a>       Density of saturated liquid water
  <a href="program/H2Oprops.m">H2Oprops</a>         Thermophysical properties of saturated liquid water
  <a href="program/addmult.m">addmult</a>          Compute sum and product of two matrices
  <a href="program/cvcon.m">cvcon</a>            Use switch construct to assign constants in curve fit for cv of fluids
  <a href="program/demoArgs.m">demoArgs</a>         Variable numbers of input and output parameters
  <a href="program/demoBreak.m">demoBreak</a>        Show how "break" causes exit from a while loop
  <a href="program/demoCopy.m">demoCopy</a>         Demonstrate vectorized copy operations
  <a href="program/demoLogic.m">demoLogic</a>        Demonstrate logical comparison operators, print results to screen
  <a href="program/demoLoop.m">demoLoop</a>         Script file to demonstrate for and while loops
  <a href="program/demoReturn.m">demoReturn</a>       Show how "return" causes exit from a function
  <a href="program/demoXcosx.m">demoXcosx</a>        Use an inline function object with the fsum function
  <a href="program/easyplot.m">easyplot</a>         Script to plot data in file xy.dat
  <a href="program/elvis.m">elvis</a>            Script file to demonstrate output from fprintf
  <a href="program/fsum.m">fsum</a>             Computes sum of f(x) values at n points in  a <= x <= b
  <a href="program/inputAbuse.m">inputAbuse</a>       Use annoying input messages to compute sum of three variables 
  <a href="program/myCon.m">myCon</a>            Defines useful constants in the workspace
  <a href="program/plotData.m">plotData</a>         Plot (x,y) data from columns of an external file
  <a href="program/plotfun.m">plotfun</a>          Plot sin(x), cos(x), and sin(x)*cos(x) for a prescribed
  <a href="program/polyGeom.m">polyGeom</a>         Compute area and perimeter of a regular polygon
  <a href="program/sincos.m">sincos</a>           Evaluates sin(x)*cos(x) for any input x
  <a href="program/takeout.m">takeout</a>          Script to display restaurant telephone numbers. 
  <a href="program/threesum.m">threesum</a>         Add three variable and returns the result
  <a href="program/trigplot.m">trigplot</a>         Script to plot sin(x), cos(x), and sin(x)*cos(x)
  <a href="program/twosum.m">twosum</a>           Add two matrices and print the result
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="rootfind">Contents of rootfind directory</a></h2>
<pre>
  <a href="rootfind/bisect.m">bisect</a>           Use bisection to find a root of the scalar equation f(x) = 0
  <a href="rootfind/brackPlot.m">brackPlot</a>        Find subintervals on x that contain sign changes of f(x)
  <a href="rootfind/demoBisect.m">demoBisect</a>       Use bisection to find the root of x - x^(1/3) - 2
  <a href="rootfind/demoNewton.m">demoNewton</a>       Use Newton's method to find the root of f(x) = x - x^(1/3) - 2
  <a href="rootfind/fx3.m">fx3</a>              Evaluates f(x) = x - x^(1/3) - 2
  <a href="rootfind/fx3n.m">fx3n</a>             Evaluate f(x) = x - x^(1/3) - 2 and dfdx for Newton algorithm
  <a href="rootfind/legs.m">legs</a>             Evaluate f(theta) for picnic leg geometry.  Used with root-finders.
  <a href="rootfind/legsn.m">legsn</a>            Evaluate f(theta) and fprime(theta) for the picnic leg problem.
  <a href="rootfind/legsnNG.m">legsnNG</a>          Evaluate f(theta) and fprime(theta) for the picnic leg problem.
  <a href="rootfind/legz.m">legz</a>             Evaluate f(theta) for picnic leg geometry.  Use pass-through
  <a href="rootfind/newton.m">newton</a>           Newton's method to find a root of the scalar equation f(x) = 0
  <a href="rootfind/newtonNG.m">newtonNG</a>         Newton's method to find a root of the scalar equation f(x) = 0,
  <a href="rootfind/secant.m">secant</a>           Secant method for finding roots of scalar f(x) = 0
  <a href="rootfind/tablen.m">tablen</a>           Use Newton's method to find dimensions of picnic table legs
  <a href="rootfind/tablenNG.m">tablenNG</a>         Use Newton's method to find dimensions of picnic table legs
  <a href="rootfind/tablez.m">tablez</a>           Use fzero to find dimensions of picnic table legs
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


<br>
<hr width="90%" align=left>
<h2><a name="utils">Contents of utils directory</a></h2>
<pre>
  <a href="utils/addpwd.m">addpwd</a>           Add the current directory to the path
  <a href="utils/chop10.m">chop10</a>           Round a floating point number to n base-10 digits.
  <a href="utils/drawPlane.m">drawPlane</a>        Draws a plane in 3D
  <a href="utils/loadColData.m">loadColData</a>      Import a file containing header text, column titles and data
  <a href="utils/loadColDateData.m">loadColDateData</a>   Import header text, column titles, date and numeric data
  <a href="utils/makeHTMLindex.m">makeHTMLindex</a>    Creates HTML file listing contents of the NMM toolbox
  <a href="utils/myArrow.m">myArrow</a>          Draw 2D arrows with filled tip(s).
  <a href="utils/myArrow3.m">myArrow3</a>         Draw 3D arrows with filled head. Size and color of
  <a href="utils/nmmCheck.m">nmmCheck</a>         Verify installation of NMM toolbox
  <a href="utils/nmmVersion.m">nmmVersion</a>       Returns version number of NMM toolbox currently installed
  <a href="utils/yesNoAnswer.m">yesNoAnswer</a>      Prompt user with a question, and return 1 for 'yes' and 0 for 'no'
</pre>
<p align=center>
<A HREF="#contents"><SMALL>Back to List of Directories</SMALL></A>
</p>


</BODY>
</HTML>