Multivariate Calculus
&
Ordinary Differential Equations

Lecture 2



1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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?

  • $\ds \dif{P}{t} = kP $                                - 1st order
  • $ m\ddot{x} = -kx $                               - 2nd order
  • $ x(y'')^2+y'y'''+4y^5 = yy' $   - 3rd order
  • $ y' = xy $                                    - 1st order


1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.1.7 Motion of Projectiles

Click to drop apples! 🍏


1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.1.7 Motion of Projectiles

Click to throw apples! 🍏



1 Ordinary Differential Equations

1.1.7 Motion of Projectiles

Water park in Southbank, Brisbane.




1 Ordinary Differential Equations

1.1.8 Extra reading: Realistic models


👀 Extra reading! 📖

Not assessable! 😃







1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.1.10 Main points

  • You should understand what a differential equation (DE) is and how they can arise.
  • You should understand what it means for a function to be a solution to an ODE.
  • You should know what a general solution is.
  • You should know what an initial value problem (IVP) is.
  • You should be able to solve simple differential equations of the form $\ds {\dif{y}{x}=f(x)}$ or $\displaystyle \frac{\dup^2 y}{\dup x^2} = g(x)$.
  • You should understand what is meant by the order of a differential equation.
  • You should be able to derive, using Newton's second law, equations of motion in a force field (such as gravity).
  • You should understand the difference between an analytical and numerical solution to an ODE.


Credits