View source for Physics 212, 2018: Lectures 7

{{PHYS212-2018}}

Most models we will work on in this class will be written mathematically in terms of systems of coupled ordinary differential equations. How do we verify if our solution is not obviously wrong? One way of doing this is first to see (by solving the system numerically) whether the solutions settle down to a steady state, also known as a stable fixed point or a point attractor. If we can verify in an independent way that the attractor is correct, this will provide a good verification for our system. How do we do this?

Let's focus for now on a system of a single differential equation, like the constrained growth problem. If the system settled down, then in the equation <math>dx/dt=f(x,a)</math> the time derivative is zero, and the equation becomes <math>f(x,a)=0</math>. It's an algebraic equation! How do we find its roots? We use the [https://en.wikipedia.org/wiki/Newton%27s_method Newton-Raphson method], which I will derive in class.

Your own work: use the scripts I provided and change them to solve a quadratic equation. Compare the solution to the output of a function that gives an exact solution, which you wrote last time. Verify that different roots can be obtained using different initial conditions. Then use the scripts to solve an equation <math>e^x-ax=0</math> for arbitrary <math>a</math>

Please [https://docs.google.com/forms/d/e/1FAIpQLSeNvrPddCgvMfHdlSA8uYDddBFqLK81Rrr4l7jfozud6J7iww/viewform submit your work].