Integrating Factor Exactness Cooling
1. **Problem 11:** Given the differential equation \((2xy + y^2) dx + (2x^2 + 3xy) dy = 0\) and an integrating factor \(\mu(x,y) = x^m y^n\), find \(m\) and \(n\) so the equation becomes exact, then solve it.
2. **Step 1:** Check if the original equation is exact.
Calculate \(M = 2xy + y^2\) and \(N = 2x^2 + 3xy\).
Compute partial derivatives:
$$\frac{\partial M}{\partial y} = 2x + 2y$$
$$\frac{\partial N}{\partial x} = 4x + 3y$$
Since \(\frac{\partial M}{\partial y} \neq \frac{\partial N}{\partial x}\), the equation is not exact.
3. **Step 2:** Multiply by integrating factor \(\mu = x^m y^n\) to get exactness.
New functions:
$$\tilde{M} = x^m y^n (2xy + y^2) = 2x^{m+1} y^{n+1} + x^m y^{n+2}$$
$$\tilde{N} = x^m y^n (2x^2 + 3xy) = 2x^{m+2} y^n + 3x^{m+1} y^{n+1}$$
4. **Step 3:** For exactness, require:
$$\frac{\partial \tilde{M}}{\partial y} = \frac{\partial \tilde{N}}{\partial x}$$
Calculate:
$$\frac{\partial \tilde{M}}{\partial y} = 2x^{m+1} (n+1) y^n + x^m (n+2) y^{n+1}$$
$$\frac{\partial \tilde{N}}{\partial x} = 2(m+2) x^{m+1} y^n + 3(m+1) x^m y^{n+1}$$
Equate coefficients of \(x^{m+1} y^n\) and \(x^m y^{n+1}\):
$$2(n+1) = 2(m+2) \implies n+1 = m+2 \implies n = m+1$$
$$n+2 = 3(m+1) \implies n+2 = 3m + 3$$
Substitute \(n = m+1\) into second:
$$m+1 + 2 = 3m + 3 \implies m + 3 = 3m + 3 \implies 0 = 2m \implies m=0$$
Then \(n = m+1 = 1\).
5. **Step 4:** Integrating factor is \(\mu = x^0 y^1 = y\).
Multiply original equation by \(y\):
$$y(2xy + y^2) dx + y(2x^2 + 3xy) dy = (2xy^2 + y^3) dx + (2x^2 y + 3x y^2) dy = 0$$
Check exactness:
$$M = 2xy^2 + y^3, \quad N = 2x^2 y + 3x y^2$$
$$\frac{\partial M}{\partial y} = 4xy + 3y^2$$
$$\frac{\partial N}{\partial x} = 4xy + 3y^2$$
Exact.
6. **Step 5:** Find potential function \(\psi(x,y)\) such that:
$$\frac{\partial \psi}{\partial x} = M = 2xy^2 + y^3$$
Integrate w.r.t. \(x\):
$$\psi = \int (2xy^2 + y^3) dx = x^2 y^2 + x y^3 + h(y)$$
Differentiate \(\psi\) w.r.t. \(y\):
$$\frac{\partial \psi}{\partial y} = 2x^2 y + 3x y^2 + h'(y)$$
Set equal to \(N = 2x^2 y + 3x y^2\), so \(h'(y) = 0\).
Thus, \(h(y) = C\).
7. **Step 6:** General solution:
$$\psi(x,y) = x^2 y^2 + x y^3 = C$$
---
8. **Problem 12:** Given
$$\mu(x,y) = x f(y)$$
and the differential equation
$$(2 \sin^3 y + 3x \sin y) dx + (4x \cos y \sin^2 y + 2x^2 \cos y) dy = 0$$
Find the integrating factor and solve.
9. **Step 1:** Let \(\mu = x f(y)\). Multiply original equation by \(\mu\):
$$M = 2 \sin^3 y + 3x \sin y$$
$$N = 4x \cos y \sin^2 y + 2x^2 \cos y$$
New:
$$\tilde{M} = x f(y) (2 \sin^3 y + 3x \sin y) = 2x f(y) \sin^3 y + 3x^2 f(y) \sin y$$
$$\tilde{N} = x f(y) (4x \cos y \sin^2 y + 2x^2 \cos y) = 4x^2 f(y) \cos y \sin^2 y + 2x^3 f(y) \cos y$$
10. **Step 2:** For exactness:
$$\frac{\partial \tilde{M}}{\partial y} = \frac{\partial \tilde{N}}{\partial x}$$
Calculate partial derivatives:
$$\frac{\partial \tilde{M}}{\partial y} = 2x (f'(y) \sin^3 y + f(y) 3 \sin^2 y \cos y) + 3x^2 (f'(y) \sin y + f(y) \cos y)$$
$$\frac{\partial \tilde{N}}{\partial x} = 8x f(y) \cos y \sin^2 y + 6x^2 f(y) \cos y$$
11. **Step 3:** Equate coefficients of powers of \(x\):
For \(x\) terms:
$$2 (f'(y) \sin^3 y + 3 f(y) \sin^2 y \cos y) = 8 f(y) \cos y \sin^2 y$$
Simplify:
$$2 f'(y) \sin^3 y + 6 f(y) \sin^2 y \cos y = 8 f(y) \cos y \sin^2 y$$
Rearranged:
$$2 f'(y) \sin^3 y = 2 f(y) \cos y \sin^2 y$$
Divide both sides by \(2 \sin^2 y\):
$$f'(y) \sin y = f(y) \cos y$$
12. **Step 4:** Solve ODE for \(f(y)\):
$$\frac{f'(y)}{f(y)} = \frac{\cos y}{\sin y} = \cot y$$
Integrate:
$$\int \frac{f'(y)}{f(y)} dy = \int \cot y dy$$
$$\ln |f(y)| = \ln |\sin y| + C$$
$$f(y) = K \sin y$$
13. **Step 5:** Integrating factor:
$$\mu = x f(y) = K x \sin y$$
Ignore constant \(K\).
14. **Step 6:** Multiply original equation by \(x \sin y\) and solve exact equation similarly (omitted for brevity).
---
15. **Problem 13:** Cooling problem.
Given:
- Room temperature \(T_r = 18\) °C
- Initial temperature \(T_0 = 70\) °C
- Temperature after 5 min \(T(5) = 57\) °C
- Find time to reach \(T = 40\) °C
16. **Step 1:** Newton's law of cooling:
$$\frac{d\theta}{dt} = -k (\theta - T_r)$$
17. **Step 2:** Solve ODE:
Separate variables and integrate:
$$\int \frac{d\theta}{\theta - T_r} = -k \int dt$$
$$\ln |\theta - T_r| = -k t + C$$
$$\theta - T_r = C_1 e^{-k t}$$
At \(t=0\), \(\theta = 70\):
$$70 - 18 = C_1 = 52$$
So:
$$\theta(t) = 18 + 52 e^{-k t}$$
18. **Step 3:** Use \(\theta(5) = 57\):
$$57 = 18 + 52 e^{-5k}$$
$$39 = 52 e^{-5k}$$
$$e^{-5k} = \frac{39}{52} = \frac{3}{4}$$
$$-5k = \ln \frac{3}{4}$$
$$k = -\frac{1}{5} \ln \frac{3}{4} = \frac{1}{5} \ln \frac{4}{3}$$
19. **Step 4:** Find time \(t\) when \(\theta = 40\):
$$40 = 18 + 52 e^{-k t}$$
$$22 = 52 e^{-k t}$$
$$e^{-k t} = \frac{22}{52} = \frac{11}{26}$$
$$-k t = \ln \frac{11}{26}$$
$$t = -\frac{1}{k} \ln \frac{11}{26} = -\frac{1}{\frac{1}{5} \ln \frac{4}{3}} \ln \frac{11}{26} = -5 \frac{\ln \frac{11}{26}}{\ln \frac{4}{3}}$$
20. **Step 5:** Calculate time elapsed after 5 minutes:
$$\Delta t = t - 5 = -5 \frac{\ln \frac{11}{26}}{\ln \frac{4}{3}} - 5 = 5 \left(- \frac{\ln \frac{11}{26}}{\ln \frac{4}{3}} - 1\right)$$
This is the additional time needed to cool from 57 °C to 40 °C.
---
**Final answers:**
- Problem 11: \(m=0, n=1\), solution \(x^2 y^2 + x y^3 = C\)
- Problem 12: Integrating factor \(\mu = x \sin y\)
- Problem 13: Additional cooling time \(\Delta t = 5 \left(- \frac{\ln \frac{11}{26}}{\ln \frac{4}{3}} - 1\right)\) minutes