Next: Contents
Cornell University
Department of Physics
Phys 480/680, Astro 690 February 18, 2002
Computational Physics, Spring 2002
Homework Assignment # 3
(Due Tuesday, February 26 at 8:40 am in computer lab.)
Agenda and readings for the week of February 12:
Readings marked NR are from Numerical Recipes: The Art of
Scientific Computing, 2
edition (in C). Readings
marked LN are from the course lecture notes to be found at
http://www.muchomas.lassp.cornell.edu/P480.
Goal: This problem set develops a fully automated solver which finds
all of the energy states for any specified external potential ``V[]''
to high precision.
Readings marked NR are from Numerical Recipes: The Art of
Scientific Computing, 2
edition (in C). Readings
marked LN are from the course lecture notes to be found at
http://www.muchomas.lassp.cornell.edu/P480.
Recap of last week:
- Lec6, 02/12 (Tue: mini-lecture, 25mins) Runge-Kutta; Overview of
ODE solvers.
Reading: NR 16.0-16.1. - Lab1, 02/12 (Tue): Change of variables; Scaling of errors for
, and for periodic functions - Lec7, 02/14 (Thu): Derivation of Kohn-Sham equations,
Schrödinger's equation for atoms (spherically symmetric
solutions); boundary value problems; standard form for root finding;
bisection
Reading: NR 17.0-17.1, 9.0-9.1
Preview of coming week:
- Lec8, 02/19 (Tue: mini-lecture, 25mins) Overview of root
finding; secant method; convergence rates
Reading: NR 9.2 - Lab2, 02/19 (Tue): Convergence of solution of ODE's; numerical
solutions of Schrödinger's equation.
- Lec9, 02/21 (Thu): Ridder's method for roots; bracketing of
roots for Schrödinger's equation; numerical interpolation
Reading: NR 9.2, 3.0-3.2 - Lec10, 02/26 (Tue: mini-lecture, 25mins): Numerical stability of
integrating ODEs
Reading: NR 17.2 - Lab3, 02/23: Spectrum of hydrogen, convergence, stability of solutions,
effectiveness of root finders
Tomas Arias
Mon Feb 18 16:02:23 EST 2002