DifferentialEquations.Rmd

# Differential Equations

##  Modeling with Differential Equations

In mathematical modeling, we often describe real-world phenomena by focusing on how quantities change over time or space.  When a model involves an unknown function and one or more of its derivatives, it is called a **differential equation**.

This makes sense because many real systems are driven by change:

- populations grow or decline,  
- objects move and accelerate,  
- temperatures change,  
- springs oscillate,  
- chemicals react.  

Differential equations allow us to **predict future behavior** using information about **current change**.

---

### Population Growth Models


<div class="note">
**Exponential Growth Model**

Let:

- \(t\) = time (independent variable)  
- \(P(t)\) = population size (dependent variable)

>  Assumption:  The rate of change of the population is proportional to the population size.

This leads to the differential equation:
\[
\frac{dP}{dt} = kP
\]
where \(k\) is a constant.
</div>


**Interpretation:**

- If \(k > 0\), the population always increases.  
- As \(P\) increases, \(\frac{dP}{dt}\) increases → **faster growth over time**.  
- Growth rate depends on current population size.  


The solutions are exponential functions:
\[
P(t) = Ce^{kt}
\]

If \(P(0) = P_0\), then:
\[
C = P_0 \quad \Rightarrow \quad P(t) = P_0 e^{kt}.
\]
This is the exponential growth model.


::: {.example}
Exponential Growth

**Model:**
\[
\frac{dP}{dt} = 0.3P, \quad P(0)=100
\]

**Solution:**
\[
P(t) = 100e^{0.3t}
\]

**Interpretation:**  
Population grows exponentially from an initial size of 100.
:::


<div class="note">
**Logistic Growth Model (Limited Resources)**

Real environments have **limited resources**, so unlimited exponential growth is unrealistic.  To model this, we introduce the **carrying capacity** \(M\):
\[
\frac{dP}{dt} = kP\left(1 - \frac{P}{M}\right).
\]
This is the **logistic differential equation**.

**Behavior of the Logistic Model:**

- If \(P\) is small:  
  \(\frac{P}{M} \approx 0 \Rightarrow \frac{dP}{dt} \approx kP\) → exponential growth  
- If \(P > M\):    
  \(1 - \frac{P}{M} < 0 \Rightarrow \frac{dP}{dt} < 0\) → population decreases  
- If \(P \to M\):    
  \(\frac{dP}{dt} \to 0\) → population levels off  


**Equilibrium Solutions:**
\[
P(t) = 0, \qquad P(t) = M
\]

These are **equilibrium solutions**:

- \(P=0\): extinction equilibrium  
- \(P=M\): stable carrying capacity equilibrium  

Solutions move **away from 0** and **toward \(M\)**.
</div>

**Logistic Interpretation**

-  If \(P(0) < M\): population increases toward \(M\)    
-  If \(P(0) > M\): population decreases toward \(M\)  

**Conclusion:** Logistic growth stabilizes naturally.

---

### Physical Model: Motion of a Spring

Consider a mass \(m\) attached to a spring.

From physics:

- **Hooke’s Law:** restoring force \(= -kx\)  
- **Newton’s Second Law:** \(F = m a = m \frac{d^2x}{dt^2}\)  

This gives the differential equation:
\[
m\frac{d^2x}{dt^2} = -kx
\]

or
\[
\frac{d^2x}{dt^2} = -\frac{k}{m}x
\]


**Interpretation**

The acceleration is proportional to displacement but in the opposite direction.  This produces oscillatory motion.  Solutions are combinations of sine and cosine:
\[
x(t) = A\cos(\omega t) + B\sin(\omega t).
\]
This matches physical intuition: springs oscillate.



---

### General Differential Equations

A **differential equation**:

- Contains an unknown function  
- Contains one or more derivatives of that function  

**Order**

- **First-order:** highest derivative is first derivative  
- **Second-order:** highest derivative is second derivative  
- etc.



<div class="definition">
A function \(f(x)\) is a **solution** if substituting it and its derivatives into the equation makes the equation true.


The general solution is a **family of solutions** with constants:
\[
y = f(x) + C
\]

A differential equation **plus** an initial condition:
\[
y(t_0) = y_0.
\]
This selects **one unique solution** from the family.
</div>



