Physics 381C: Non-Linear Equations & Newton's Method
Please report all errors/typos. etc to
matt@infeld.ph.utexas.edu
- General Reference
- Lectures 2 and 3 (Feb 8, 11)
- newt2.f:
Fortran program for computing root of simple set of non-linear
equations (d=2) discussed in class.
Makefile
and sample
output on SGIs.
A check on the
computation using Maple's numerical root finding capabilities
(fsolve).
- dgesv.f:
LAPACK driver routine for solving general linear system via
LU decomposition with partial pivoting. Used in Newton
iteration to determine unknown update vector from the
Jacobian matrix and residual vector.
- Building and running
newt2 on the Cray J90, charon.cc.utexas.edu.
Illustrates use of script
f77transcray,
which converts real*8 SGI Fortran to single precision
(also 8-byte) Cray Fortran.
-
LAPACK User's Guide (html). We will discuss
and use LAPACK further in weeks to come.
-
LAPACK Source Code (browse directory)
- Supplementary Material
- A very short blurb on
Raphson.
Thanks to Ilian
for the reference.