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

Course Notes for Math 465

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