Ordinary Differential Equations Vocabulary

Overview of Differential Equations

  • Definition: A differential equation is an equation where the unknown is a function, and both the function and its derivatives may appear in the equation.

  • Differential equations are essential for physical theories, including:

    • Newton’s and Lagrange’s equations (Classical Mechanics).

    • Maxwell’s equations (Electromagnetism).

    • Schrödinger’s equation (Quantum Mechanics).

    • Einstein’s equation (General Relativity).

  • Ordinary Differential Equation (ODE): The unknown function depends on a single independent variable, usually represented as tt.

    • Examples:

      • Radioactive Decay: dudt(t)=ku(t)\frac{du}{dt}(t) = -k u(t), where k > 0 is a constant.

      • Newton's Law: md2xdt2(t)=f(t,x(t),dxdt(t))m \frac{d^2 x}{dt^2}(t) = f(t, x(t), \frac{dx}{dt}(t)).

  • Partial Differential Equation (PDE): The unknown function depends on two or more independent variables (e.g., t,x,y,zt, x, y, z).

    • Examples:

      • Heat Equation: Tt(t,x)=k(2Tx2+2Ty2+2Tz2)\frac{\partial T}{\partial t}(t, x) = k (\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2}).

      • Wave Equation: 2ut2=v2(2ux2+2uy2+2uz2)\frac{\partial^2 u}{\partial t^2} = v^2 (\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}).

  • Order of an Equation: The highest derivative order appearing in the equation.

First Order Linear Equations

  • Definition: A first order ODE is linear iff the source function ff is linear in its second argument:

    • y(t)=a(t)y+b(t)y'(t) = a(t)y + b(t).

  • Constant Coefficients: Linear iff both aa and bb are constants.

  • Variable Coefficients: Linear iff aa or bb depend on tt.

  • General Solution (Constant Coefficients): For y=ay+by' = ay + b with a0a \neq 0:

    • y(t)=ceatbay(t) = c e^{at} - \frac{b}{a}, where cRc \in \mathbb{R}.

  • Integrating Factor Method:

    • Transformation into a total derivative: μ(t)yaμ(t)y=μ(t)b(t)\mu(t)y' - a\mu(t)y = \mu(t)b(t).

    • The integrating factor μ(t)\mu(t) is chosen such that μ(t)=a(t)μ(t)\mu'(t) = -a(t)\mu(t), resulting in μ(t)=eA(t)\mu(t) = e^{-A(t)} where A(t)=a(t)dtA(t) = \int a(t) dt.

    • Final form: (μy)=μb(\mu y)' = \mu b.

  • Initial Value Problem (IVP): Consists of the ODE and an initial condition y(t0)=y0y(t_0) = y_0.

    • Unique Solution for Constant Coefficients: y(t)=(y0+ba)ea(tt0)bay(t) = (y_0 + \frac{b}{a}) e^{a(t-t_0)} - \frac{b}{a}.

Linear Variable Coefficient Equations and Bernoulli Equations

  • Theorem 1.2.1 (Variable Coefficients): If a,ba, b are continuous on (t1,t2)(t_1, t_2), then y=a(t)y+b(t)y' = a(t)y + b(t) has fundamental solution:

    • y(t)=ceA(t)+eA(t)eA(t)b(t)dty(t) = c e^{A(t)} + e^{A(t)} \int e^{-A(t)}b(t) dt

    • Where A(t)=a(t)dtA(t) = \int a(t) dt.

  • Bernoulli Equation: A first order nonlinear equation of the form:

    • y=p(t)y+q(t)yny' = p(t)y + q(t)y^n, where nRn \in \mathbb{R}.

    • Transformation: For n1n \neq 1, define v=1y(n1)v = \frac{1}{y^{(n-1)}}.

    • Linearized Equation for vv: v=(n1)p(t)v(n1)q(t)v' = -(n-1)p(t)v - (n-1)q(t).

Separable and Euler Homogeneous Equations

  • Separable Equations: Defined as h(y)y=g(t)h(y)y' = g(t).

    • Solution via direct integration: h(y)dy=g(t)dt+c\int h(y) dy = \int g(t) dt + c.

    • Implicit vs Explicit form: Solutions are often roots of algebraic equations (e.g., cubic polynomials).

  • Euler Homogeneous Equations: Equations of the form y=F(yt)y' = F(\frac{y}{t}).

    • Scale Invariance: F(cyct)=F(yt)F(\frac{cy}{ct}) = F(\frac{y}{t}).

    • Solving Strategy: Substitute v=ytv = \frac{y}{t}. This yields y=tvy = tv and y=v+tvy' = v + tv'.

    • Resulting Separable Equation: v=F(v)vtv' = \frac{F(v)-v}{t}, or vF(v)v=1t\frac{v'}{F(v)-v} = \frac{1}{t}.