::: {.example}
General Solution Check

Given:
\[
y = \frac{1+ce^t}{1-ce^t}
\]

This family satisfies:
\[
y' = \frac{1}{2}(y^2 - 1)
\]

Every value of \(c\) gives a valid solution → **general solution**.
:::



::: {.example}
Initial-Value Problem

Given:
\[
y' = \frac{1}{2}(y^2 - 1), \qquad y(0)=2
\]

From the general solution:
\[
y = \frac{1+ce^t}{1-ce^t}
\]

Apply the initial condition:
\[
2 = \frac{1+c}{1-c}
\Rightarrow c = \frac{1}{3}
\]

**Final solution:**
\[
y(t) = \frac{3+e^t}{3-e^t}
\]
:::

---

### Putting It All Together

Differential equations are the language of change.  They connect mathematics to reality by turning:

> “How something changes” into “What will happen next.”

They form the foundation of modeling in:

- biology (population, epidemics),  
- physics (motion, forces, energy),  
- chemistry (reaction rates),  
- economics (growth models),  
- engineering (systems and control),  
- data science and ML (dynamics and optimization).  

They are not just equations — they are predictive models of the real world.

#### Conceptual Takeaways

- Differential equations model **change-driven systems**  
- Models come from **physical laws** or **empirical assumptions**  
- Realistic models must include **constraints** (e.g. carrying capacity)  
- Differential equations describe **dynamics**, not just values  
- Many systems naturally lead to:  
  - exponential growth/decay    
  - logistic stabilization  
  - oscillatory motion  
- Initial conditions are essential for **prediction**  
- One equation can describe infinitely many behaviors  

---

#### Skills You Should Be Able to Do

After this section, you should be able to:

- Explain what a **differential equation** is  
- Identify the **order** of a differential equation  
- Interpret differential equations in real-world terms  
- Build simple models from assumptions  
- Recognize:  
  - exponential growth models  
  - logistic growth models  
  - oscillation models  
- Understand equilibrium solutions  
- Distinguish between:  
  - general solutions  
  - particular solutions  
- Solve simple **initial-value problems**  
- Interpret solutions physically and graphically  
- Explain why differential equations are used in modeling  



---

####  Problems



















##  Direction Fields and Euler’s Method


In most real-world applications, differential equations cannot be solved explicitly in closed form. This means we often cannot write solutions as formulas.  However, we can still understand solution behavior using:

- **Graphical methods** → direction (slope) fields   
- **Numerical methods** → Euler’s method  

These approaches allow us to:

- Visualize solution curves  
- Predict long-term behavior  
- Approximate solutions at specific points  
- Interpret physical and biological models  

---

### Direction Fields (Slope Fields)


A first-order differential equation has the form
\[
y' = F(x,y).
\]

At each point \((x,y)\), the value \(F(x,y)\) gives the **slope** of the solution curve at that point.  A **direction field** is a collection of short line segments drawn at many points \((x,y)\), each with slope \(F(x,y)\).  These line segments show the *direction* in which solution curves must flow.  This allows us to sketch solution curves without solving the equation.


::: {.example}
Basic Direction Field

\[
y' = x + y, \qquad y(0) = 1
\]


At any point \((x,y)\), the slope is:
\[
\text{slope} = x + y.
\]

So:

- At \((0,1)\): slope = \(1\)  
- At \((1,2)\): slope = \(3\)  
- At \((-1,1)\): slope = \(0\)  

**Solution Behavior**

1. Draw short line segments with slope \(x+y\) at many points.  
2. Start at the initial point \((0,1)\).  
3. Trace a curve that stays **parallel to the local line segments**.  

This produces a smooth solution curve following the direction field.  See the textbook.
:::


---

### Autonomous Differential Equations

An equation of the form
\[
y' = f(y)
\]
is called **autonomous** (no explicit \(x\) dependence).

**Key Properties:**

- Slopes depend only on \(y\)  
- Direction field has **horizontal translation symmetry**  
- Solutions can be shifted left/right in \(x\)  
- Constant solutions \(y = C\) where \(f(C)=0\) are called equilibrium solutions  


::: {.example}
Equilibrium Solutions


