% SIMIODE-TeacherStudent-Template.tex and simiode.cls for LaTeX2e
% Version 22 November 2016.
% 
% This template is to be used with the  class file "simiode.cls". 
%  
% The latest version of this class and template files 
% can be found at
% 
%       <http://www.SIMIODE.org> in the Author section.
%
% A "Quick Start Guide" appears below.
%
% This file must be processed with version LaTeX2e or higher.
% Appreciation to Mark Frantz, IUPUI; John Thoo, Yuba College; and Eric Sullivan, Carroll College.
%
\documentclass{simiode}
%
% With this version of the tex file and the associated SIMIODE.cls file you can produce BOTH a Student Version and a Teacher Version of your modeling scenario or technique narrative.
%
% Student Version contains the title STUDENT VERSION along with the NAME OF SCENARIO, with no abstract, no keywords, no tags, no author identifiers, and no COMMENTS section. The Student Version contains only the STATEMENT of the modeling scenario only with appropriate figures, data, instructions, and perhaps a bibliography appropriate for student resourcing, but nothing in the bibliography which would "give away" the nature of the modeling scenario, etc.
%
%The Teacher Version contains the material in the Student Version, namely the title NAME OF SCENARIO, now headed with   TEACHER VERSION, but also full author  identification, abstract, key words, tag words, and, following the STATEMENT, a section of COMMENTS which can be discussion of presentation  ideas, "solutions," techniques, etc.
%
% Authors can prepare a full Teacher Version using the notions immediately below  to create a pdf for review, but in review please eliminate all information about the author as we use a double-blind referee system. Such information will be restored when accepted and published in SIMIODE at www.simiode.org.
%
%***********   Toggle the following line to turn the teacher version on (1) or off (0). ***********
%
\def\TeacherEdition{0} % Toggle this line's value.
\newcommand{\teacher}[1]{\ifnum\TeacherEdition=1 #1 \fi}
%
% The \teacher{ } command wraps any text that will appear only in the teacher version of
% the document. See this file for examples. This has been done for you in the document below.   
% Rename the pdf after you compile each final version.
%
\usepackage{amsmath,amssymb}
\usepackage{color}
\usepackage{graphicx}
%\usepackage{hyperref}
\usepackage{url}
%
\renewcommand{\baselinestretch}{1.3}
%
% NOTE regarding macros %%
%
% Any macros defined for your paper should be contained in 
% the top matter. Likewise, any environment definitions such 
% as \newtheorem or \newenvironment should also go in the 
% top matter. 
% 
% Do not \input your macros as separate files within your 
% final version because more files make it harder for the 
% editors and users to keep track of and/or modify and 
% append your paper.
%
% Top matter is above and the beginning of the document below.
%
\begin{document}

\newpage
%
% You only need to change the NAME OF SCENARIO in the next line and delete the footnote. 
% See above for instructions on toggling from Teacher Version to Student Version.
%
\title{\ifnum\TeacherEdition=1 TEACHER VERSION  \else STUDENT VERSION \fi\\Modeling with Differential Equations}


\markboth{Modeling with Differential Equations}{Modeling with Differential Equations}

\author{Michael Grayling\\
University of Cambridge\\
Cambridge  CB2 1TN  United Kingdom\\
michael.grayling@mrc-bsu.cam.ac.ukx
}


   
\makeStitlePDFLaTex


\begin{abstract}
Several models using first order differential equations are offered with some questions on formulating a differential equations model with solutions  provided.
\end{abstract}


\section*{SCENARIO DESCRIPTION}

The following is material reproduced from a larger mathematics project. All internal references to course work are retained.

The world around us poses many questions that can in theory be better understood using mathematics, and the process of writing down an equation describing how a variable of interest changes according to time or space, is called mathematical modelling. Unfortunately though, we are usually not able to write down such an equation straight away. But, we can often make simplifying assumptions that allow us to write down its rate of change with time: this gives us a differential equation. Solving this differential equation then gives us the solution to the original problem. For example; if we wished to describe how the number of bacteria, $y$, in a particular culture grew with time, $t$, proceeding to immediately write down a functional form for $y$'s dependence upon $t$, $y(t)$, is no easy feat. However, if we were to make the assumption that the number of bacteria grows at a rate proportional to its current size we can write down a differential equation for $y$:

$$\frac{dy}{dt} = \beta y\, ,$$

\noindent with $ \beta$ a constant of proportionality. We could then proceed to use integration in order to find the form for $y$ and we're done!

Now for STEP you need to be able to formulate your own simple differential equations, to then continue to solve. You should have gained some practice in this from your A-Level modules, however a good plan to keep in mind is the following four point procedure:

\begin{enumerate}
\item {\bf Identify the problems variables}: What are the dependent and independent variables we need to include?
\item {\bf Construct the differential equation}: What is the rate of change of the dependent variable given by the problem?
\item {\bf Solve the differential equation}: Use one of the standard methods, covered in more detail in the other articles in this module, to find the solution.
\item {\bf Interpret the mathematical solution}: Does the solution make sense in relation to the original problem? Are there any particularly interesting elements of the solution?
\end{enumerate}

One thing that the above does not mention though, that we should give greater consideration to is using initial or boundary conditions to find our exact solution of interest. Having integrated the differential equation, as usual a constant of integration will be introduced. You need to be able to interpret the original problem to write down the implied initial or boundary condition. Hopefully initial conditions will be familiar from A-Level, but boundary conditions may not. These imply a particular behaviour for a solution as the independent variable tends to some limit. For example if we integrated to find:

$$y(t)=e^{-t}+C\, ,$$
\noindent and our original problem implies $y\to 0$  as $t\to\infty$, then necessarily we must have $C=0$.

So the above covers the basics of formulating a differential equation, but to get good at this you need to gain much experience. Therefore, we'll now proceed to go through a few more examples before posing some problems for you to attempt!


\section*{Example 1}

 An object of mass $m=5$ kg is falling towards earth, with drag proportional to its velocity with a drag coefficient of 2 kg/sec. Formulate a differential equation for the velocity $v$.
 
\subsection*{Solution}
 Here we need a little bit of knowledge from mechanics, to known that we can write down a differential equation for $v$ using $F=m\cdot a=m\frac{dv}{dt}$. Then, all we need to do is realise the forcing pulling the object towards earth is going to be mg, and the drag slowing it down will be $2v$ using the information in the question. Putting this all together we have our differential equation:

$$m\frac{dv}{dt}=mg - 2v \, . $$

\section*{Example 2}

 Newton's Law of Cooling states that the rate of cooling of an object is proportional to the temperature difference between the room the object is in, and its own temperature. Write down a differential equation which would describe the temperature of the object as a function of time.

\subsection*{Solution}

 Clearly here our dependent variable will be the temperature of the object $T$, which depends on the independent variable $t$. Additionally here though we'll need a constant to describe the temperature of the room that we assume doesn't change with $t$. We'll call this $T_r$. Now we having everything we need to write down our differential equation using the information in the problem:

$$\frac{dT}{dt} =k(T-T_r)\, .$$

\section*{Exercise 1}

 A drug is administered to a patient as a constant rate $c$. As it is administered though it is converted by the patient’s body to other substances at a rate proportional to its current concentration. Formulate a differential equation for how the concentration of the drug, $D$, changes with time, $t$.


\section*{Exercise 2}
 The growth of a particular species of algae, $A$, in time, $t$, is assumed to be jointly proportional to both $\sqrt{A}$ and $10-A$. Write down the implied differential equation for $A$.

 
\section*{Summary}

OK, so that's the basics of mathematical modelling using differential equations! 


\teacher{
\section*{NOTES FOR TEACHERS}

The material in the STATEMENT above is reproduced from \cite{Grayling2014} as a part of a larger mathematics project. All internal references to course work are retained. We render the solutions offered in the original work here.

\subsection*{Exercise 1 Solution}


$\frac{dD}{dt} = c−kD$, $k$ a constant of proportionality.



\subsection*{Exercise 2 Solution}
$\frac{dA}{dt} =k\sqrt{A}(10−A)\, .$

\begin{thebibliography}{98}

\bibitem{Grayling2014} Grayling, Michael. 2014. First Order Differential Equations. University of Cambridge, Millenium Mathematics Project.  NRICH. \url{nrich.maths.org/11053}. 

\end{thebibliography}

}


\end{document}
