ODE
\[ y'''(x)=x^2+y(x) \] ODE Classification
[[_3rd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.154347 (sec), leaf count = 59
\[\left \{\left \{y(x)\to -x^2+c_1 e^x+c_2 e^{-x/2} \cos \left (\frac {\sqrt {3} x}{2}\right )+c_3 e^{-x/2} \sin \left (\frac {\sqrt {3} x}{2}\right )\right \}\right \}\]
Maple ✓
cpu = 0.027 (sec), leaf count = 40
\[\left [y \left (x \right ) = -x^{2}+\textit {\_C1} \,{\mathrm e}^{x}+\textit {\_C2} \,{\mathrm e}^{-\frac {x}{2}} \cos \left (\frac {\sqrt {3}\, x}{2}\right )+\textit {\_C3} \,{\mathrm e}^{-\frac {x}{2}} \sin \left (\frac {\sqrt {3}\, x}{2}\right )\right ]\] Mathematica raw input
DSolve[y'''[x] == x^2 + y[x],y[x],x]
Mathematica raw output
{{y[x] -> -x^2 + E^x*C[1] + (C[2]*Cos[(Sqrt[3]*x)/2])/E^(x/2) + (C[3]*Sin[(Sqrt[3]*x)/2])/E^(x/2)}}
Maple raw input
dsolve(diff(diff(diff(y(x),x),x),x) = x^2+y(x), y(x))
Maple raw output
[y(x) = -x^2+_C1*exp(x)+_C2*exp(-1/2*x)*cos(1/2*3^(1/2)*x)+_C3*exp(-1/2*x)*sin(1/2*3^(1/2)*x)]