ODE
\[ x y'(x) \left (\text {a0}+\text {b0} x^k\right )+y(x) \left (\text {a1}+\text {b1} x^k+\text {c1} x^{2 k}\right )+x^2 y''(x)=0 \] ODE Classification
[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.280749 (sec), leaf count = 410
\[\left \{\left \{y(x)\to 2^{\frac {1}{2} \left (\frac {\sqrt {k^2 \left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right )}}{k^2}+1\right )} x^{\frac {1}{2} (-\text {a0}-k+1)} \left (x^k\right )^{\frac {1}{2} \left (\frac {\sqrt {k^2 \left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right )}}{k^2}+1\right )} e^{-\frac {\left (\sqrt {\text {b0}^2-4 \text {c1}}+\text {b0}\right ) x^k}{2 k}} \left (c_1 U\left (\frac {\left (k^2+\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}\right ) \text {b0}^2+\sqrt {\text {b0}^2-4 \text {c1}} k (\text {a0}+k-1) \text {b0}-2 \text {b1} \sqrt {\text {b0}^2-4 \text {c1}} k-4 \text {c1} \left (k^2+\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}\right )}{2 \left (\text {b0}^2-4 \text {c1}\right ) k^2},\frac {k^2+\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}}{k^2},\frac {\sqrt {\text {b0}^2-4 \text {c1}} x^k}{k}\right )+c_2 L_{-\frac {\left (k^2+\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}\right ) \text {b0}^2+\sqrt {\text {b0}^2-4 \text {c1}} k (\text {a0}+k-1) \text {b0}-2 \text {b1} \sqrt {\text {b0}^2-4 \text {c1}} k-4 \text {c1} \left (k^2+\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}\right )}{2 \left (\text {b0}^2-4 \text {c1}\right ) k^2}}^{\frac {\sqrt {\left (\text {a0}^2-2 \text {a0}-4 \text {a1}+1\right ) k^2}}{k^2}}\left (\frac {\sqrt {\text {b0}^2-4 \text {c1}} x^k}{k}\right )\right )\right \}\right \}\]
Maple ✓
cpu = 1.076 (sec), leaf count = 167
\[\left [y \left (x \right ) = \textit {\_C1} \,x^{-\frac {\mathit {a0}}{2}+\frac {1}{2}-\frac {k}{2}} {\mathrm e}^{-\frac {\mathit {b0} \,x^{k}}{2 k}} \WhittakerM \left (-\frac {\left (\mathit {a0} -1+k \right ) \mathit {b0} -2 \mathit {b1}}{2 \sqrt {\mathit {b0}^{2}-4 \mathit {c1}}\, k}, \frac {\sqrt {\mathit {a0}^{2}-2 \mathit {a0} -4 \mathit {a1} +1}}{2 k}, \frac {\sqrt {\mathit {b0}^{2}-4 \mathit {c1}}\, x^{k}}{k}\right )+\textit {\_C2} \,x^{-\frac {\mathit {a0}}{2}+\frac {1}{2}-\frac {k}{2}} {\mathrm e}^{-\frac {\mathit {b0} \,x^{k}}{2 k}} \WhittakerW \left (-\frac {\left (\mathit {a0} -1+k \right ) \mathit {b0} -2 \mathit {b1}}{2 \sqrt {\mathit {b0}^{2}-4 \mathit {c1}}\, k}, \frac {\sqrt {\mathit {a0}^{2}-2 \mathit {a0} -4 \mathit {a1} +1}}{2 k}, \frac {\sqrt {\mathit {b0}^{2}-4 \mathit {c1}}\, x^{k}}{k}\right )\right ]\] Mathematica raw input
DSolve[(a1 + b1*x^k + c1*x^(2*k))*y[x] + x*(a0 + b0*x^k)*y'[x] + x^2*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> (2^((1 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]/k^2)/2)*x^((1 - a0 - k)/2)*(x^k)^((1 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]/k^2)/2)*(C[1]*HypergeometricU[(-2*b1*Sqrt[b0^2 - 4*c1]*k + b0*Sqrt[b0^2 - 4*c1]*k*(-1 + a0 + k) + b0^2*(k^2 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]) - 4*c1*(k^2 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]))/(2*(b0^2 - 4*c1)*k^2), (k^2 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2])/k^2, (Sqrt[b0^2 - 4*c1]*x^k)/k] + C[2]*LaguerreL[-1/2*(-2*b1*Sqrt[b0^2 - 4*c1]*k + b0*Sqrt[b0^2 - 4*c1]*k*(-1 + a0 + k) + b0^2*(k^2 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]) - 4*c1*(k^2 + Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]))/((b0^2 - 4*c1)*k^2), Sqrt[(1 - 2*a0 + a0^2 - 4*a1)*k^2]/k^2, (Sqrt[b0^2 - 4*c1]*x^k)/k]))/E^(((b0 + Sqrt[b0^2 - 4*c1])*x^k)/(2*k))}}
Maple raw input
dsolve(x^2*diff(diff(y(x),x),x)+x*(a0+b0*x^k)*diff(y(x),x)+(a1+b1*x^k+c1*x^(2*k))*y(x) = 0, y(x))
Maple raw output
[y(x) = _C1*x^(-1/2*a0+1/2-1/2*k)*exp(-1/2/k*b0*x^k)*WhittakerM(-1/2*((a0-1+k)*b0-2*b1)/(b0^2-4*c1)^(1/2)/k,1/2*(a0^2-2*a0-4*a1+1)^(1/2)/k,(b0^2-4*c1)^(1/2)/k*x^k)+_C2*x^(-1/2*a0+1/2-1/2*k)*exp(-1/2/k*b0*x^k)*WhittakerW(-1/2*((a0-1+k)*b0-2*b1)/(b0^2-4*c1)^(1/2)/k,1/2*(a0^2-2*a0-4*a1+1)^(1/2)/k,(b0^2-4*c1)^(1/2)/k*x^k)]