Contents

• Extending ``rk4'' and ``derivs''
  • Coding
  • Debugging
• Solving Poisson's equation
  • Kernel for Poisson's equation
  • Poisson's equation solver
  • Solve for the potential of the 1s state of hydrogen
• HW#2 ADDENDUM
  • Double precision
    • Integration routines
    • ``rk4ind''
    • Poisson equation solver
  • ``User friendly'' code
    • Scaling with Nu
    • Poisson's equation for the 1s state of hydrogen
  • Schrödinger's equation
• About this document ...

Note: For the programs in this problem set we face a small notational difficulty. We have been using the variable n(r) to represent the number of electrons per unit volume throughout space. It would make sense then to define an array n[] to hold the values of the electron density. However, n is so commonly used as an indexing variable (particularly in Numerical Recipes) that this is likely to cause difficulties. Thus, we shall use the array rho[] to represent the electron density n(r).