与えられた微分方程式 $(2xy \log y)dx + (x^2 + 3y^3)dy = 0$ に対して、$x^\alpha y^\beta$ 型の積分因子を見つけて、一般解を求める。

解析学微分方程式積分因子完全微分方程式
2025/5/19
## 問題7の解答

1. 問題の内容

与えられた微分方程式
(2xylogy)dx+(x2+3y3)dy=0(2xy \log y)dx + (x^2 + 3y^3)dy = 0
に対して、xαyβx^\alpha y^\beta 型の積分因子を見つけて、一般解を求める。

2. 解き方の手順

まず、与えられた微分方程式を
M(x,y)dx+N(x,y)dy=0M(x,y)dx + N(x,y)dy = 0
とすると、
M(x,y)=2xylogyM(x,y) = 2xy\log y
N(x,y)=x2+3y3N(x,y) = x^2 + 3y^3
である。
この微分方程式が完全微分形であるかどうか調べる。
My=2xlogy+2x\frac{\partial M}{\partial y} = 2x\log y + 2x
Nx=2x\frac{\partial N}{\partial x} = 2x
MyNx\frac{\partial M}{\partial y} \neq \frac{\partial N}{\partial x}
であるから、完全微分形ではない。
積分因子を xαyβx^\alpha y^\beta と仮定する。
xαyβM(x,y)dx+xαyβN(x,y)dy=0x^\alpha y^\beta M(x,y)dx + x^\alpha y^\beta N(x,y)dy = 0
すなわち
(xα+1yβ+1logy)dx+(xα+2yβ+3xαyβ+3)dy=0(x^{\alpha+1}y^{\beta+1}\log y)dx + (x^{\alpha+2}y^\beta + 3x^\alpha y^{\beta+3})dy = 0
が完全微分形となるように α\alphaβ\beta を定める。
(xα+1yβ+1logy)y=(β+1)xα+1yβlogy+xα+1yβ\frac{\partial (x^{\alpha+1}y^{\beta+1}\log y)}{\partial y} = (\beta+1)x^{\alpha+1}y^\beta \log y + x^{\alpha+1}y^\beta
(xα+2yβ+3xαyβ+3)x=(α+2)xα+1yβ\frac{\partial (x^{\alpha+2}y^\beta + 3x^\alpha y^{\beta+3})}{\partial x} = (\alpha+2)x^{\alpha+1}y^\beta
(xα+1yβ+1logy)y=(xα+2yβ+3xαyβ+3)x\frac{\partial (x^{\alpha+1}y^{\beta+1}\log y)}{\partial y} = \frac{\partial (x^{\alpha+2}y^\beta + 3x^\alpha y^{\beta+3})}{\partial x}
(β+1)xα+1yβlogy+xα+1yβ=(α+2)xα+1yβ(\beta+1)x^{\alpha+1}y^\beta \log y + x^{\alpha+1}y^\beta = (\alpha+2)x^{\alpha+1}y^\beta
(β+1)logy+1=α+2(\beta+1)\log y + 1 = \alpha+2
(β+1)logy=α+1(\beta+1)\log y = \alpha+1
この式が任意の yy に対して成り立つためには、β+1=0\beta+1=0 かつ α+1=0\alpha+1=0 でなければならない。
したがって β=1\beta=-1, α=1\alpha=-1 である。
よって、積分因子は x1y1=1xyx^{-1}y^{-1} = \frac{1}{xy} である。
積分因子をかけると
(2logy)dx+(xy+3y2x)dy=0(2\log y)dx + (\frac{x}{y} + \frac{3y^2}{x})dy = 0
積分因子をかけた後の微分方程式が完全微分形になっているかを確認する。
(2logy)y=2y\frac{\partial (2\log y)}{\partial y} = \frac{2}{y}
(xy+3y2x)x=1y3y2x2\frac{\partial (\frac{x}{y} + \frac{3y^2}{x})}{\partial x} = \frac{1}{y} - \frac{3y^2}{x^2}
積分因子が間違っていたようです。
もう一度計算し直します。
(β+1)xα+1yβlogy+xα+1yβ=(α+2)xα+1yβ(\beta+1)x^{\alpha+1}y^\beta \log y + x^{\alpha+1}y^\beta = (\alpha+2)x^{\alpha+1}y^\beta
(β+1)xα+1yβ(logy+1β+1)=(α+2)xα+1yβ(\beta+1)x^{\alpha+1}y^\beta (\log y + \frac{1}{\beta + 1}) = (\alpha+2)x^{\alpha+1}y^\beta
M(x,y)=2xylogyM(x,y) = 2xy\log y
N(x,y)=x2+3y3N(x,y) = x^2 + 3y^3
My=2xlogy+2x\frac{\partial M}{\partial y} = 2x\log y + 2x
Nx=2x\frac{\partial N}{\partial x} = 2x
MyNxN=2xlogy+2x2xx2+3y3=2xlogyx2+3y3\frac{\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}}{N} = \frac{2x\log y + 2x - 2x}{x^2+3y^3} = \frac{2x\log y}{x^2+3y^3}
NxMyM=2x2xlogy2x2xylogy=2xlogy2xylogy=1y\frac{\frac{\partial N}{\partial x} - \frac{\partial M}{\partial y}}{M} = \frac{2x - 2x\log y - 2x}{2xy \log y} = \frac{-2x\log y}{2xy \log y} = -\frac{1}{y}
μ(y)=e1ydy=elogy=1y\mu(y) = e^{\int -\frac{1}{y}dy} = e^{-\log y} = \frac{1}{y}
積分因子 μ(y)=1y\mu(y) = \frac{1}{y} をかけると
2xlogydx+(x2y+3y2)dy=02x\log y dx + (\frac{x^2}{y} + 3y^2) dy = 0
(2xlogy)y=2xy\frac{\partial (2x\log y)}{\partial y} = \frac{2x}{y}
(x2y+3y2)x=2xy\frac{\partial (\frac{x^2}{y} + 3y^2)}{\partial x} = \frac{2x}{y}
よって、完全微分形である。
F(x,y)=2xlogydx=x2logy+g(y)F(x,y) = \int 2x\log y dx = x^2 \log y + g(y)
Fy=x2y+g(y)=x2y+3y2\frac{\partial F}{\partial y} = \frac{x^2}{y} + g'(y) = \frac{x^2}{y} + 3y^2
g(y)=3y2g'(y) = 3y^2
g(y)=y3+Cg(y) = y^3 + C
F(x,y)=x2logy+y3=CF(x,y) = x^2 \log y + y^3 = C

