Computational Physics at Carleton University Index
Course Information Numerical Methods Monte Carlo Methods Statistics for Physicists Special Topics
Root Finding


Root Finding

The general problem of solving a nonlinear equation can be reexpressed as a problem of finding points where a function is zero. This section studies various methods for zeroing in on the root of the problem as quickly as possible. In one dimension, bracketing of the root is possible, and Ridders and Brent methods are good general purpose methods that converge on the solution very quickly. Roots of polynomial equations can be found using Laguerre's method.

The Newton-Raphson method requires knowledge of the function and its first derivatives. A diversion into fractals generated with the approach is included. The method can be used for multidimensional problems.

The exercise involves finding the maximum of a blackbody radiation spectrum.

Reference:

Numerical Recipies, chapter 9.

By clicking the above thumbnail page images, you can peruse the slides. You can also download the (1.0M) colour postscript file that contains these slides.

The entire set of slides for this course is available in colour postscript or in black and white postscript.