\[
\frac{dI}{dt} = 15 - 3I
\]

**Equilibrium Condition**
\[
0 = 15 - 3I \Rightarrow I = 5
\]

So:
\[
I(t) = 5
\]
is an equilibrium solution.

**Interpretation**

All nearby solutions move toward \(I=5\).  This represents stable equilibrium behavior.
:::

---

### Euler’s Method


Euler’s method is a numerical algorithm for approximating solutions of:
\[
y' = F(x,y), \qquad y(x_0)=y_0
\]
Instead of solving exactly, we approximate step-by-step using tangent lines.


<div class="note">
**Euler’s Method Formula**

Let step size be \(h\). Define points:
\[
x_n = x_0 + nh
\]
Then the recursion formula is:
\[
y_{n} = y_{n-1} + h\,F(x_{n-1}, y_{n-1})
\]
</div>

This means:

> New value = old value + (step size × slope)


::: {.example}
Euler Approximation

Approximate $y$ near $x = 0.5$ given
\[
y' = x + y, \qquad y(0)=1, \qquad h=0.1
\]

So:
\[
F(x,y) = x+y.
\]


**Step 1:**
\[
y_1 = 1 + 0.1(0+1) = 1.1
\]

**Step 2:**
\[
y_2 = 1.1 + 0.1(0.1+1.1) = 1.22
\]

**Step 3:**
\[
y_3 = 1.22 + 0.1(0.2+1.22) = 1.362
\]

**Table**

| \(n\) | \(x_n\) | \(y_n\) |
|------|--------|--------|
| 0 | 0.0 | 1.0000 |
| 1 | 0.1 | 1.1000 |
| 2 | 0.2 | 1.2200 |
| 3 | 0.3 | 1.3620 |
| 4 | 0.4 | 1.5282 |
| 5 | 0.5 | 1.7210 |

So,
\[
y(0.5) \approx 1.721
\]
:::



**Accuracy and Step Size**

- Smaller \(h\) → better accuracy  
- Larger \(h\) → faster computation, lower accuracy  
- As \(h \to 0\), Euler’s method approaches the true solution

Euler’s method is the foundation of modern numerical solvers used in scientific computing.

---

### Putting It All Together

> Differential equations describe change.  Direction fields show **how change behaves**.  Euler’s method shows **how change accumulates**.  

Together, they provide a complete framework for understanding systems when exact solutions are impossible.

#### Conceptual Takeaways

- Most differential equations **cannot be solved explicitly**  
- Direction fields show **global behavior** of solutions  
- Euler’s method gives **local numerical approximations**  
- Direction fields → qualitative understanding    
- Euler’s method → quantitative approximation  
- Equilibrium solutions model **steady states**  
- Autonomous equations exhibit **horizontal shift symmetry**  

---

#### Skills You Should Be Able to Do

After this section, you should be able to:

**Direction Fields**

