On a cold day, it doesn't take long for a cup of coffee to become tepid; if our refrigerator breaks, our frozen items tend to warm up.
These day-to-day experiences result from heat moving about, a topic that greatly interested Sir Isaac Newton. Newton is famous for inventing calculus and rules of motion, but he also discovered a law for cooling:
An object's temperature changes at a rate proportional to the difference between it and the surrounding temperature.
This will be our starting point for a quick pass through the essential ideas of differential equations. No experience with physics is required; just a little single-variable calculus and a dash of vector arithmetic will do.
Let's cook up a differential equation based on Newton's Cooling Law:
“An object's temperature changes at a rate proportional to the difference between it and the surrounding temperature.”
If the surrounding air has constant temperature what equation can we infer for the object's temperature?
In all options, is a constant.
So Newton's Cooling Law leads us to the following differential equation: Right now we don't know anything about other than it's constant. For the wrong choice of our equation could predict something weird.
For what values of does the differential equation match our intuition about the real world?
Hint: Think about different cases for the initial temperature and what the equation implies about these scenarios.
We'll learn many techniques for solving differential equations, but it sometimes happens that an equation is just too difficult to solve.
That doesn't mean we give up and go home; we do what we can to glean as much information as possible from the equation itself.
To see how this works, let's start with a fact that's true for all equations in our course: solution curves never cross!
Taking this for granted, what can we say about from More than one option may be true!
But how can always increase when and not cross over The answer lies in the concavity, or "bendiness," of
From the figure above, we see that if the graph is concave up, while means it's concave down.
The second derivative is (see the solution for the details). What can we say about 's concavity from this?
Let's summarize what we've found without solving the equation:
If the function is strictly decreasing (increasing) and its graph is always concave up (down).
So the graph of depends on the initial condition which isn't surprising given what we uncovered in the last quiz.
Each of the four plots below shows two graphs that could be potential solution curves, one with initial condition above the other below.
Select the only option that displays correct graphs for The dashed line is in all plots.
There's another way of getting a quick sketch of the solution curves without finding explicitly.
Since also gives us the slope of lines tangent to solution curves, we lay down a grid on the -plane and sketch a short arrow with slope at every grid point
This results in a slope or direction field: starting at a point on the -axis, we follow the arrows to trace the shape of a solution with
The interactive below shows the direction field for our equation with adjustable and sliders. It also displays the solution.
Notice how the curve follows the arrows for all choices of
By following the direction field, we find an estimate for the solution: Here, measures how far we follow any arrow in the field.
The Python code below computes for and compares the results to the exact answer we'll find soon.
You can set the hop size the initial temperature “To”, and the air temperature “Ta”. We set
The larger the relative error, the worse the approximation; the smaller the relative error, the better! Play around with different values for the constants and then select the correct option.
from math import exp h = 10 # Hop size To = 0.5 # Initial temperature Ta= 1 # Air temperature n = 1 while n <= 10: difference = To - Ta exponential_h = exp(-n * abs(h) ) numerator_c = abs (Ta+(difference) * (1- (abs(h) ** n)) - (Ta+difference) * exponential_h) denominator_c= abs( Ta + (difference)* exponential_h) c = numerator_c / denominator_c print(c) n += 1 # Increase n by 1
Choose the size of your hop!
Set your initial temperature!
Set the ambient temperature!
This is our counter; please don't change it!
This measures the percent error of our approximation.
The method we just used is attributed to Euler, but beyond this brief shoutout, we won't have very much to say about numerical methods in differential equations. That's a course unto itself!
The main takeaway of the last problem is this: the smaller the value of the better the approximation tends to be!
The figure shows the first three approximations we get by following the direction field arrows (green) starting from a point (red) on a solution curve (blue) with small hop lengths.
The green points are pretty close to the curve, so this scheme works pretty well. In fact, this process of hopping along direction field arrows is one of the simplest ways of solving a differential equation with a computer, which can be programmed to do the tedious algebra.
To finish, let's find the actual solution for by rewriting the estimate as
According to the last page, this becomes exact if we take the limit as which is equivalent to taking
What is the result of this limit and therefore the exact value of
Good to know: Euler's number is defined by
This quiz gives the gist of differential equations, but the topic is vast, and we have a long road ahead of us.
While we aim to be as self-contained as possible, some parts of our course rely heavily on results and techniques in other areas.
Below is a guide to some of the more important supporting topics and where they pop up in our course:
The rules of differentiation (up to the chain rule) and integration (up to simple -substitutions) are used in all chapters: advanced integration methods are not needed at all.
Chapter 1 on first-order separable equations only uses single-variable calculus topics.
Chapter 2 on advanced first-order equations uses a few ideas from multivariable calculus and vector calculus, specifically the gradient, the curl, and level sets. These topics will recur in later chapters.
Chapters 3 and 4 on linear systems and higher-order equations rely on topics from linear algebra, including vector and matrix arithmetic and properties, eigenvalues, eigenvectors, and matrix diagonalizability.