Math 5334, Spring 2009
Numerical Analysis for Scientists & Engineers

Previous Announcements:

• 4/10/09: The deadline for HW6 has been moved to Thursday 4/16 by 8:00 a.m.

• 4/10/09: Hint on HW6.
In the bird chirp problem, you can play the chirps (and several other sound files) with the demo "xpsound". If you "edit xpsound" and search for "chirp", you can see how to load the chirp data in order to do your FFT analysis. In particular:

load('chirp','Fs','y'); % to get chirp data and sampling freq.
y=y/max(abs(y));t=(0:length(y)-1)/Fs;
plot(t,y) % to display the chirps in time
sound(y,Fs) % to listen to all of the chirps

You will need to break up the chirp data so that you can do an FFT on each chirp separately.

• 4/6/09: Homework 6 is posted below. Due Tuesday 4/14 by 8:00 a.m.

• 3/31/09:Suggestions for HW5
• B&F 11.3 #3: in the linear FD algorithm in B&F, notice that steps 4 through 8 are solving the tridiagonal system. I recommend replacing these steps with backslash or tridisolve.
• B&F 11.4 #1: the nonlinear FD algorithm in B&F claims to need N >=2. I don't see any reason it won't work with N=1.
As in the previous problem, I recommend replacing steps 8 through 12 of this algorithms with backslash or tridisolve.

• 3/31/09: Requests for homework assignments in general:
• Make sure your name appears in every file you send me. When I print out your files to grade, your name should be (at least) on the front page of every printout. In Matlab code, include your name in a comment near the beginning of the code.
• Please make every effort to limit the length of the HW assignment you turn in. It is not necessary to print out every vector or matrix produced by Matlab unless that vector happens to be the answer to the problem, etc. (Use the semicolon in Matlab to suppress output unless that output is the answer.) Do not just dump all of your Matlab output into a file and send me 45 pages or more. Format your answers and make it a reasonable sized document. The problem numbers and your answers should be clearly marked and easy to find when I'm grading. The trees and my printer cartridges will thank you.
• Please put the problems in order.

• 3/31/09: For the projects, note that it is OK to work in small groups as long as each member contributes a substantial amount to the project. Let me know who you are working with when you send me your topics on Wednesday.

• 3/31/09:Update to course grading: I said below that homework is 60% and the projects are 40% of your grade. I am amending that to have whichever (of HW and the project) is your highest grade counts 60% and whichever is your lowest counts 40%. So you can make up for some lower HW scores if you need to by extra work on your project.

• 3/25/09:
• HW 5 is posted below. It is due by Thursday, 4/2, at 8:00 a.m.

• Here is information on the course project.
Note that since we are not having the in-class exams originally listed in the syllabus that homework is 60% of the grade and the project is 40%.
Please email me with your choice of project topic by Wednesday, 4/1.

• 3/9/09: No take-home exam: After further consideration, I decided that a take home exam is not necessary since there is already sufficient homework (plus the upcoming projects) in the course. So we do not have a take home exam this week.

• 3/9/09: HW4 solutions are posted below.

• 3/6/09: Recall that we have our take-home midterm exam next week. Since we voted in class several times to have the exam before rather than during Spring break, I will hand out the exam in class on Monday, 3/9, and it will be due in class Friday, 3/13.

The exam will cover course material up to (but not including) quadrature.

HW3 solutions are posted below. I am still receiving HW4 submissions, so I will post solutions to HW4 after the 5 possible late days have passed. So I will post HW4 solutions this weekend.

I will have HW3 and HW4 to return to you on Monday.

• 2/22/09: HW4 is posted below. Due Tuesday 3/2 by 8:00 a.m.

• 2/14/09: HW1 is graded and will be returned Monday. Solutions for HW2 are posted below.

• 2/13/09: HW3 deadline has been extended. Please submit HW3 by 8:00 a.m. Tuesday 2/17.
• Note that the Burden & Faires problems are generally meant to be done by hand, not with Matlab (and definitely not with backslash). You are certainly free to use Matlab as a calculator or to use a calculator.
• In B & F 6.1 #4, if you prefer, you can write your own Matlab function to do Gaussian Elimination (LU) without pivoting and forcing each calculation to use only two-digits of accuracy. However, for a 3x3 problem, it is probably faster to do it by hand (with a calculator and round each step to two digit accuracy).
• In B & F 2.2 #6, you can do the problem by hand or write a Matlab function to do fixed point iteration on the given function.
• In all cases, if you write Matlab code, turn in the code with your homework solutions.

• 2/5/09: HW3 is posted below.

• 2/3/09: HW1 solutions are posted below.

• 1/25/09: I removed the Burden & Faires problems from HW2 since the book is not available to everyone yet. For Friday, just turn in the problems from the Moler book.
Please check Varsity books and Barnes & Noble for the B&F book on Tuesday.
Monday we will finish up sparsity (Moler section 2.10) and begin interpolation (Moler chapter 3).
For anyone who joined our class late, Monday (1/26) midnight is the deadline for turning in HW1. After that I will post the solutions.

• 1/23/09: HW2 is posted below. Due next Friday, 1/30, by midnight.

• 1/16/09: HW 1 deadline has been extended until midnight tonight. Please email me your homework solutions by midnight (unless you are one of the newly added students that I gave an extension)

• 1/14/09: For your homeworks, please write up your solutions as a text file, or if you prefer something like Word, please save it as a PDF file (I don't have Word on my computer). You can email me your HW. Please include copies of any matlab files (.m files) that you wrote. Thank you.

• 1/13/09: Don't forget HW1 due this Friday by 4:00.
You should also be reading the Moler book (chapters 1 & 2) and the first couple of sections of Driscoll's "A Crashcourse in MATLAB"

Last updated: 4/29/2009