2.4.1.5 Maple. Time used: 0.001 (sec). Leaf size: 14
ode:=diff(y(x),x)+y(x)*x = x; 
dsolve(ode,y(x), singsol=all);
 
\begin{align*} y = 1+{\mathrm e}^{-\frac {x^{2}}{2}} c_1 \end{align*}

Maple trace

Methods for first order ODEs: 
--- Trying classification methods --- 
trying a quadrature 
trying 1st order linear 
<- 1st order linear successful
 

Maple step by step

\[ \begin {array}{lll} & {} & \textrm {Let's solve}\hspace {3pt} \\ {} & {} & \frac {d}{d x}y \left (x \right )+x y \left (x \right )=x \\ \bullet & {} & \textrm {Highest derivative means the order of the ODE is}\hspace {3pt} 1 \\ {} & {} & \frac {d}{d x}y \left (x \right ) \\ \bullet & {} & \textrm {Solve for the highest derivative}\hspace {3pt} \\ {} & {} & \frac {d}{d x}y \left (x \right )=-x y \left (x \right )+x \\ \bullet & {} & \textrm {Separate variables}\hspace {3pt} \\ {} & {} & \frac {\frac {d}{d x}y \left (x \right )}{y \left (x \right )-1}=-x \\ \bullet & {} & \textrm {Integrate both sides with respect to}\hspace {3pt} x \\ {} & {} & \int \frac {\frac {d}{d x}y \left (x \right )}{y \left (x \right )-1}d x =\int -x d x +\mathit {C1} \\ \bullet & {} & \textrm {Evaluate integral}\hspace {3pt} \\ {} & {} & \ln \left (y \left (x \right )-1\right )=-\frac {x^{2}}{2}+\mathit {C1} \\ \bullet & {} & \textrm {Solve for}\hspace {3pt} y \left (x \right ) \\ {} & {} & y \left (x \right )={\mathrm e}^{-\frac {x^{2}}{2}+\mathit {C1}}+1 \\ \bullet & {} & \textrm {Redefine the integration constant(s)}\hspace {3pt} \\ {} & {} & y \left (x \right )=\mathit {C1} \,{\mathrm e}^{-\frac {x^{2}}{2}}+1 \end {array} \]