2.3   ODE No. 3

  1. Problem in Latex
  2. Mathematica input
  3. Maple input

\[ a y(x)-b \sin (c x)+y'(x)=0 \] Mathematica : cpu = 0.0729384 (sec), leaf count = 40

\[\left \{\left \{y(x)\to \frac {b (a \sin (c x)-c \cos (c x))}{a^2+c^2}+c_1 e^{-a x}\right \}\right \}\] Maple : cpu = 0.035 (sec), leaf count = 38

\[\left \{y \left (x \right ) = c_{1} {\mathrm e}^{-a x}-\frac {\left (-a \sin \left (c x \right )+c \cos \left (c x \right )\right ) b}{a^{2}+c^{2}}\right \}\]

Hand solution

\begin {equation} \frac {dy}{dx}+ay\left ( x\right ) =b\sin \left ( cx\right ) \tag {1} \end {equation}

Integrating factor \(\mu =e^{\int adx}=e^{ax}\). Hence (1) becomes

\begin {align*} \frac {d}{dx}\left ( \mu y\left ( x\right ) \right ) & =\mu b\sin \left ( cx\right ) \\ \mu y\left ( x\right ) & =b\int \mu \sin \left ( cx\right ) dx+C \end {align*}

Replacing \(\mu \) by \(e^{ax}\)

\begin {equation} y\left ( x\right ) =be^{-ax}\int e^{ax}\sin \left ( cx\right ) dx+Ce^{-ax}\tag {2} \end {equation}

Using \(\sin \left ( cx\right ) =\frac {e^{icx}-e^{-icx}}{2i}\) then \begin {align*} \int e^{ax}\sin \left ( cx\right ) dx & =\int \frac {e^{\left ( ic+a\right ) x}-e^{\left ( -ic+a\right ) x}}{2i}dx\\ & =\frac {1}{2i}\left ( \frac {e^{\left ( ic+a\right ) x}}{ic+a}-\frac {e^{\left ( -ic+a\right ) x}}{-ic+a}\right ) \\ & =\frac {1}{2i}e^{ax}\left ( \frac {e^{icx}}{ic+a}-\frac {e^{-icx}}{-ic+a}\right ) \\ & =\frac {1}{2i}e^{ax}\left ( \frac {e^{icx}\left ( -ic+a\right ) -e^{-icx}\left ( ic+a\right ) }{\left ( ic+a\right ) \left ( -ic+a\right ) }\right ) \\ & =\frac {1}{2i}e^{ax}\left ( \frac {-ice^{icx}+ae^{icx}-ice^{-icx}-ae^{-icx}}{\left ( c^{2}+a^{2}\right ) }\right ) \\ & =\frac {1}{2i}e^{ax}\left ( \frac {-ic\left ( e^{icx}+e^{-icx}\right ) +a\left ( e^{icx}-e^{-icx}\right ) }{\left ( c^{2}+a^{2}\right ) }\right ) \\ & =\frac {e^{ax}}{\left ( c^{2}+a^{2}\right ) }\left ( \frac {-ic\left ( e^{icx}+e^{-icx}\right ) }{2i}+\frac {a\left ( e^{icx}-e^{-icx}\right ) }{2i}\right ) \\ & =\frac {e^{ax}}{\left ( c^{2}+a^{2}\right ) }\left ( -c\cos cx+a\sin cx\right ) \end {align*}

Therefore (2) becomes

\begin {align*} y\left ( x\right ) & =be^{-ax}\left [ \frac {e^{ax}}{\left ( c^{2}+a^{2}\right ) }\left ( -c\cos cx+a\sin cx\right ) \right ] +Ce^{-ax}\\ & =\frac {b}{\left ( c^{2}+a^{2}\right ) }\left ( -c\cos cx+a\sin cx\right ) +Ce^{-ax} \end {align*}