ODE
\[ (\text {a0}+\text {b0} x) y'(x)+y(x) (\text {a1}+\text {b1} x)+y''(x)=0 \] ODE Classification
[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.0522915 (sec), leaf count = 132
\[\left \{\left \{y(x)\to e^{-\text {a0} x+\frac {\text {b1} x}{\text {b0}}-\frac {\text {b0} x^2}{2}} \left (c_2 \, _1F_1\left (\frac {\text {b0}^3-\text {a1} \text {b0}^2+\text {a0} \text {b1} \text {b0}-\text {b1}^2}{2 \text {b0}^3};\frac {1}{2};\frac {\left (x \text {b0}^2+\text {a0} \text {b0}-2 \text {b1}\right )^2}{2 \text {b0}^3}\right )+c_1 H_{\frac {-\text {a0} \text {b0} \text {b1}+\text {a1} \text {b0}^2-\text {b0}^3+\text {b1}^2}{\text {b0}^3}}\left (\frac {\text {a0} \text {b0}+\text {b0}^2 x-2 \text {b1}}{\sqrt {2} \text {b0}^{3/2}}\right )\right )\right \}\right \}\]
Maple ✓
cpu = 0.074 (sec), leaf count = 98
\[ \left \{ y \left ( x \right ) ={{\rm e}^{-{\frac {{\it b1}\,x}{{\it b0}}}}} \left ( {{\sl U}\left ({\frac {-{\it a0}\,{\it b0}\,{\it b1}+{\it a1}\,{{\it b0}}^{2}+{{\it b1}}^{2}}{2\,{{\it b0}}^{3}}},\,{\frac {1}{2}},\,-{\frac { \left ( x{{\it b0}}^{2}+{\it a0}\,{\it b0}-2\,{\it b1} \right ) ^{2}}{2\,{{\it b0}}^{3}}}\right )}{\it \_C2}+{{\sl M}\left ({\frac {-{\it a0}\,{\it b0}\,{\it b1}+{\it a1}\,{{\it b0}}^{2}+{{\it b1}}^{2}}{2\,{{\it b0}}^{3}}},\,{\frac {1}{2}},\,-{\frac { \left ( x{{\it b0}}^{2}+{\it a0}\,{\it b0}-2\,{\it b1} \right ) ^{2}}{2\,{{\it b0}}^{3}}}\right )}{\it \_C1} \right ) \right \} \] Mathematica raw input
DSolve[(a1 + b1*x)*y[x] + (a0 + b0*x)*y'[x] + y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> E^(-(a0*x) + (b1*x)/b0 - (b0*x^2)/2)*(C[1]*HermiteH[(a1*b0^2 - b0^3 - a0*b0*b1 + b1^2)/b0^3, (a0*b0 - 2*b1 + b0^2*x)/(Sqrt[2]*b0^(3/2))] + C[2]*Hypergeometric1F1[(-(a1*b0^2) + b0^3 + a0*b0*b1 - b1^2)/(2*b0^3), 1/2, (a0*b0 - 2*b1 + b0^2*x)^2/(2*b0^3)])}}
Maple raw input
dsolve(diff(diff(y(x),x),x)+(b0*x+a0)*diff(y(x),x)+(b1*x+a1)*y(x) = 0, y(x),'implicit')
Maple raw output
y(x) = exp(-b1/b0*x)*(KummerU(1/2*(-a0*b0*b1+a1*b0^2+b1^2)/b0^3,1/2,-1/2/b0^3*(b0^2*x+a0*b0-2*b1)^2)*_C2+KummerM(1/2*(-a0*b0*b1+a1*b0^2+b1^2)/b0^3,1/2,-1/2/b0^3*(b0^2*x+a0*b0-2*b1)^2)*_C1)