- Interpret \(y' = F(x,y)\) as a slope rule  
- Construct and read direction fields  
- Sketch solution curves using slope fields  
- Identify equilibrium solutions  
- Predict long-term behavior from direction fields  

**Euler’s Method**

- Apply the Euler update formula  
- Build numerical solution tables  
- Approximate \(y(x)\) values  
- Understand the effect of step size \(h\)  
- Interpret numerical solutions physically  



---

####  Problems



















##  Separable Equations

So far, we have studied differential equations using:

- **Geometric methods** → direction fields  
- **Numerical methods** → Euler’s method  

These approaches help us visualize and approximate solutions. However, in some special cases, we can go further and find explicit symbolic formulas for solutions.  One of the most important solvable classes is the class of **separable differential equations**.


<div class="definition">
A first-order differential equation is called **separable** if it can be written in the form
\[
\frac{dy}{dx} = t(x)\,f(y)
\]
</div>

That is, the right-hand side of a separable DE factors into:

- a function of \(x\) only, and  
- a function of \(y\) only.  


If the equation is separable, we can rewrite it as
\[
h(y)\,dy = t(x)\,dx
\]

so that:

- all \(y\)-terms are on one side,  
- all \(x\)-terms are on the other side.  

Then we integrate:
\[
\int h(y)\,dy = \int t(x)\,dx
\]
This produces an **implicit solution**.  Sometimes we can solve explicitly for \(y\); sometimes we cannot.


<div class="definition">
**Solving a Separable Equation**

**Step 1:** Write the equation in differential form  
\[
\frac{dy}{dx} = t(x)f(y)
\]

**Step 2:** Separate variables  
\[
\frac{1}{f(y)}\,dy = t(x)\,dx
\]

**Step 3:** Integrate both sides  
\[
\int \frac{1}{f(y)}\,dy = \int t(x)\,dx
\]

**Step 4:** Add constant of integration \(C\)

**Step 5:** (If possible) solve explicitly for \(y\)
</div>



::: {.example}
Basic Separable Equation

Solve
\[
\frac{dy}{dx} = \frac{x^2}{y^2}, \qquad y(0)=2
\]

**Solution:**

**Step 1: Separate variables**
\[
y^2\,dy = x^2\,dx
\]

**Step 2: Integrate**
\[
\int y^2\,dy = \int x^2\,dx
\]
\[
\frac{1}{3}y^3 = \frac{1}{3}x^3 + C
\]

**Step 3: Solve for \(y\)**
\[
y^3 = x^3 + K
\]
\[
y = \sqrt[3]{x^3 + K}
\]

**Step 4: Apply initial condition \(y(0)=2\)**
\[
2 = \sqrt[3]{K} \Rightarrow K=8
\]

**Final Solution**
\[
y = \sqrt[3]{x^3 + 8}
\]
:::



::: {.example}
Implicit Solution

Solve
\[
\frac{dy}{dx} = \frac{6x^2}{2y+\cos y}
\]

**Solution:**

**Separation**
\[
(2y+\cos y)\,dy = 6x^2\,dx
\]

**Integration**
\[
\int (2y+\cos y)\,dy = \int 6x^2\,dx
\]
\[
y^2 + \sin y = 2x^3 + C
\]

**Result (Implicit Solution)**
\[
y^2 + \sin y = 2x^3 + C
\]

> This equation cannot be solved explicitly for \(y\), but it still defines the solution curve.

:::



::: {.example}
Exponential Growth/Decay Structure

Solve
\[
y' = x^2 y
\]

**Solution:**

**Separation**
\[
\frac{1}{y}\,dy = x^2\,dx
\]

**Integration**
\[
\ln|y| = \frac{x^3}{3} + C
\]

**Solve for \(y\)**
\[
y = Ae^{x^3/3}
\]

**General Solution**
\[
y = Ae^{x^3/3}
\]
:::



---

### Physical Application Example


::: {.example}
Electric Circuit Model

Solve
\[
\frac{dI}{dt} = 15 - 3I, \qquad I(0)=0
\]

**Solution:**

**Separation**
\[
\frac{dI}{15-3I} = dt
\]

**Integration**
\[
-\frac{1}{3}\ln|15-3I| = t + C
\]

**Solve for \(I(t)\)**
\[
I(t) = 5 - 5e^{-3t}
\]

**Limiting Value**
\[
\lim_{t\to\infty} I(t) = 5
\]

**Interpretation**
The current increases over time and approaches a **steady-state equilibrium** of **5 amperes**.
:::


---

### Putting It All Together


> Direction fields show **how solutions behave**.  Euler’s method shows **how solutions accumulate**.  Separable equations show **why solutions exist symbolically**

Together, these form the three core perspectives of differential equations: geometric, numerical, and analytic.

#### Conceptual Takeaways

- Some differential equations can be solved **symbolically**  
- Separable equations allow **exact analytic solutions**  
- The key structure is:  
  \[
  y' = t(x)f(y)
  \]  
- Separation transforms differential equations into **integral equations**  
- Solutions may be:  
  - explicit (\(y=f(x)\)), or    
  - implicit (\(F(x,y)=C\))  
- Equilibrium values arise naturally in separable models  
- Many physical systems are separable:  
  - population models   
  - cooling laws  
  - circuits  
  - mixing problems  
  - drug concentration models  

---

#### Skills You Should Be Able To Do

After this section, you should be able to:

- Recognize when a differential equation is separable  
- Factor equations into \(t(x)f(y)\) form  
- Separate variables correctly  
- Integrate both sides  
- Apply constants of integration properly  
- Solve for \(y\) when possible  
- Work with implicit solutions  
- Solve initial-value problems  
- Interpret equilibrium behavior   
- Model real systems symbolically  
- Connect direction fields to analytic solutions  
- Understand limiting behavior (\(t \to \infty\))  


---

####  Problems



















##  Models for Population Growth

Population growth is one of the most important real-world applications of differential equations. In this section, we develop mathematical models that describe how populations change over time, starting with simple exponential growth and progressing to more realistic models that incorporate environmental limits and resource constraints.

We focus on two core models:

- **The Law of Natural Growth (Exponential Growth Model)**  
- **The Logistic Growth Model**  

These models help us understand how assumptions about resources, space, and environment shape long-term population behavior.

---

### The Law of Natural Growth (Exponential Model)


The law of natural growth assumes that:

> The rate of change of the population is proportional to the size of the population.

Mathematically, this is written as:

\[
\frac{dP}{dt} = kP
\]

where:

- \( P(t) \) = population size at time \( t \)  
- \( k \) = growth constant   
  - \( k > 0 \): population increases  
  - \( k < 0 \): population decreases  

This means larger populations grow faster, and smaller populations grow more slowly, but the *relative growth rate* remains constant.


This differential equation is separable and solves to:
\[
P(t) = P_0 e^{kt}
\]
where:

- \( P_0 = P(0) \) is the initial population.  

This model produces **exponential growth** (or decay).


**Interpretation**

\[
\frac{1}{P}\frac{dP}{dt} = k
\]

This shows that the relative growth rate (growth per individual) is constant.  So exponential growth occurs when each individual contributes equally and independently to growth, and resources are unlimited.


::: {.example}
Exponential Growth

A bacterial population grows according to:
\[
\frac{dP}{dt} = 0.3P, \qquad P(0) = 500.
\]
Find the population function and the population after 5 hours.

**Solution**
\[
P(t) = 500e^{0.3t}
\]

At \( t = 5 \):

\[
P(5) = 500e^{1.5} \approx 2240
\]
:::



Exponential growth assumes:

- unlimited resources  
- unlimited space  
- no competition  
- no environmental constraints  

This makes it realistic only for short-term growth or early population stages.



---

### The Logistic Growth Model


Real populations do not grow exponentially forever. As populations increase:

- resources become limited  
- competition increases  
- growth slows  
- populations stabilize  

To model this, we introduce a **carrying capacity**.


\[
\frac{dP}{dt} = kP\left(1 - \frac{P}{M}\right)
\]

where:

- \( M \) = carrying capacity (maximum sustainable population)  
- \( k \) = intrinsic growth rate  


**Interpretation:**

- If \( P \ll M \):  
  \[
  \frac{dP}{dt} \approx kP
  \]  
  → growth is approximately exponential.  

- If \( P \to M \):  
  \[
  \frac{dP}{dt} \to 0
  \]  
  → growth slows and stabilizes.  

- If \( P > M \):  
  \[
  \frac{dP}{dt} < 0
  \]  
  → population decreases.  


The logistic equation has two equilibrium populations:
\[
P = 0 \quad \text{and} \quad P = M
\]

- \( P = 0 \): unstable equilibrium  
- \( P = M \): stable equilibrium (carrying capacity)  


The solution to the logistic differential equation is:
\[
P(t) = \frac{M}{1 + Ae^{-kt}},
\]
where:
\[
A = \frac{M - P_0}{P_0}
\]
and \( P_0 = P(0) \).


::: {.example}
Logistic Growth

Solve
\[
\frac{dP}{dt} = 0.08P\left(1 - \frac{P}{1000}\right), \qquad P(0) = 100
\]

**Solution**

\[
A = \frac{1000 - 100}{100} = 9
\]

\[
P(t) = \frac{1000}{1 + 9e^{-0.08t}}
\]

:::




The logistic model has maximum growth at:
\[
P = \frac{M}{2}.
\]
This is where the solution curve changes concavity and the population grows fastest.


---

### Comparing the Models

| Feature | Exponential Model | Logistic Model |
|------|------|------|
| Resources | Unlimited | Limited |
| Long-term behavior | Unbounded growth | Stabilizes |
| Realism | Short-term only | Long-term realistic |
| Equilibrium | None | Carrying capacity |
| Shape | Exponential | S-shaped (sigmoid) |


Extensions of the logistic model include:

- **Harvesting models**:
  \[
  \frac{dP}{dt} = kP\left(1 - \frac{P}{M}\right) - c
  \]  
- **Threshold extinction models**:
  \[
  \frac{dP}{dt} = kP\left(1 - \frac{P}{M}\right)\left(1 - \frac{m}{P}\right)
  \]  

These models account for:

- constant removal (fishing, hunting)  
- extinction thresholds  
- minimum viable populations  

---

### Putting It All Together

Population growth models connect calculus to real biological systems by translating assumptions into differential equations:

- **Exponential growth** models unlimited environments and constant relative growth rates.  
- **Logistic growth** introduces environmental limits through carrying capacity.  
- Direction fields reveal long-term behavior visually.  
- Equilibrium solutions describe stable and unstable population states.  
- Inflection points identify where growth is fastest.  
- Real populations typically follow **logistic-like behavior**, not exponential growth.  

These models demonstrate how simple equations can capture complex real-world dynamics when built on meaningful assumptions.

---

#### Conceptual Takeaways

- Growth models depend on assumptions about resources and space  
- Differential equations encode biological laws  
- Long-term behavior is more important than short-term predictions  
- Equilibria represent natural population limits  
- Carrying capacity defines sustainability  
- Mathematical structure reflects biological reality  

---

#### Skills You Should Be Able to Do

After this section, you should be able to:

- Interpret population growth differential equations  
- Solve separable differential equations  
- Identify exponential vs logistic models  
- Compute carrying capacity from equations  
- Find equilibrium solutions  
- Analyze stability of equilibria  
- Interpret direction fields  
- Identify inflection points  
- Model real data using growth equations  
- Compare model predictions to observations  
- Explain biological meaning of parameters  
- Choose appropriate models based on context  


---

####  Problems



















##  Linear Equations

First-order linear differential equations form one of the most important classes of differential equations in calculus. They appear naturally in physics, biology, economics, engineering, and chemistry, especially in models involving growth, decay, circuits, mixing, and population dynamics.


<div class="definition">
A **first-order linear differential equation** is any equation that can be written in the form
\[
\frac{dy}{dx} + P(x)y = Q(x),
\]
where:

- \( P(x) \) and \( Q(x) \) are continuous functions on an interval,  
- \( y \) is the unknown function,  
- the equation is **linear in \( y \)** (no powers of \( y \), no products like \( y^2 \), \( y\sin y \), etc.).  
</div>



The equation is called **linear** because:

- \( y \) appears only to the first power  
- \( y \) is not inside nonlinear functions  
- coefficients depend only on \( x \)  

This makes the structure predictable and solvable using a general method.



Many real systems naturally produce linear equations:

- population harvesting models  
- electric circuits  
- cooling/heating laws  
- mixing problems  
- growth with forcing terms  
- economic adjustment models  

They are often **not separable**, but still solvable using a systematic technique.

---

### The Integrating Factor Method

To solve
\[
\frac{dy}{dx} + P(x)y = Q(x),
\]
we use an **integrating factor** \( I(x) \), defined by:
\[
I(x) = e^{\int P(x)\,dx}
\]

Multiplying the equation by \( I(x) \) transforms the left-hand side into the derivative of a product:

\[
I(x)\frac{dy}{dx} + I(x)P(x)y 
= \frac{d}{dx}\big(I(x)y\big)
\]

So the equation becomes:

\[
\frac{d}{dx}\big(I(x)y\big) = I(x)Q(x)
\]

If we integrate both sides:

\[
I(x)y = \int I(x)Q(x)\,dx + C
\]

and solve for \( y \):

\[
y(x) = \frac{1}{I(x)}\left(\int I(x)Q(x)\,dx + C\right)
\]


<div class="note">
To solve a first-order linear DE:

1. Put in standard form:
   \[
   \frac{dy}{dx} + P(x)y = Q(x)
   \]

2. Compute the integrating factor:
   \[
   I(x) = e^{\int P(x)\,dx}
   \]

3. Multiply the entire equation by \( I(x) \)

4. Rewrite left side as a product derivative:
   \[
   \frac{d}{dx}(I(x)y)
   \]

5. Integrate both sides

6. Solve for \( y(x) \)
</div>



::: {.example}
Basic Linear Equation

Solve
\[
\frac{dy}{dx} + 3x^2 y = 6x^2
\]

**Solution**

Standard form is already given:

- \( P(x) = 3x^2 \)  
- \( Q(x) = 6x^2 \)  

Integrating factor:
\[
I(x) = e^{\int 3x^2 dx} = e^{x^3}
\]

Multiply through:
\[
e^{x^3}\frac{dy}{dx} + 3x^2 e^{x^3} y
= 6x^2 e^{x^3}
\]

Left-hand side becomes:
\[
\frac{d}{dx}(e^{x^3}y) = 6x^2 e^{x^3}
\]

Integrate:
\[
e^{x^3}y = \int 6x^2 e^{x^3} dx
\]

Let \( u = x^3 \), \( du = 3x^2 dx \), so:
\[
e^{x^3}y = 2e^{x^3} + C
\]

Solve:
\[
y = 2 + Ce^{-x^3}
\]
:::



::: {.example}
Initial Value Problem

Solve
\[
x^2 y' + xy = 1, \qquad y(1) = 2, \quad x>0
\]

**Solution**

Divide by \( x^2 \):
\[
y' + \frac{1}{x}y = \frac{1}{x^2}
\]

Integrating factor:
\[
I(x) = e^{\int \frac{1}{x}dx} = e^{\ln x} = x
\]

Multiply through:
\[
x y' + y = \frac{1}{x}
\]

\[
\frac{d}{dx}(xy) = \frac{1}{x}
\]

Integrate:
\[
xy = \ln x + C
\]

Solve:
\[
y = \frac{\ln x + C}{x}
\]

Apply initial condition \( y(1)=2 \):
\[
2 = \frac{0 + C}{1} \Rightarrow C = 2
\]

Therefore,
\[
y = \frac{\ln x + 2}{x}
\]
:::

---

### Application Example: Electric Circuits (RL Circuit)

A simple electric circuit with:

- resistance \( R \),  
- inductance \( L \),  
- voltage source \( E(t) \)  

satisfies:
\[
L\frac{dI}{dt} + RI = E(t)
\]
This is a first-order linear differential equation in \( I(t) \).

Solution structure:
\[
I(t) = \frac{E}{R}\big(1 - e^{-Rt/L}\big)
\]

This shows:

- current starts at 0   
- increases smoothly  
- approaches a limiting value \( \frac{E}{R} \)  

---

### Putting It All Together

First-order linear differential equations form a powerful modeling framework:

- Their structure guarantees solvability.  
- The integrating factor transforms structure into exact derivatives.  
- Solutions naturally separate into:  
  - **transient parts** (decaying exponentials)  
  - **steady-state parts** (long-term behavior)  
- Long-term system behavior is often independent of initial conditions.  
- They unify many applied models into one mathematical form.  

This reveals a deep idea in calculus:

> Structure determines solvability, and assumptions determine behavior.

---

#### Conceptual Takeaways

- Linearity refers to structure, not graph shape  
- Integrating factors convert equations into exact derivatives  
- Exponential terms control stability and decay  
- Long-term behavior dominates system dynamics  
- Transient behavior fades over time  
- Steady-state solutions govern equilibrium  
- Mathematical form encodes physical laws  
- One method solves many real-world systems  

---

#### Skills You Should Be Able To Do

After this section, you should be able to:

- Identify linear vs nonlinear differential equations  
- Convert equations to standard form  
- Compute integrating factors  
- Apply the integrating factor method  
- Solve initial value problems  
- Interpret physical meaning of solutions  
- Analyze long-term behavior  
- Identify transient vs steady-state components  
- Solve applied modeling problems  
- Work with forcing functions  
- Solve circuit models  
- Recognize when equations are linear but not separable  
- Connect structure to solution behavior  


---

####  Problems



















##  Predator-Prey Systems