기본 콘텐츠로 건너뛰기

라벨이 Separable Equations인 게시물 표시

[matplotlib]quiver()함수

비선형 1차 미분방정식(Separable Equations)

비선형 1차 미분방정식(Separable Equations) 비선형 1차 미분방정식의 첫번째 유형은 분리가능한 형태(Separable form) 입니다. $$N(y)\frac{dy}{dx}=M(x)$$ 위 식은 y와 x가 분리되어 있는 형태입니다. 이 식의 해를 결정하기 위해 양변을 적분합니다. $$\begin{aligned}&\int N(y)\frac{dy}{dx}\,, dx=\int M(x)\,, dx\\ &\int N(y)\, dy = \int M(x)\, dx \end{aligned}$$ 예) $\frac{dy}{dx}=6y^2x$의 해를 결정합니다. 이 식의 초기 조건은 $y(1)=\frac{1}{25}$입니다. $$\begin{aligned}&\frac{1}{6y^2}\frac{dy}{dx}=x\\ &\int \frac{1}{6y^2}\frac{dy}{dx}\, dx=\int x\, dx \\ &\int \frac{1}{6y^2}\, dy=\int x\, dx\\ &-\frac{1}{6}\frac{1}{y}=\frac{1}{2}x^2+C\\ &\frac{1}{y}=-3x^2+C\end{aligned}$$ 위 식에 초기조건을 대입하여 C를 결정합니다. $$\begin{aligned}&\frac{1}{1/25}=-3+C\\&C=28\end{aligned}$$ 그러므로 해는 다음과 같습니다. $$\begin{aligned}\frac{1}{y}=-3x^2+28 \\ y=\frac{1}{28-3x^2} \end{aligned}$$ x=symbols('x') y=Function('y')(x) eq=Eq(y.diff(x), 6*y**2*x) sol=dsolve(eq, y, ics={y.subs(x, 1): 1/25}); sol $\quad\color{navy}{\scriptstyle y{\left(x \right)} = - \frac{1}{3 x^...