Syllabus:
syllabus
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)
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