# Course Materials for Math 465: Numerical Analysis II

These materials are for Math 465,
taught by Anne Greenbaum in the Winter term of 2009 at the
University of Washington.
Syllabus:
syllabus

## Assignments and handouts.

HW1: Page 23, Problem 1; pages 38-43, Problems 2(a,b,e,f), 6, 11, and 12. (Due Fri., Jan. 16)

HW2: Pages 69-70, Problems 1, 2, 5, and 7. (Due Fri., Jan. 23)

HW3: Page 220, Problems 1, 2, 3, and 4. (Due Fri., Jan. 30)

HW4: Page 279, Problems 1, 2, 3, and 4. (Due Fri., Feb. 6)

HW5: Page 279, Problems 5, 6, 7, and 13. (Due Fri., Feb. 13)

HW6: Page 279-281, Problems 8, 10, 11, and 14. (Due Mon., Feb. 23)

Information on projects:
projects

Project proposals due: Wed., Feb. 18.

MIDTERM: Wed., Feb. 25.

HW7: Chapter on More Numerical Linear Algebra (p. 264 in evpages,
p. 305 in entire text), Problems 1, 2, 3, and 4. (Due Fri., Mar. 6)

HW8: Chapter on More Numerical Linear Algebra (p. 264 in evpages,
p. 305 in entire text), Problem 5. (Due Fri., Mar. 13)

## Sample MATLAB codes.

Matlab code for Yoda: yoda.m .

Data for Yoda: yoda_low.mat .

Runge-Kutta code for solving planetary motion problem 13 on p. 280-281.
rk_hw5.m

To solve problem using ode45, download function fty.m:
fty.m

and in Matlab type: [tout,yout] = ode45('fty', [tmin tmax], y0) ,
where tmin is the initial time, tmax is the final time, and y0 is the vector
of initial values. Output array yout will contain 4 columns, with the 4
components of y at each of the times in tout.

Following are some solution plots from ode45:

hw5_plot1.pdf

hw5_plot2.pdf

Runge-Kutta code for solving stiff system in problem 14 on p. 281-282.
rk_hw6.m

Following is some sample output:

hw6.eps