Fall 2023: Differential Equations II (Math 332)

Overview

Lecture sketches and homework

To help you study for this class, I am posting lecture sketches. These are not a substitute for your personal lecture notes or coming to class (for instance, lots of details and motivation are not included in the sketches). I hope that they are useful to you for revisiting the material and for preparing for exams.

About the homework

• Homework problems are posted for each unit. Homework is submitted online, and you have an unlimited number of attempts. Only the best score is used for your grade.

  Most problems have a random component (which allows you to continue practicing throughout the semester without putting your scores at risk).

• Aim to complete the problems well before the posted due date.

  A 15% penalty applies if homework is submitted late.

• Collect a bonus point for each mathematical typo you find in the lecture notes (that is not yet fixed online), or by reporting mistakes in the homework system. Each bonus point is worth 1% towards a midterm exam.

  The homework system is written by myself in the hope that you find it beneficial. Please help make it as useful as possible by letting me know about any issues!

Exams and practice material

The following material will help you prepare for the exams.

• Midterm Exam 1:
  midterm01-practice.pdf, midterm01-practice-solution.pdf
  midterm01.pdf, midterm01-solution.pdf
• Midterm Exam 2:
  midterm02-practice.pdf, midterm02-practice-solution.pdf
  midterm02.pdf, midterm02-solution.pdf
• Final Exam:
  final-practice.pdf, final-practice-solution.pdf

Sage

For more involved calculations, we will explore the open-source free computer algebra system Sage.

If you just want to run a handful quick computations (without saving your work), you can use the text box below.

A convenient way to use Sage more seriously is https://cocalc.com. This free cloud service does not require you to install anything, and you can access your files and computations from any computer as long as you have internet. To do computations, once you are logged in and inside a project, you will need to create a "Sage notebook" as a new file.

Here are some other things to try:

• Sage can solve (symbolically!) simple initial value problems. For instance, to solve Example 22 in Lecture 4:

  y = function('y')(x)
  desolve(diff(y,x,2) - diff(y,x) - 2*y == 0, y, ics=[0,4,5])
  

  Unfortunately, this code currently only works for differential equations of first and second order. To likewise solve a third-order differential equation, we can add the instructions algorithm='fricas':

  y = function('y')(x)
  desolve(diff(y,x,3) == 3*diff(y,x,2) - 4*y, y, ics=[0,1,-2,3], algorithm='fricas')
  

  The output shows that the function $ y(x) = \frac13(3x - 2) e^{2x} + \frac53 e^{-x} $ is the unique solution to the initial value problem $ y''' = 3 y'' -4 y $ with $ y(0) = 1 $, $ y' (0) = -2 $ and $ y'' (0) = 3 $.
• Sage is happy to compute eigenvalues and eigenvectors. For instance, to solve Example 57 in Lecture 10:

  A = matrix([[8,-10],[5,-7]])
  A.eigenvectors_right()
  

• Sage can even compute general powers of a matrix. For instance, to solve Example 64 in Lecture 12:

  n = var('n')
  M = matrix([[0,1],[2,1]])
  M^n
  

• Similarly, Sage can compute matrix exponentials. For instance, to solve Example 72 in Lecture 13:

  M = matrix([[8,-10],[5,-7]])
  exp(M*x)
  

  (Note that x is predefined as a symbolic variable in Sage; that's why we don't need x = var('x') as in the case of n above.)
• We can plot phase portraits as follows:

  x,y = var('x y');
  streamline_plot((x*(y-1),y*(x-1)), (x,-3,3), (y,-3,3))
  

• Using matrix exponentials, we can easily solve inhomogeneous systems of differential equations. For instance, to solve Example 94 in Lecture 20:

  t = var('t')
  A = matrix([[2,3],[2,1]])
  y = exp(A*t)*integrate(exp(-A*t)*vector([0,-2*e^(3*t)]), t)
  y.simplify_full()
  

  Or, to solve Example 95 in Lecture 20:

  s, t = var('s, t')
  A = matrix([[1,2],[-1,4]])
  y = exp(A*t)*vector([1,2]) + exp(A*t)*integrate(exp(-A*s)*vector([0,2*e^s]), s,0,t)
  y.simplify_full()
  

• The following plots the solution to the heat equation with initial temperature profile equal to 1 and boundary conditions of zero on both ends.

  x, t = var('x,t')
  uxt = sum(4/(pi*n) * exp(-pi^2*n^2*t) * sin(pi*n*x) for n in range(1,20,2))
  density_plot(uxt, (x,0,1), (t,0,0.4), plot_points=200, cmap='hot')