Professor Gilbert Strang
Office: E17420
email: gs@math.mit.edu
Lecture:
 Tuesday and Thursday, 1011AM in 1190
 Tuesday, 34PM in 54100
Teaching Assistants:
 Jason Choi (choij@mit.edu)
 JanChristian Huetter (huetter@mit.edu)
 Florent Bekerman (bekerman@math.mit.edu)
 Nathan Reid Harman (nharman@mit.edu)
Exam Dates:
 Tuesday, Sep. 30th (26100)
 Thursday, Oct. 30th (Walker 50340)
 Thursday, Dec. 4th (Walker 50340)
Help Sessions:
 Monday 34PM in E17401Q (Office Hour with Florent Bekerman)
 Wednesday 45PM in 66168 (Help Session)
Announcements:
 The first help session will take place on Wednesday, September 10, from 45PM in 66168. Please bring your questions! In addition, Florent will offer an office hour from 34PM every Monday, from 34PM.
 The course Piazza is now up! Please register by following this link: http://piazza.com/mit/fall2014/18085
 First homework has been posted! This will be due next Thursday, Sept. 11 at the start of class. Please write your NAME neatly on the assignment. Include any necessary compute printouts. Collaboration is encouraged, but copying is not allowed.
Homework: [Fall 2014]
#1 (Due: Thursday, Sep. 11)
 A = [1 3 2] is a 1 by 3 matrix and A' = 3 by 1 transpose
 Find K = A' * A and explain why this K is singular (many answers, all OK)
 Find all solutions to Ku=0
 Compute K^4 = A'AA'AA'AA'A by hand in a smart order
 Page 16 of the text multiplies our 2nd difference matrix K (make it infinite to avoid boundary effects) by vectors ones = 1,1,1,.. linear 0,1,2,3,.. squares 0,1,4,9,...
 Do the same for cubes 0,1,8,27,.. and quartic 0,1,16,81,..
 Which ones give the wrong derivative and what is the error
 Multiply these vectors by the FORWARD 1st difference matrix D+ and the CENTERED first difference matrix D_0 (infinite) and compare to the derivatives of x^2, x^3, x^4.
 Solve u'' = delta(xa) with u(0)=0 and u'(1)=0 (free). Graph u(x) and explain what is happening at x=a between 0,1
 Problem 1.2.9 on page 24 is about the 4th difference with coefficients 1,4,6,4,1. Where did those numbers come from and should you divide by a power of h=delta x ?
 Problem 1.2.19 on page 25 is about u'' + u' = 1. Choose a finite difference approximation. Solve for U.
 Invert our T matrix with T_11=1 and diagonals 1,2,1 for n=3 and n=4. From the pattern guess a formula for any n. Can you show that your guess is the correct inverse (any n)?
Homework Notes:
 PLEASE PRINT YOUR NAME CLEARLY
 Good to underline your last name.
 The Class List with class numbers will be created from these homeworks.
 Homeworks are for learning. Discussion is OK. Write clear answers for the graders: please.
 Problems will often come from the 18.085 textbook on CSE. (Printouts / Graphs for MATLAB questions.)
 Section 1.1 is available on math.mit.edu/cse
 A printing of 1.1.27 has a typo (my name for a mistake)

These are the simplest corrections: Verify that
 K comes from A0 A0 ' (prime = transpose)
 T comes from A1 A1 '
 B comes from A2 A2 '
 (Those have 2 then 1 then 0 boundary conditions. They change from positive definite to semidefinite.)
 problem set 1.2 from CSE textbook
Notes from Class:
 Solution to Ku=F is u = K\F backslash in Matlab, K first !
 Join Piazza for this course: http://piazza.com/mit/fall2014/18085
 Notes from Tuesday afternoon
 Notes from Tuesday afternoon 17 Sep
Course Topics:
 Applied Linear Algebra
 Applied Differential Equations
 Fourier Methods
 Algorithms
 Course outline
Additional Information:
 Quotes collected by my 18.085 class
 Goals for the Course: See applications of calculus, ODE, linear algebra, and discrete methods without going into too much proof.
 Textbook: Computational Science and Engineering (WellesleyCambridge, 2007).
 Grades: Homework 25%, 3 quizzes 75%, no final exam. Please email Prof. Strang about conflicts with quiz dates.
 Use of MATLAB for tedious calculations is encouraged, however you need to know how to do the basic algorithms taught in the course by hand (at least for small matrices) for the quizzes.
Class Resources:
 Movie of elimination: moe.m , realmmd.m
 Code to create K,T,B,C as sparse matrices
 MATLAB's backslash command to solve Ax = b: ps, pdf
 Getting started with Matlab: http://ocw.mit.edu/OcwWeb/Mathematics/1806Spring2005/RelatedResources/