Exact Differential Equations

  • Definition: An equation N(t,y)y+M(t,y)=0N(t, y)y' + M(t, y) = 0 is exact iff tN=yM\partial_t N = \partial_y M.

  • Potential Function (Poincaré Theorem): If exact, there exists ψ(t,y)\psi(t, y) such that:

    • yψ=N\partial_y \psi = N and tψ=M\partial_t \psi = M.

    • The differential equation behaves as dψdt(t,y(t))=0\frac{d\psi}{dt}(t, y(t)) = 0.

    • The solutions satisfy the level curve equation ψ(t,y(t))=c\psi(t, y(t)) = c.

  • Semi-Exact Equations: Equations that are not exact but become exact after multiplying by an integrating factor μ\mu.

    • If h=yMtNNh = \frac{\partial_y M - \partial_t N}{N} depends only on tt, then μ(t)=eH(t)\mu(t) = e^{H(t)} where H=hdtH = \int h dt.

    • If =yMxNM\ell = -\frac{\partial_y M - \partial_x N}{M} depends only on yy, then μ(y)=eL(y)\mu(y) = e^{L(y)} where L=dyL = \int \ell dy.

Applications of Linear Equations

  • Exponential Decay: dNdt=kN\frac{dN}{dt} = -kN.

    • Solution: N(t)=N0ektN(t) = N_0 e^{-kt}.

    • Half-life (\tau): The time for half the material to decay. Relationship: kτ=ln(2)k\tau = \ln(2).

  • Newton’s Cooling Law: The rate of change of temperature is proportional to the difference between the object and the medium temperature.

    • ΔT=kΔT\Delta T' = -k \Delta T, where ΔT(t)=T(t)Ts\Delta T(t) = T(t) - T_s.

    • Solution: T(t)=(T0Ts)ekt+TsT(t) = (T_0 - T_s) e^{-kt} + T_s.

  • Mixing Problems:

    • Volume balance: V(t)=riroV'(t) = r_i - r_o.

    • Salt balance: Q(t)=riqi(t)roQ(t)V(t)Q'(t) = r_i q_i(t) - r_o \frac{Q(t)}{V(t)}.

    • Model: Q(t)=a(t)Q(t)+b(t)Q'(t) = a(t)Q(t) + b(t).

