# Course Materials for Math/AMath 596: Numerical Solution of Integral Equations

These materials are for Math/AMath 596, taught by Anne Greenbaum in the Spring term of 2001 at the University of Washington.

Syllabus: postscript format, or html

## Assignments and handouts.

First homework assignment due Wed., Apr. 18: postscript format, or html

Second homework assignment due Wed., Apr. 25: postscript format, or html

Third homework assignment due Wed., May 9: postscript format, or html

## Sample MATLAB codes.

Code to solve Laplace's equation inside an ellipse. lap2d_int.m
Uses boundary points that are equally spaced in theta, not arclength, and therefore is only second order accurate. Solves linear system directly.

Code to solve Laplace's equation inside an ellipse. lap2d_int_equal.m
Uses functions: inteval.m and spacing.m
Uses boundary points that are equally spaced in arclength (to within about 1.e-8), and therefore converges superalgebraically to this tolerance. Solves linear system directly.

Code to solve Laplace's equation outside an ellipse. lap2d_ext.m
Uses boundary points that are equally spaced in theta, not arclength, and therefore is only second order accurate. Solves linear system using MATLAB routine GMRES.

Code to implement Mayo's method of computing the solution to Laplace's equation on points throughout a lattice, by first solving an integral equation to determine jumps in the solution and its derivatives and then using a fast Poisson solver on an embedding rectangle. Code written by Michal Skokan. The main code is p1.m . It uses the following files: f.m , fprime.m , f_sec.m , mgres.m , mu_nul.m , mu_prime.m , mu_sec.m , mu_xx.m .