Lecture 2
1.1.4 Initial value problems
An initial value problem (IVP) is the problem of solving an ODE subject to some initial value problem initial conditions of the form $y(t_0)=a,$ $y'(t_0)=b,$ etc.
The solution to an initial value problem no longer contains arbitrary constants from the general solution to the ODE - these are determined by the initial conditions of the problem at hand.
1.1.4 Initial value problems
Example: A flow-meter in a pipeline measures flow-through as $2+\sin t $ litres/second. How much fluid passes through the pipeline from time zero up to time $T$?
Flow through is Volume per unit of time.
$\ds\frac{dV}{dt}$ $=2 +\sin t$ $\Ra V(t) = 2t - \cos t + C.$
$V(0) = 2 (0) - \cos (0)+C $ $ = C-1 $ $ = 0$ $\,\Ra C= 1.$
So $V(T) = 2T - \cos T + 1.$
1.1.5 The order of an ODE
The order of the highest-order derivative in an ODE defines the order of the ODE.
Exercise: What is the order of each of these ODEs?
1.1.6 Riding your bike at constant speed
Find the position of your bike (at time $t$) if you are travelling at a constant speed of 60 km/h along a perfectly straight road.
If $x=x(t)$ is the distance you have travelled at time $t$ then the corresponding ODE is \[ \dif{x}{t}=60. \qquad \qquad (*) \]
We can directly integrate to find $x(t):$
$\ds\dif{x}{t}=60 \quad$ $\Ra\quad\ds \int \dup x=\int 60 \dup t.$
1.1.6 Riding your bike at constant speed
We can directly integrate to find $x(t):$
$\ds\dif{x}{t}=60 \quad$ $\Ra\quad\ds\int \dup x=\int 60 \dup t.$
This yields the general solution to the ODE \[ x(t)=60t+C. \] To determine the constant $C$ we need more information, such as your initial position at time $0$, in which case we would be solving an IVP.
1.1.6 Riding your bike at constant speed
For different values of $C$ we get different solutions, and below we graph some of these. If $C=0$ then $x=60t.$ All the other solutions are parallel to this line.
1.1.6 Riding your bike at constant speed
The curves $x=60t+C$ are called solution curves to the ODE $\ds\dif{x}{t}=60$.
Note that in this particular case all curves are straight lines with slope 60.
1.1.7 Motion of Projectiles
Example: Consider an apple falling under gravity. Find an expression for the height $y$ of the apple at time $t.$
Newton's 2nd law of motion is $F= ma.$
The velocity of the apple is $\dot{y}.$
The acceleration of the apple is $\ddot{y}.$
Force on the apple is $F = -mg.$
$\qquad \Ra \;-mg = m \ddot{y}$ $\quad \Ra\; \ddot{y} = -g.$
1.1.7 Motion of Projectiles
Example: Consider an apple falling under gravity. Find an expression for the height $y$ of the apple at time $t$.
👉 $\large \ddot{y} = -g$
Integrating we obtain $\dot{y} = -gt + C.$
Integrating again we get $y =\ds -\frac{g}{2}t^2+ Ct + D$
$\dot{y}(0) = -g(0)+ C $ $= C.$ ($C$ is the initial velocity)
$y(0) = \ds -\frac{g}{2}(0)^2+ C(0) + D$ $= D.$ (Position at $t=0$)
1.1.7 Motion of Projectiles
Click to drop apples! 🍏
1.1.7 Motion of Projectiles
Example: Now suppose you are throwing apples over your fence to your neighbour. Find an expression for the position $(x(t),y(t))$ of the apple if you assume the initial position is $x(0)= 0, y(0)=0$ and the initial velocity is $\dot{x}(0)=u,$ $\dot{y}(0)=v$.
We have two equations of motion:
$\ds\frac{\dup^2y}{\dup t^2}=-g\;$ and $\;\ds\frac{d^2x}{dt^2}=0.$
1.1.7 Motion of Projectiles
Example: Now suppose you are throwing apples over your fence to your neighbour. Find an expression for the position $(x(t),y(t))$ of the apple if you assume the initial position is $x(0)= 0, y(0)=0$ and the initial velocity is $\dot{x}(0)=u,$ $\dot{y}(0)=v$.
From previous example:
$y(t) = \ds - \frac{1}{2} gt^2+ Ct +D$ $= \ds - \frac{1}{2} gt^2+ vt,$
since $\,y(0)= 0\,$ and $\,\dot{y}(0) = v.$
1.1.7 Motion of Projectiles
Example: Now suppose you are throwing apples over your fence to your neighbour. Find an expression for the position $(x(t),y(t))$ of the apple if you assume the initial position is $x(0)= 0, y(0)=0$ and the initial velocity is $\dot{x}(0)=u$, $\dot{y}(0)=v$.
Now we need to solve: $\;\ds\frac{d^2x}{dt^2}=0.$
In this example gravity acts only in the $y$ direction.
1.1.7 Motion of Projectiles
Since $\ds\frac{d^2x}{dt^2}=0,$ we have that $\dot{x} = \alpha$ $= u$ ($\dot x(0)=u$)
$\Ra x(t) = ut + \beta. \,$ And $ x(0) = \beta $ $ = 0\,$ $\,\Ra x(t) = ut . $
Thus we have \[ \left(x(t), y(t)\right) = \left(ut, \ds - \frac{1}{2} gt^2+ vt\right) \]
1.1.7 Motion of Projectiles
Thus we have \[ \left(x(t), y(t)\right) = \left(ut, \ds - \frac{1}{2} gt^2+ vt\right) \]
Note: Set $t = x/u,$ then \[ y= - \frac{1}{2} g\left(\frac{x}{u}\right)^2+ v\left(\frac{x}{u}\right) \] \[ = - \frac{g}{2u^2} x^2+\frac{v}{u}x\quad \,\, \]
1.1.7 Motion of Projectiles
Click to throw apples! 🍏
1.1.7 Motion of Projectiles
Water park in Southbank, Brisbane.
1.1.8 Extra reading: Realistic models
👀 Extra reading! 📖
Not assessable! 😃
1.1.9 Analytical and numerical solutions
To solve an ODE (or IVP) analytically means to give a solution curve in terms of continuous functions defined over a specified interval, where the solution is obtained exactly by analytic means (e.g., by integration). The solution satisfies the ODE (and initial conditions) on direct substitution.
To solve an ODE (or IVP) numerically means to use an algorithm to generate a sequence of points which approximates a solution curve.
1.1.9 Analytical and numerical solutions
Important remark: As we have already seen, the ODE $$\dif{y}{t}=f(t)$$ can be solved analytically. The solution simply is \[ y=\int f(t) \dup t. \]
1.1.9 Analytical and numerical solutions
Important remark: As we have already seen, the ODE $\ds \dif{y}{t}=f(t)$ can be solved analytically. The solution simply is $ \ds y=\int f(t) \dup t. $
Now this may seem like a cop-out because all we are saying is that the solution is given by the anti-derivative of $f(t),$ and the above solution is as informative as the actual ODE. In practice one hopes to be able perform the above integral to get a more explicit form of the solution. Depending on the form of $f$ this may either be done exactly or by numerical means.
1.1.10 Main points