Nonlinear Equations: Picard-Lindelöf and Picard Iteration

  • Picard-Lindelöf Theorem: Concerns the IVP y=f(t,y),y(t0)=y0y' = f(t, y), y(t_0) = y_0.

    • If ff is continuous and Lipschitz continuous in yy on a domain DaD_a, there exists a unique solution in a neighborhood of t0t_0.

  • Picard Iteration: Constructs a sequence of approximate solutions:

    • y0(t)=y0y_0(t) = y_0.

    • yn+1(t)=y0+t0tf(s,yn(s))dsy_{n+1}(t) = y_0 + \int_{t_0}^t f(s, y_n(s)) ds.

  • Existence vs Uniqueness: Non-uniqueness can occur if the Lipschitz condition is not met (e.g., y=y1/3y' = y^{1/3} has multiple solutions for y(0)=0y(0) = 0).

  • Direction Fields: Graphical representation of slope=f(t,y)\text{slope} = f(t, y), used for qualitative analysis.

Second Order Linear Equations

  • General Form: y+a1(t)y+a0(t)y=b(t)y'' + a_1(t) y' + a_0(t) y = b(t).

  • Linearity and Superposition: If y1,y2y_1, y_2 solve the homogeneous equation (b(t)=0b(t)=0), then c1y1+c2y2c_1 y_1 + c_2 y_2 also solves it.

  • Wronskian: W12(t)=y1y2y1y2W_{12}(t) = y_1 y'_2 - y'_1 y_2.

    • y1,y2y_1, y_2 are linearly independent iff W120W_{12} \neq 0.

  • Abel’s Theorem: The Wronskian of two solutions satisfies W+a1(t)W=0W' + a_1(t) W = 0, so W(t)=W(t0)ea1(s)dsW(t) = W(t_0) e^{-\int a_1(s) ds}.

  • Homogeneous Constant Coefficients: Characteristic polynomial p(r)=r2+a1r+a0=0p(r) = r^2 + a_1 r + a_0 = 0.

    • Roots r1,r2r_1, r_2 distinct real: y(t)=c1er1t+c2er2ty(t) = c_1 e^{r_1 t} + c_2 e^{r_2 t}.

    • Roots repeated real (r0r_0): y(t)=(c1+c2t)er0ty(t) = (c_1 + c_2 t) e^{r_0 t}.

    • Roots complex (α±iβ\alpha \pm i\beta): Real fundamental solutions are eαtcos(βt)e^{\alpha t}\cos(\beta t) and eαtsin(βt)e^{\alpha t}\sin(\beta t).

  • Reduction of Order: If y1y_1 is known, y2(t)=y1(t)eA1(t)y12(t)dty_2(t) = y_1(t) \int \frac{e^{-A_1(t)}}{y_1^2(t)} dt.

  • Euler Equidimensional Equation: (tt0)2y+a1(tt0)y+a0y=0(t - t_0)^2 y'' + a_1 (t - t_0) y' + a_0 y = 0.

    • Indicial polynomial: r(r1)+a1r+a0=0r(r-1) + a_1 r + a_0 = 0.

    • Solutions involve (tt0)r(t-t_0)^r or (tt0)rln(tt0)(t-t_0)^r \ln(t-t_0).

  • Nonhomogeneous Equations:

    • Undetermined Coefficients: Guess ypy_p based on the shape of the source term b(t)b(t).

    • Variation of Parameters: yp=u1y1+u2y2y_p = u_1 y_1 + u_2 y_2.

      • u1=y2fW12u'_1 = -\frac{y_2 f}{W_{12}}, u2=y1fW12u'_2 = \frac{y_1 f}{W_{12}}.

Power Series Solutions

  • Regular Points: Points where coefficients p(x),q(x)p(x), q(x) are analytic.

  • Method: Assume y(x)=n=0an(xx0)ny(x) = \sum_{n=0}^{\infty} a_n (x-x_0)^n. Plug into the ODE and solve the recurrence relation for coefficients ana_n.

  • Legendre Equation: (1x2)y2xy+(+1)y=0(1-x^2)y'' - 2xy' + \ell(\ell+1)y = 0. Solutions lead to Legendre Polynomials Pn(x)P_n(x).

  • Regular Singular Points: Point x0x_0 is a regular singular point iff (xx0)p(x)(x-x_0) p(x) and (xx0)2q(x)(x-x_0)^2 q(x) are analytic at x0x_0.

  • Frobenius Method: Assume y(x)=n=0an(xx0)n+ry(x) = \sum_{n=0}^{\infty} a_n (x-x_0)^{n+r}.

    • Define the indicial equation to find possible values for rr.

    • Bessel Equation: x2y+xy+(x2α2)y=0x^2 y'' + xy' + (x^2 - \alpha^2)y = 0. Solution defines Bessel functions Jα(x)J_\alpha(x).

Laplace Transform Method

  • Definition: L[f](s)=0estf(t)dtL[f](s) = \int_0^\infty e^{-st} f(t) dt.

  • Key Properties:

    • Linearity: L[af+bg]=aL[f]+bL[g]L[af + bg] = aL[f] + bL[g].

    • Translation (I): L[u(tc)f(tc)]=ecsL[f]L[u(t-c) f(t-c)] = e^{-cs} L[f].

    • Translation (II): L[ectf(t)]=L[f](sc)L[e^{ct} f(t)] = L[f](s-c).

    • Derivative to Multiplication: L[f]=sL[f]f(0)L[f'] = s L[f] - f(0).

    • Higher Derivatives: L[f]=s2L[f]sf(0)f(0)L[f''] = s^2 L[f] - s f(0) - f'(0).

  • Common Transforms:

    • L[1]=1sL[1] = \frac{1}{s}

    • L[eat]=1saL[e^{at}] = \frac{1}{s-a}

    • L[tn]=n!sn+1L[t^n] = \frac{n!}{s^{n+1}}

    • L[sin(at)]=as2+a2L[\sin(at)] = \frac{a}{s^2+a^2}, L[cos(at)]=ss2+a2L[\cos(at)] = \frac{s}{s^2+a^2}.

  • Dirac Delta Functions (Generalized Sources): δ(tc)\delta(t-c).

    • L[δ(tc)]=ecsL[\delta(t-c)] = e^{-cs} for c0c \geq 0.

    • Integral property: abf(t)δ(tc)dt=f(c)\int_a^b f(t) \delta(t-c) dt = f(c).

  • Convolution: (fg)(t)=0tf(τ)g(tτ)dτ(f*g)(t) = \int_0^t f(\tau) g(t-\tau) d\tau.

    • Theorem: L[fg]=L[f]L[g]L[f*g] = L[f]L[g].

Systems of Linear Differential Equations

  • Representation: x=A(t)x+b(t)\mathbf{x}' = A(t)\mathbf{x} + \mathbf{b}(t).

  • General Solution (Homogeneous): xgen(t)=cix(i)(t)\mathbf{x}_{gen}(t) = \sum c_i \mathbf{x}^{(i)}(t), where x(i)\mathbf{x}^{(i)} are fundamental solutions.

  • Matrix Exponential Solution: For constant matrix AA, IVP x=Ax,x(t0)=x0\mathbf{x}' = A\mathbf{x}, \mathbf{x}(t_0) = \mathbf{x}_0:

    • x(t)=eA(tt0)x0\mathbf{x}(t) = e^{A(t-t_0)} \mathbf{x}_0.

  • Diagonalizable Case: If A=PDP1A = PDP^{-1}, then eAt=PeDtP1e^{At} = P e^{Dt} P^{-1}. Solution decomposes into eigenmodes: x(t)=cieλitv(i)\mathbf{x}(t) = \sum c_i e^{\lambda_i t} \mathbf{v}^{(i)}.

  • Fundamental Matrix (X(t)X(t)): Matrix formed by fundamental solution vectors.

  • Nonhomogeneous Systems: x(t)=eA(tt0)x0+t0teA(tτ)b(τ)dτ\mathbf{x}(t) = e^{A(t-t_0)} \mathbf{x}_0 + \int_{t_0}^t e^{A(t-\tau)} \mathbf{b}(\tau) d\tau.

Autonomous Systems and Stability

  • Definition: System y=f(y)y' = f(y) where ff is time-independent.

  • Critical Points: xc\mathbf{x}_c such that f(xc)=0f(\mathbf{x}_c) = 0.

  • Stability Analysis:

    • Attractor (Sink): Solutions flow toward it.

    • Repeller (Source): Solutions flow away.

    • Linearization (Jackobian Matrix): Df(xc)D\mathbf{f}(\mathbf{x}_c). Stability governed by eigenvalues.

  • Phase Portraits (2D):

    • Nodes (real distinct eigenvalues, same sign).

    • SADDLE points (real distinct eigenvalues, opposite signs).

    • Spirals/Centers (complex eigenvalues).

  • Competition Species Model (Lotka-Volterra):

    • x1=x1(3x12x2)x'_1 = x_1(3 - x_1 - 2x_2), x2=x2(2x2x1)x'_2 = x_2(2 - x_2 - x_1).

    • Basins of attraction determine which species survives.

Boundary Value Problems and Heat Equation

  • Eigenfunction Problems: Find λ\lambda and non-zero yy for y+λy=0y'' + \lambda y = 0 with BC like y(0)=0,y(L)=0y(0)=0, y(L)=0.

    • Solutions: Eigenvalues λn=(nπL)2\lambda_n = (\frac{n\pi}{L})^2, Eigenfunctions yn(x)=sin(nπxL)y_n(x) = \sin(\frac{n\pi x}{L}).

  • Fourier Series: Expands f(x)f(x) on [L,L][-L, L] as a02+[ancos(nπxL)+bnsin(nπxL)]\frac{a_0}{2} + \sum [a_n \cos(\frac{n\pi x}{L}) + b_n \sin(\frac{n\pi x}{L})].

  • Heat Equation Application: tu=kx2u\partial_t u = k \partial^2_x u.

    • Separation of Variables: u(t,x)=v(t)w(x)u(t, x) = v(t) w(x).

    • Result for Dirichlet BC (u(t,0)=u(t,L)=0u(t,0)=u(t,L)=0): u(t,x)=cnek(nπL)2tsin(nπxL)u(t, x) = \sum c_n e^{-k(\frac{n\pi}{L})^2 t} \sin(\frac{n\pi x}{L}) where cnc_n are Fourier coefficients of initial data f(x)f(x).