3. 最終的な答え

x2logy+y3=Cx^2 \log y + y^3 = C

「解析学」の関連問題

与えられた関数の $n$ 次導関数を求めます。 (1) $x^m$ ($m < 0$, $n \le m$, $0 \le m < n$ に場合分けせよ) (2) $\frac{1}{1+x}$ (3...

微分導関数ガンマ関数高階微分
2025/6/7

実数列 $\{a_n\}_{n=1}^{\infty}$ を $a_n = (-1)^n a$ で定義する。ここで $a$ は正の実数である。この数列が収束しないことを示すための証明の空欄(ア〜エ)を...

数列収束極限証明
2025/6/7

実数列 $\{a_n\}_{n=1}^\infty$ を $a_n = (-1)^n a$ で定義する。この数列が収束しないことを示す証明における空欄ア~エを埋める問題。

数列収束証明実数列
2025/6/7

正の実数 $a$ に対して、数列 $\{a_n\}_{n=1}^{\infty}$ を $a_n = (-1)^n a$ で定める。この数列が収束しないことを示す証明の穴埋め問題。

数列収束証明絶対値
2025/6/7

正の実数 $a$ に対し、数列 $\{a_n\}_{n=1}^\infty$ を $a_n = (-1)^n a$ で定義する。このとき、数列 $\{a_n\}_{n=1}^\infty$ が収束しな...

数列収束極限絶対値三角不等式証明
2025/6/7

正の実数 $a$ に対して、数列 $\{a_n\}_{n=1}^\infty$ を $a_n = (-1)^n a$ で定義する。この数列が収束しないことを示すための証明の空欄ア~エに適切な語句を選択...

数列収束極限三角不等式
2025/6/7

問題は2つあります。 (1) 関数 $y = x^3 + ax^2 + 2ax + 4$ が極値を持つような $a$ の値の範囲を求める問題です。 (2) 次の接線の方程式を求める問題です。 ...

微分極値接線微分方程式
2025/6/7

曲線 $y = \frac{1}{x} + 1$、直線 $y = 0$、$x = 1$、$x = 3$ で囲まれた部分を $x$ 軸の周りに回転してできる立体の体積を求め、画像の空欄AからMを埋める。

積分回転体の体積定積分対数関数
2025/6/7

曲線 $y = x^2 - 1$ と $x$軸で囲まれた部分を $x$軸の周りに回転させてできる立体の体積を求める問題です。積分を用いて体積を計算する過程が示されており、空欄を埋める形式になっています...

積分回転体の体積定積分偶関数
2025/6/7

曲線 $y = \sqrt{x+2}$、x軸、y軸で囲まれた図形の面積を求める問題です。積分を用いて面積を計算します。

積分面積定積分置換積分関数のグラフ
2025/6/7